Physical Sciences and Mathematics Commons^™

On Two-Stage Hypothesis Testing Procedures Via Asymptotically Independent Statistics, James Y. Dai, Charles Kooperberg, Michael Leblanc, Ross L. Prentice

UW Biostatistics Working Paper Series

Kooperberg08 proposed a two-stage testing procedure to screen for significant interactions in genome-wide association (GWA) studies by a soft threshold on marginal associations (MA), though its theoretical properties and generalization have not been elaborated. In this article, we discuss conditions that are required to achieve strong control of the Family-Wise Error Rate (FWER) by such procedures for low or high-dimensional hypothesis testing. We provide proof of asymptotic independence of marginal association statistics and interaction statistics in linear regression, logistic regression, and Cox proportional hazard models in a randomized clinical trial (RCT) with a rare event. In case-control studies nested within …

Multi-State Life Tables, Equilibrium Prevalence, And Baseline Selection Bias, Paula Diehr, David Yanez

UW Biostatistics Working Paper Series

Consider a 3-state system with one absorbing state, such as Healthy, Sick, and Dead. If the system satisfies the 1-step Markov conditions, the prevalence of the Healthy state will converge to a value that is independent of the initial distribution. This equilibrium prevalence and its variance are known under the assumption of time homogeneity, and provided reasonable estimates in the time non-homogeneous systems studied. Here, we derived the equilibrium prevalence for a system with more than three states. Under time homogeneity, the equilibrium prevalence distribution was shown to be an eigenvector of a partition of the matrix of transition probabilities. …

Model-Robust Regression And A Bayesian `Sandwich' Estimator, Adam A. Szpiro, Kenneth M. Rice, Thomas Lumley

UW Biostatistics Working Paper Series

The published version of this paper in Annals of Applied Statistics (Vol. 4, No. 4 (2010), 2099–2113) is available from the journal web site at http://dx.doi.org/10.1214/10-AOAS362.

We present a new Bayesian approach to model-robust linear regression that leads to uncertainty estimates with the same robustness properties as the Huber-White sandwich estimator. The sandwich estimator is known to provide asymptotically correct frequentist inference, even when standard modeling assumptions such as linearity and homoscedasticity in the data-generating mechanism are violated. Our derivation provides a compelling Bayesian justification for using this simple and popular tool, and it also clarifies what is being estimated …

Asymptotic Properties Of The Sequential Empirical Roc And Ppv Curves, Joseph S. Koopmeiners, Ziding Feng

UW Biostatistics Working Paper Series

The receiver operating characteristic (ROC) curve, the positive predictive value (PPV) curve and the negative predictive value (NPV) curve are three common measures of performance for a diagnostic biomarker. The independent increments covariance structure assumption is common in the group sequential study design literature. Showing that summary measures of the ROC, PPV and NPV curves have an independent increments covariance structure will provide the theoretical foundation for designing group sequential diagnostic biomarker studies. The ROC, PPV and NPV curves are often estimated empirically to avoid assumptions about the distributional form of the biomarkers. In this paper we derive asymptotic theory …

Nonparametric And Semiparametric Analysis Of Current Status Data Subject To Outcome Misclassification, Victor G. Sal Y Rosas, James P. Hughes

UW Biostatistics Working Paper Series

In this article, we present nonparametric and semiparametric methods to analyze current status data subject to outcome misclassification. Our methods use nonparametric maximum likelihood estimation (NPMLE) to estimate the distribution function of the failure time when sensitivity and specificity may vary among subgroups. A nonparametric test is proposed for the two sample hypothesis testing. In regression analysis, we apply the Cox proportional hazard model and likelihood ratio based confidence intervals for the regression coefficients are proposed. Our methods are motivated and demonstrated by data collected from an infectious disease study in Seattle, WA.

Panel Count Data Regression With Informative Observation Times, Petra Buzkova

UW Biostatistics Working Paper Series

When patients are monitored for potentially recurrent events such as infections or tumor metastases, it is common for clinicians to ask patients to come back sooner for follow-up based on the results of the most recent exam. This means that subjects’ observation times will be irregular and related to subject-specific factors. Previously proposed methods for handling such panel count data assume that the dependence between the events process and the observation time process is time-invariant. This article considers situations where the observation times are predicted by time-varying factors, such as the outcome observed at the last visit or cumulative exposure. …

Bio-Creep In Non-Inferiority Clinical Trials, Siobhan P. Everson-Stewart, Scott S. Emerson

UW Biostatistics Working Paper Series

After a non-inferiority clinical trial, a new therapy may be accepted as effective, even if its treatment effect is slightly smaller than the current standard. It is therefore possible that, after a series of trials where the new therapy is slightly worse than the preceding drugs, an ineffective or harmful therapy might be incorrectly declared efficacious; this is known as “bio-creep.” Several factors may influence the rate at which bio-creep occurs, including the distribution of the effects of the new agents being tested and how that changes over time, the choice of active comparator, the method used to model the …

Estimates Of Information Growth In Longitudinal Clinical Trials, Abigail Shoben, Kyle Rudser, Scott S. Emerson

UW Biostatistics Working Paper Series

In group sequential clinical trials, it is necessary to estimate the amount of information present at interim analysis times relative to the amount of information that would be present at the final analysis. If only one measurement is made per individual, this is often the ratio of sample sizes available at the interim and final analyses. However, as discussed by Wu and Lan (1992), when the statistic of interest is a change over time, as with longitudinal data, such an approach overstates the information. In this paper, we discuss other problems that can result in overestimating the information, such as …

Robustness Of Approaches To Roc Curve Modeling Under Misspecification Of The Underlying Probability Model, Sean Devlin, Elizabeth Thomas, Scott S. Emerson

UW Biostatistics Working Paper Series

The receiver operating characteristic (ROC) curve is a tool of particular use in disease status classification with a continuous medical test (marker). A variety of statistical regression models have been proposed for the comparison of ROC curves for different markers across covariate groups. A full parametric modeling of the marker distribution has been generally found to be overly reliant on the strong parametric assumptions. Pepe (2003) has instead developed parametric models for the ROC curve that induce a semi-parametric model for the marker distributions. The estimating equations proposed for use in these ROC-GLM models may differ from commonly used estimating …

Exploring The Benefits Of Adaptive Sequential Designs In Time-To-Event Endpoint Settings, Sarah C. Emerson, Kyle Rudser, Scott S. Emerson

UW Biostatistics Working Paper Series

Sequential analysis is frequently employed to address ethical and financial issues in clinical trials. Sequential analysis may be performed using standard group sequential designs, or, more recently, with adaptive designs that use estimates of treatment effect to modify the maximal statistical information to be collected. In the general setting in which statistical information and clinical trial costs are functions of the number of subjects used, it has yet to be established whether there is any major efficiency advantage to adaptive designs over traditional group sequential designs. In survival analysis, however, statistical information (and hence efficiency) is most closely related to …

Pragmatic Estimation Of A Spatio-Temporal Air Quality Model With Irregular Monitoring Data, Paul D. Sampson, Adam A. Szpiro, Lianne Sheppard, Johan Lindström, Joel D. Kaufman

UW Biostatistics Working Paper Series

Statistical analyses of the health effects of air pollution have increasingly used GIS-based covariates for prediction of ambient air quality in “land-use” regression models. More recently these regression models have accounted for spatial correlation structure in combining monitoring data with land-use covariates. The current paper builds on these concepts to address spatio-temporal prediction of ambient concentrations of particulate matter with aerodynamic diameter less than 2.5 μm (PM2.5) on the basis of a model representing spatially varying seasonal trends and spatial correlation structures. Our hierarchical methodology provides a pragmatic approach that fully exploits regulatory and other supplemental monitoring data which jointly …

Is Survival The Only Or Even The Right Outcome For Evaluating Treatments For Out-Of-Hospital Cardiac Arrest? A Proposed Test Based On Both An Intermediate And Ultimate Outcome., Al Hallstrom

UW Biostatistics Working Paper Series

It is generally agreed that the goal of resuscitation is survival with neurological and physiological status similar to that preceding the cardiac arrest. Previously I have argued that the lack of improvement in outcome from resuscitation over the past 3 to 4 decades, as compared to the substantial progress made in treatment of ischemic heart disease, is a consequence of the absence of randomized clinical trials of new interventions and the use of intermediate endpoints such as return of spontaneous circulation or admittance to hospital. Proponents of these intermediate endpoints have argued that those involved in the resuscitation have no …

Robustness Of Semiparametric Efficiency In Nearly-Correct Models For Two-Phase Samples, Thomas Lumley

UW Biostatistics Working Paper Series

Augmented inverse-probability weighted (AIPW) estimators for incomplete-data models typically do not have full semiparametric efficiency, but do have model-robustness properties not shared by the efficient estimator. We examine the performance of efficient and AIPW estimators when the complete-data model is nearly correctly specified, in the sense that the misspecification is not reliably detectable from the data by any possible diagnostic or test. Asymptotic results for these nearly true models are obtained by representing them as sequences of misspecified models that are mutually contiguous with a correctly specified model. For some least favorable direction of model misspecification the bias in the …

Nonparametric And Semiparametric Estimation Of The Three Way Receiver Operating Characteristic Surface, Jialiang Li, Xiao-Hua Zhou

UW Biostatistics Working Paper Series

In many situations the diagnostic decision is not limited to a binary choice. Binary statistical tools such as receiver operating characteristic (ROC) curve and area under the ROC curve (AUC) need to be expanded to address three-category classification problem. Previous authors have suggest various ways to model the extension of AUC but not the ROC surface. Only simple parametric approaches are proposed for modeling the ROC measure under the assumption that test results all follow normal distributions. We study the estimation methods of three dimensional ROC surfaces with nonparametric and semiparametric estimators. Asymptotical results are provided as a basis for …

Evaluating Markers For Treatment Selection Based On Survival Time, Xiao Song, Xiao-Hua Zhou

UW Biostatistics Working Paper Series

For many medical conditions several treatment options may be available for treating patients. We consider evaluating markers based on a simple treatment selection policy that incorporates information on the patient's marker value exceeding a threshold. For example, colon cancer patients may be treated by surgery alone or surgery plus chemotherapy. The c-myc gene expression level may be used as a biomarker for treatment selection. Although traditional regression methods may assess the effect of the marker and treatment on outcomes, it is appealing to quantify more directly the potential impact on the population of using the marker to select treatment. A …

Interval Estimation For The Difference In Paired Areas Under The Roc Curves In The Absence Of A Gold Standard Test, Hsin-Neng Hsieh, Hsiu-Yuan Su, Xiao-Hua Zhou

UW Biostatistics Working Paper Series

Receiver operating characteristic (ROC) curves can be used to assess the accuracy of tests measured on ordinal or continuous scales. The most commonly used measure for the overall diagnostic accuracy of diagnostic tests is the area under the ROC curve (AUC). A gold standard test on the true disease status is required to estimate the AUC. However, a gold standard test may sometimes be too expensive or infeasible. Therefore, in many medical research studies, the true disease status of the subjects may remain unknown. Under the normality assumption on test results from each disease group of subjects, using the expectation-maximization …

A Semi-Parametric Two-Part Mixed-Effects Heteroscedastic Transformation Model For Correlated Right-Skewed Semi-Continuous Data, Huazhen Lin, Xiao-Hua Zhou

UW Biostatistics Working Paper Series

In longitudinal or hierarchical structure studies, we often encounter a semi-continuous variable that has a certain proportion of a single value and a continuous and skewed distribution among the rest of values. In the paper, we propose a new semi-parametric two-part mixed-effects transformation model to fit correlated skewed semi-continuous data. In our model, we allow the transformation to be non-parametric. Fitting the proposed model faces computational challenges due to intractable numerical integrations. We derive the estimates for the parameter and the transformation function based on an approximate likelihood, which has high order accuracy but less computational burden. We also propose …

Relaxing Latent Ignorability In The Itt Analysis Of Randomized Studies With Missing Data And Noncompliance, L Taylor, Xiao-Hua Zhou

UW Biostatistics Working Paper Series

Abstract: In this paper we consider the problem in causal inference of estimating the local complier average causal effect (CACE) parameter in the setting of a randomized clinical trial with a binary outcome, cross-over noncompliance, and unintentional missing data on the responses. We focus on the development of a moment estimator that relaxes the assumption of latent ignorability and incorporates sensitivity parameters that represent the relationship between potential outcomes and associated potential response indicators. If conclusions are insensitive over a range of logically possible values of the sensitivity parameters, then the number of interpretations of the data is reduced, and …

Multiple Imputation Methods For Treatment Noncompliance And Nonresponse In Randomized Clinical Trials, Leslie Taylor, Xiao-Hua (Andrew) Zhou

UW Biostatistics Working Paper Series

Summary: Randomized clinical trials are a powerful tool for investigating causal treatment effects, but in human trials there are oftentimes problems of noncompliance which standard analyses, such as the intention-to-treat or as-treated analysis, either ignore or incorporate in such a way that the resulting estimand is no longer a causal effect. One alternative to these analyses is the complier average causal effect (CACE) which estimates the average causal treatment effect among a subpopulation that would comply under any treatment assigned. We focus on the setting of a randomized clinical trial with crossover treatment noncompliance (e.g., control subjects could receive the …

Semiparametric Two-Part Models With Proportionality Constraints: Analysis Of The Multi-Ethnic Study Of Atherosclerosis (Mesa), Anna Liu, Richard Kronmal, Xiao-Hua Zhou, Shuangge Ma

UW Biostatistics Working Paper Series

SUMMARY. In this article, we analyze the coronary artery calcium (CAC) score in the Multi-Ethnic Study of Atherosclerosis (MESA), where about half of the CAC scores are zero and the rest are continuously distributed. When the observed data has a mixture distribution, two-part models can be the natural choice. With a two-part model, there are two covariate effects, with one in each part of the model. Determination of whether the two covariate effects are proportional can provide more insights into the process underlying development and progression of CAC. In this study, we model the CAC score using a semiparametric two-part …

Pooled Nucleic Acid Testing To Identify Antiretroviral Treatment Failure During Hiv Infection, Susanne May, Anthony Gamst, Richard Haubrich, Constance Benson, Davey Smith

UW Biostatistics Working Paper Series

Abstract Background: Pooling strategies have been used to reduce the costs of polymerase chain reaction based screening for acute HIV infection in populations where the prevalence of acute infection is low (<1%). Only limited research has been done for conditions where the prevalence of screening positivity is higher (>1%). Methods and Results: We present data on a variety of pooling strategies that incorporate the use of PCR-based quantitative measures to monitor for virologic failure among HIV-infected patients receiving antiretroviral therapy. For a prevalence of virologic failure between 1% and 25%, we demonstrate relative efficiency and accuracy of various strategies. These results could be used to choose the best strategy based on the requirements of individual laboratory …

Measures To Summarize And Compare The Predictive Capacity Of Markers, Wen Gu, Margaret Pepe

UW Biostatistics Working Paper Series

The predictive capacity of a marker in a population can be described using the population distribution of risk (Huang et al., 2007; Pepe et al., 2008a; Stern, 2008). Virtually all standard statistical summaries of predictability and discrimination can be derived from it (Gail and Pfeiffer, 2005). The goal of this paper is to develop methods for making inference about risk prediction markers using summary measures derived from the risk distribution. We describe some new clinically motivated summary measures and give new interpretations to some existing statistical measures. Methods for estimating these summary measures are described along with distribution theory that …

Synthesis Analysis Of Regression Models With A Continuous Outcome, Andrew Zhou, Nan Hu, Guizhou Hu, Martin Root

UW Biostatistics Working Paper Series

Synthesis Analysis of Regression Models with a Continuous Outcome Xiao-Hua Zhou 1,2, Nan Hu 2, Guizhou Hu3, and Martin Root3 1 HSR&D Center of Excellence, VA Puget Sound Health Care System, Seattle, WA 98101. 2 Department of Biostatistics, University of Washington, Seattle, WA 98195. 3 BioSignia, Inc., 1822 East NC Highway 54, Suite 350, Durham, NC 27713 To estimate the multivariate regression model from multiple individual studies, it would be challenging to obtain results if the input from individual studies only provide univariate or incomplete multivariate regression information. Samsa et al [1] proposed a simple method to combine coefficients from …

Predicting Intra-Urban Variation In Air Pollution Concentrations With Complex Spatio-Temporal Interactions, Adam A. Szpiro, Paul D. Sampson, Lianne Sheppard, Thomas Lumley, Sara D. Adar, Joel Kaufman

UW Biostatistics Working Paper Series

We describe a methodology for assigning individual estimates of long-term average air pollution concentrations that accounts for a complex spatio-temporal correlation structure and can accommodate unbalanced observations. This methodology has been developed as part of the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air), a prospective cohort study funded by the U.S. EPA to investigate the relationship between chronic exposure to air pollution and cardiovascular disease. Our hierarchical model decomposes the space-time field into a “mean” that includes dependence on covariates and spatially varying seasonal and long-term trends and a “residual” that accounts for spatially correlated deviations from the …

Trading Bias For Precision: Decision Theory For Intervals And Sets, Kenneth M. Rice, Thomas Lumley, Adam A. Szpiro

UW Biostatistics Working Paper Series

Interval- and set-valued decisions are an essential part of statistical inference. Despite this, the justification behind them is often unclear, leading in practice to a great deal of confusion about exactly what is being presented. In this paper we review and attempt to unify several competing methods of interval-construction, within a formal decision-theoretic framework. The result is a new emphasis on interval-estimation as a distinct goal, and not as an afterthought to point estimation. We also see that representing intervals as trade-offs between measures of precision and bias unifies many existing approaches -- as well as suggesting interpretable criteria to …

Using Longitudinal Data To Estimate The Effect Of Starting To Exercise On The Health Of Sedentary Older Adults, Paula Diehr, Calvin Hirsch

UW Biostatistics Working Paper Series

Background It is difficult to estimate the effect of exercise on future health from observational data because exercising may be both a cause and an effect of health status. Unadjusted analyses suffer from selection bias (healthier persons more likely to exercise), while adjusted analyses may adjust away some of the benefits of exercise.

Objective To obtain a "low-bias" interpretable estimate of the effect of exercise on future health.

Methods We used data from the Cardiovascular Health Study, a longitudinal study of 5,888 older adults. The number of blocks walked in the previous week, collected annually, were classified as Sedentary (less …

Estimation For Arbitrary Functionals Of Survival, Kyle Rudser, Michael L. Leblanc, Scott S. Emerson

UW Biostatistics Working Paper Series

No abstract provided.

Semiparametric And Nonparametric Methods For Evaluating Risk Prediction Markers In Case-Control Studies, Ying Huang, Margaret Pepe

UW Biostatistics Working Paper Series

The performance of a well calibrated risk model, Risk(Y)=P(D=1|Y), can be characterized by the population distribution of Risk(Y) and displayed with the predictiveness curve. Better performance is characterized by a wider distribution of Risk(Y), since this corresponds to better risk stratification in the sense that more subjects are identified at low and high risk for the outcome D=1. Although methods have been developed to estimate predictiveness curves from cohort studies, most studies to evaluate novel risk prediction markers employ case-control designs. Here we develop semiparametric and nonparametric methods that accommodate case-control data and assume apriori knowledge of P(D=1). Large and …

Accounting For Errors From Predicting Exposures In Environmental Epidemiology And Environmental Statistics, Adam A. Szpiro, Lianne Sheppard, Thomas Lumley

UW Biostatistics Working Paper Series

PLEASE NOTE THAT AN UPDATED VERSION OF THIS RESEARCH IS AVAILABLE AS WORKING PAPER 350 IN THE UNIVERSITY OF WASHINGTON BIOSTATISTICS WORKING PAPER SERIES (http://www.bepress.com/uwbiostat/paper350).

In environmental epidemiology and related problems in environmental statistics, it is typically not practical to directly measure the exposure for each subject. Environmental monitoring is employed with a statistical model to assign exposures to individuals. The result is a form of exposure misspecification that can result in complicated errors in the health effect estimates if the exposure is naively treated as known. The exposure error is neither “classical” nor “Berkson”, so standard regression calibration methods …

Semiparametric Methods For Evaluating The Covariate-Specific Predictiveness Of Continuous Markers In Matched Case-Control Studies, Ying Huang, Margaret S. Pepe

UW Biostatistics Working Paper Series

To assess the value of a continuous marker in predicting the risk of a disease, a graphical tool called the predictiveness curve has been proposed. It characterizes the marker's predictiveness, or capacity to risk stratify the population by displaying the population distribution of risk endowed by the marker. Methods for making inference about the curve and for comparing curves in a general population have been developed. However, knowledge about a marker's performance in the general population only is not enough. Since a marker's effect on the risk model and its distribution can both differ across subpopulations, its predictiveness may vary …