Physical Sciences and Mathematics Commons^™

Population Intervention Models In Causal Inference, Alan E. Hubbard, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

Marginal structural models (MSM) provide a powerful tool for estimating the causal effect of a] treatment variable or risk variable on the distribution of a disease in a population. These models, as originally introduced by Robins (e.g., Robins (2000a), Robins (2000b), van der Laan and Robins (2002)), model the marginal distributions of treatment-specific counterfactual outcomes, possibly conditional on a subset of the baseline covariates, and its dependence on treatment. Marginal structural models are particularly useful in the context of longitudinal data structures, in which each subject's treatment and covariate history are measured over time, and an outcome is recorded at …

Gauss-Seidel Estimation Of Generalized Linear Mixed Models With Application To Poisson Modeling Of Spatially Varying Disease Rates, Subharup Guha, Louise Ryan

Harvard University Biostatistics Working Paper Series

Generalized linear mixed models (GLMMs) provide an elegant framework for the analysis of correlated data. Due to the non-closed form of the likelihood, GLMMs are often fit by computational procedures like penalized quasi-likelihood (PQL). Special cases of these models are generalized linear models (GLMs), which are often fit using algorithms like iterative weighted least squares (IWLS). High computational costs and memory space constraints often make it difficult to apply these iterative procedures to data sets with very large number of cases.

This paper proposes a computationally efficient strategy based on the Gauss-Seidel algorithm that iteratively fits sub-models of the GLMM …

Computational Techniques For Spatial Logistic Regression With Large Datasets, Christopher J. Paciorek, Louise Ryan

Harvard University Biostatistics Working Paper Series

In epidemiological work, outcomes are frequently non-normal, sample sizes may be large, and effects are often small. To relate health outcomes to geographic risk factors, fast and powerful methods for fitting spatial models, particularly for non-normal data, are required. We focus on binary outcomes, with the risk surface a smooth function of space. We compare penalized likelihood models, including the penalized quasi-likelihood (PQL) approach, and Bayesian models based on fit, speed, and ease of implementation.

A Bayesian model using a spectral basis representation of the spatial surface provides the best tradeoff of sensitivity and specificity in simulations, detecting real spatial …

Is The Number Of Sick Persons In A Cohort Constant Over Time?, Paula Diehr, Ann Derleth, Anne Newman, Liming Cai

UW Biostatistics Working Paper Series

Objectives: To estimate the number of persons in a cohort who are sick, over time.

Methods: We calculated the number of sick persons in the Cardiovascular Health Study (CHS), a cohort study of older adults followed up to 14 years, using eight definitions of “healthy” and “sick”. We projected the number in each health state over time for a birth cohort.

Results: The number of sick persons in CHS was approximately constant for 14 years, for all definitions of “sick”. The estimated number of sick persons in the birth cohort was approximately constant from ages 55-75, after which it decreased. …

A Nonstationary Negative Binomial Time Series With Time-Dependent Covariates: Enterococcus Counts In Boston Harbor, E. Andres Houseman, Brent Coull, James P. Shine

Harvard University Biostatistics Working Paper Series

Boston Harbor has had a history of poor water quality, including contamination by enteric pathogens. We conduct a statistical analysis of data collected by the Massachusetts Water Resources Authority (MWRA) between 1996 and 2002 to evaluate the effects of court-mandated improvements in sewage treatment. Motivated by the ineffectiveness of standard Poisson mixture models and their zero-inflated counterparts, we propose a new negative binomial model for time series of Enterococcus counts in Boston Harbor, where nonstationarity and autocorrelation are modeled using a nonparametric smooth function of time in the predictor. Without further restrictions, this function is not identifiable in the presence …

Direct Effect Models, Mark J. Van Der Laan, Maya L. Petersen

U.C. Berkeley Division of Biostatistics Working Paper Series

The causal effect of a treatment on an outcome is generally mediated by several intermediate variables. Estimation of the component of the causal effect of a treatment that is mediated by a given intermediate variable (the indirect effect of the treatment), and the component that is not mediated by that intermediate variable (the direct effect of the treatment) is often relevant to mechanistic understanding and to the design of clinical and public health interventions. Under the assumption of no-unmeasured confounders for treatment and the intermediate variable, Robins & Greenland (1992) define an individual direct effect as the counterfactual effect of …

Development Of Dose Conversion Coefficients For Radionuclides Produced In Spallation Neutron Sources Quarterly Progress Report 4/01/05 – 6/30/05, Phillip W. Patton, Mark Rudin

Transmutation Sciences Physics (TRP)

The research consortium comprised of representatives from several universities and national laboratories has successfully generated internal and external dose conversion coefficients for twenty radionuclides produced in spallation neutron sources. These dose coefficients fill data gaps exist in Federal Guide Report No. 11 and in Publications 68 and 72 of the International Commission on Radiological Protection (ICRP), and two articles containing the data have been accepted for publication in the Journal of Health Physics. Currently, more nuclear data is needed for the rare radionuclides produced from a mercury target. While attempting to develop a workable plan to acquire this missing data, …

Attributable Risk Function In The Proportional Hazards Model, Ying Qing Chen, Chengcheng Hu, Yan Wang

UW Biostatistics Working Paper Series

As an epidemiological parameter, the population attributable fraction is an important measure to quantify the public health attributable risk of an exposure to morbidity and mortality. In this article, we extend this parameter to the attributable fraction function in survival analysis of time-to-event outcomes, and further establish its estimation and inference procedures based on the widely used proportional hazards models. Numerical examples and simulations studies are presented to validate and demonstrate the proposed methods.

Relationships Between Environmental Factors And Fungi On Occupants' Perceptions Of Indoor Air Quality, Monica Rodriguez

Masters Theses & Specialist Projects

In recent decades, concerns about potential health effects resulting from exposure to contaminants that cause indoor air pollution have dramatically increased. The purpose of this study was to assess the indoor air quality of three buildings at Western Kentucky University and to examine and characterize indoor levels of basic comfort parameters, carbon dioxide, and fungi as well as occupants' perceptions of poor indoor air quality and the role of fungi on reported health symptoms. The three buildings included in the study were: Tate Page Hall (TPH), Jones Jagger Hall (JJH) and Science and Technology Hall (STH). Fifty-three questionnaires were completed …

Exposure Assessment Of Mercury In Fish From Kentucky Surface Waters, Najla Gubari

Masters Theses & Specialist Projects

Consumption of fish contaminated with mercury is the primary exposure pathway by which humans are exposed to mercury. Mercury is known to be a neurological toxin that can cause attention and language deficits, impaired memory, and impaired visual and motor function, especially in children under the age of six. In addition, mercury exposure can lead to other health problems in adults, such as damage to the kidney and increased risk of coronary heart disease. This study assesses the exposure to mercury through consumption of fish from Kentucky surface waters. The exposure assessment is based on data collected by the Kentucky …

Causal Inference In Longitudinal Studies With History-Restricted Marginal Structural Models, Romain Neugebauer, Mark J. Van Der Laan, Ira B. Tager

U.C. Berkeley Division of Biostatistics Working Paper Series

Causal Inference based on Marginal Structural Models (MSMs) is particularly attractive to subject-matter investigators because MSM parameters provide explicit representations of causal effects. We introduce History-Restricted Marginal Structural Models (HRMSMs) for longitudinal data for the purpose of defining causal parameters which may often be better suited for Public Health research. This new class of MSMs allows investigators to analyze the causal effect of a treatment on an outcome based on a fixed, shorter and user-specified history of exposure compared to MSMs. By default, the latter represents the treatment causal effect of interest based on a treatment history defined by the …

Development Of Dose Conversion Coefficients For Radionuclides Produced In Spallation Neutron Sources Quarterly Progress Report 1/1/05 – 3/31/05, Phillip W. Patton, Mark Rudin

Transmutation Sciences Physics (TRP)

The research consortium comprised of representatives from several universities and national laboratories has successfully generated internal and external dose conversion coefficients for twenty radionuclides produced in spallation neutron sources. In addition, the group has identified twenty radionuclide that are missing electron capture files and eighteen additional radionuclides missing substantial physical data. The goal for the current year is to develop a methodology that will allow for producing the dose coefficients for these radionuclides with missing data. Methods to obtain these data are being investigated.

Against 'Individual Risk': A Sympathetic Critique Of Risk Assessment, Matthew D. Adler

All Faculty Scholarship

"Individual risk" currently plays a major role in risk assessment and in the regulatory practices of the health and safety agencies that employ risk assessment, such as EPA, FDA, OSHA, NRC, CPSC, and others. Risk assessors use the term "population risk" to mean the number of deaths caused by some hazard. By contrast, "individual risk" is the incremental probability of death that the hazard imposes on some particular person. Regulatory decision procedures keyed to individual risk are widespread. This is true both for the regulation of toxic chemicals (the heartland of risk assessment), and for other health hazards, such as …

Insights Into Latent Class Analysis, Margaret S. Pepe, Holly Janes

UW Biostatistics Working Paper Series

Latent class analysis is a popular statistical technique for estimating disease prevalence and test sensitivity and specificity. It is used when a gold standard assessment of disease is not available but results of multiple imperfect tests are. We derive analytic expressions for the parameter estimates in terms of the raw data, under the conditional independence assumption. These expressions indicate explicitly how observed two- and three-way associations between test results are used to infer disease prevalence and test operating characteristics. Although reasonable if the conditional independence model holds, the estimators have no basis when it fails. We therefore caution against using …

Standardizing Markers To Evaluate And Compare Their Performances, Margaret S. Pepe, Gary M. Longton

UW Biostatistics Working Paper Series

Introduction: Markers that purport to distinguish subjects with a condition from those without a condition must be evaluated rigorously for their classification accuracy. A single approach to statistically evaluating and comparing markers is not yet established.

Methods: We suggest a standardization that uses the marker distribution in unaffected subjects as a reference. For an affected subject with marker value Y, the standardized placement value is the proportion of unaffected subjects with marker values that exceed Y.

Results: We apply the standardization to two illustrative datasets. In patients with pancreatic cancer placement values calculated for the CA 19-9 marker are smaller …

Combining Predictors For Classification Using The Area Under The Roc Curve, Margaret S. Pepe, Tianxi Cai, Zheng Zhang, Gary M. Longton

UW Biostatistics Working Paper Series

No single biomarker for cancer is considered adequately sensitive and specific for cancer screening. It is expected that the results of multiple markers will need to be combined in order to yield adequately accurate classification. Typically the objective function that is optimized for combining markers is the likelihood function. In this paper we consider an alternative objective function -- the area under the empirical receiver operating characteristic curve (AUC). We note that it yields consistent estimates of parameters in a generalized linear model for the risk score but does not require specifying the link function. Like logistic regression it yields …

Impact Of The Introduction Of Pneumococcal Conjugate Vaccine On Immunization Coverage Among Infants, Nancy D. Lin, Ken Kleinman, K Arnold Chan, Xian-Jie Yu, Eric K. France, Stanley Xu, Feifei Wei, John Mullooly, Jeanne Santoli, Tracey A. Lieu, Vaccine Safety Datalink Group

Public Health Department Faculty Publication Series

Background

The introduction of pneumococcal conjugate vaccine (PCV) to the U.S. recommended childhood immunization schedule in the year 2000 added three injections to the number of vaccinations a child is expected to receive during the first year of life. Surveys have suggested that the addition of PCV has led some immunization providers to move other routine childhood vaccinations to later ages, which could increase the possibility of missing these vaccines. The purpose of this study was to evaluate whether introduction of PCV affected immunization coverage for recommended childhood vaccinations among 13-month olds in four large provider groups.

Methods

In this …

Development Of Dose Coefficients For Radionuclides Produced In Spallation Targets, Phillip W. Patton, Mark Rudin

Transmutation Sciences Physics (TRP)

Dose coefficients permit simple determination of radiation dose associated with various exposure scenarios, and ultimately permit radiation safety personnel to assess the health risks to workers in a nuclear facility. Specifically, radiation safety personnel use dose coefficients to determine the radiation dose incurred to a tissue or organ system from a given exposure. These parameters are often expressed in terms of Annual Limits on Intake (ALIs) and Derived Air Concentrations (DACs).

The research consortium comprised of representatives from several universities and national laboratories has successfully generated internal and external dose conversion coefficients for twenty radionuclides produced in spallation neutron sources. …

Population And Sustainability, Geoffrey Mcnicoll

Poverty, Gender, and Youth

Sustainability refers to the preservation of human-valued natural capital—the resources that provide environmental services—at a level sufficient to assure the well-being of future generations. This Population Council working paper assess the effects on the total and per capita availability of those services. The degradation of environmental services—exemplified by the overuse of aquifers or (at a global level) of the atmospheric carbon sink—is a significant threat to sustainable development, one that is often exacerbated by population growth. The critical management issue in such cases is the design of effective governing institutions to restrain service demand and safeguard supply. Uncertainties arising from …

Group Versus Individual Academic Detailing To Improve The Use Of Antihypertensive Medications In Primary Care: A Cluster-Randomized Controlled Trial, Steven R. Simon, Sumit R. Majumdar, Lisa A. Prosser, Susanne Salem-Schatz, Cheryl Warner, Ken Kleinman, Irina Miroshnik, Stephen B. Soumerai

Public Health Department Faculty Publication Series

Purpose

To compare group versus individual academic detailing to increase diuretic or β-blocker use in hypertension.

Methods

We conducted a cluster-randomized controlled trial in a large health maintenance organization. Subjects (N=9820) were patients with newly treated hypertension in the year preceding the intervention (N=3692), the 9 months following the intervention (N=3556), and the second year following intervention (N=2572). We randomly allocated 3 practice sites to group detailing (N=227 prescribers), 3 to individual detailing (N=235 prescribers), and 3 to usual care (N=319 prescribers). Individual detailing entailed a physician-educator meeting individually with clinicians to address barriers to prescribing guideline-recommended medications. The group …