Search PubMed⌕ Search

Biomedical subjects

Andrew Gelman

Publications and source records attributed to Andrew Gelman.

8 recordsLinked to original sources

Multiple imputation for model checking: completed-data plots with missing and latent data.

In problems with missing or latent data, a standard approach is to first impute the unobserved data, then perform all statistical analyses on the completed dataset--corresponding to the observed data and imputed unobserved data--using standard procedures for complete-data inference. Here, we extend this approach to model checking by demonstrating the advantages of the use of completed-data model diagnostics on imputed completed datasets. The approach is set in the theoretical framework of Bayesian posterior predictive checks (but, as with missing-data imputation, our methods of missing-data model checking can also be interpreted as "predictive inference" in a non-Bayesian context). We consider the graphical diagnostics within this framework. Advantages of the completed-data approach include: (1) One can often check model fit in terms of quantities that are of key substantive interest in a natural way, which is not always possible using observed data alone. (2) In problems with missing data, checks may be devised that do not require to model the missingness or inclusion mechanism; the latter is useful for the analysis of ignorable but unknown data collection mechanisms, such as are often assumed in the analysis of sample surveys and observational studies. (3) In many problems with latent data, it is possible to check qualitative features of the model (for example, independence of two variables) that can be naturally formalized with the help of the latent data. We illustrate with several applied examples.

Animals↗

Extension of the isobolographic approach to interactions studies between more than two drugs: illustration with the convulsant interaction between pefloxacin, norfloxacin, and theophylline in rats.

This work proposes a model to characterize the additivity or the nonadditivity of combinations of more than two agents. Using a Bayesian framework, we modeled the variability between experimental subjects, the errors that occurred during data collection, and the relationship between effects and concentrations of agents at the effect site. The model was used to characterize the additivity (or non-additivity) of norfloxacin, pefloxacin, and theophylline in causing maximal seizures in male Sprague Dawley rats. Animals received the drugs separately or in various combinations. Drug infusion was stopped at the onset of maximal seizures, and cerebrospinal fluid samples were collected for determination of drug concentration by high-performance liquid chromatography. The model was fitted to concentration data using Markov Chain Monte Carlo techniques. Results showed that induction of seizures by mixtures of theophylline and pefloxacin were additive. Seizure induction by mixtures of norfloxacin and pefloxacin or norfloxacin and theophylline were not additive and, given the model, these drugs interacted negatively. There was no triple interaction effect between the drugs. This study demonstrates the ease with which mixtures of more than two drugs can be analyzed with the proposed model.

Animals↗

Bayesian analysis of serial dilution assays.

In a serial dilution assay, the concentration of a compound is estimated by combining measurements of several different dilutions of an unknown sample. The relation between concentration and measurement is nonlinear and heteroscedastic, and so it is not appropriate to weight these measurements equally. In the standard existing approach for analysis of these data, a large proportion of the measurements are discarded as being above or below detection limits. We present a Bayesian method for jointly estimating the calibration curve and the unknown concentrations using all the data. Compared to the existing method, our estimates have much lower standard errors and give estimates even when all the measurements are outside the "detection limits." We evaluate our method empirically using laboratory data on cockroach allergens measured in house dust samples. Our estimates are much more accurate than those obtained using the usual approach. In addition, we develop a method for determining the "effective weight" attached to each measurement, based on a local linearization of the estimated model. The effective weight can give insight into the information conveyed by each data point and suggests potential improvements in design of serial dilution experiments.

Allergens↗

Using image and curve registration for measuring the goodness of fit of spatial and temporal predictions.

Conventional measures of model fit for indexed data (e.g., time series or spatial data) summarize errors in y, for instance by integrating (or summing) the squared difference between predicted and measured values over a range of x. We propose an approach which recognizes that errors can occur in the x-direction as well. Instead of just measuring the difference between the predictions and observations at each site (or time), we first "deform" the predictions, stretching or compressing along the x-direction or directions, so as to improve the agreement between the observations and the deformed predictions. Error is then summarized by (a) the amount of deformation in x, and (b) the remaining difference in y between the data and the deformed predictions (i.e., the residual error in y after the deformation). A parameter, lambda, controls the tradeoff between (a) and (b), so that as lambda-->infinity no deformation is allowed, whereas for lambda=0 the deformation minimizes the errors in y. In some applications, the deformation itself is of interest because it characterizes the (temporal or spatial) structure of the errors. The optimal deformation can be computed by solving a system of nonlinear partial differential equations, or, for a unidimensional index, by using a dynamic programming algorithm. We illustrate the procedure with examples from nonlinear time series and fluid dynamics.

Algorithms↗

Direct data manipulation for local decision analysis as applied to the problem of arsenic in drinking water from tube wells in Bangladesh.

A wide variety of tools are available, both parametric and nonparametric, for analyzing spatial data. However, it is not always clear how to translate statistical inferences into decision recommendations. This article explores the possibilities of estimating the effects of decision options using very direct manipulation of data, bypassing formal statistical analysis. We illustrate with the application that motivated this research, a study of arsenic in drinking water in nearly 5,000 wells in a small area in rural Bangladesh. We estimate the potential benefits of two possible remedial actions: (1) recommendations that people switch to nearby wells with lower arsenic levels; and (2) drilling new community wells. We use simple nonparametric clustering methods and estimate uncertainties using cross-validation.

Algorithms↗

Severe preservation injury induces Il-6/STAT3 activation with lack of cell cycle progression after partial liver graft transplantation.

Partial liver graft transplantation is a surgical advance developed to overcome severe donor shortage. Survival of these grafts involves recovery from cold ischemia and reperfusion (CIR) injury, immediate regeneration and maintenance of function. Here we examined the outcome of partial liver grafts in comparison to whole grafts following CIR injury. Lewis rats subjected to orthotopic liver transplantation (OLT) with whole grafts preserved in Viaspan were compared to rats receiving 50% and 30% grafts. Outcome was analyzed by survival and regeneration. Transplantation was associated with 100% survival for all grafts, whereas 16 h preservation resulted in 100%, 20% and 0% survival in animals receiving whole, 50% and 30% grafts, respectively. CIR induced increased IL-6 levels in 50% and 30% grafts, and activation of STAT3. Cell cycle progression (cyclin D1) and regeneration (BrdU) was initiated in all livers preserved for 1 or 8 h, but not in partial grafts preserved for 16 h. In conclusion, partial grafts recover from CIR injury through similar molecular pathways to whole grafts. Partial grafts with severe injury fail to achieve cellular proliferation despite the early initiating signals. This failure could be attributed to the impaired ability of the parenchyma to respond to initiating signals for regeneration.

Animals↗

Assessing the implementation and effects of a trauma-focused intervention for youths in residential treatment.

This paper describes methods being used to implement and assess the effects of a trauma-focused intervention in residential treatment programs for youths with emotional and behavioral problems, and histories of maltreatment and exposure to family or community violence. Preliminary baseline profiles of the therapeutic environments and youths are also presented. The intervention, referred to as the Sanctuary Model (Bloom, 1997), is based in social psychiatry, trauma theories, therapeutic community philosophy, and cognitive-behavioral approaches. Within the context of safe, supportive, stable, and socially responsible therapeutic communities, a trauma recovery treatment framework is used to teach youths effective adaptation and coping skills to replace nonadaptive cognitive, social, and behavioral strategies that may have emerged earlier as means of coping with traumatic life experiences.

Adolescent↗

Promotion of well-switching to mitigate the current arsenic crisis in Bangladesh.

OBJECTIVE: To survey tube wells and households in Araihazar upazila, Bangladesh, to set the stage for a long-term epidemiological study of the consequences of chronic arsenic exposure. METHODS: Water samples and household data were collected over a period of 4 months in 2000 from 4997 contiguous tube wells serving a population of 55000, the position of each well being determined to within +/- 30 m using Global Positioning System receivers. Arsenic concentrations were determined by graphite-furnace atomic-absorption spectrometry. In addition, groundwater samples collected every 2 weeks for an entire year from six tube wells were analysed for arsenic by high-resolution inductively coupled plasma-mass spectrometry. FINDINGS: Half of the wells surveyed in Araihazar had been installed in the previous 5 years; 94% were privately owned. Only about 48% of the surveyed wells supplied water with an arsenic content below 50 micro g/l, the current Bangladesh standard for drinking-water. Similar to other regions of Bangladesh and West Bengal, India, the distribution of arsenic in Araihazar is spatially highly variable (range: 5-860 micro g/l) and therefore difficult to predict. Because of this variability, however, close to 90% of the inhabitants live within 100 m of a safe well. Monitoring of six tube wells currently meeting the 50 micro g/l standard showed no indication of a seasonal cycle in arsenic concentrations coupled to the hydrological cycle. This suggests that well-switching is a viable option in Araihazar, at least for the short term. CONCLUSIONS: Well-switching should be more systematically encouraged in Araihazar and many other parts of Bangladesh and West Bengal, India. Social barriers to well-switching need to be better understood and, if possible, overcome.

Arsenic↗