Digital Commons Network^™

Asymptotically Optimal Model Selection Method With Right Censored Outcomes, Sunduz Keles, Mark J. Van Der Laan, Sandrine Dudoit

U.C. Berkeley Division of Biostatistics Working Paper Series

Over the last two decades, non-parametric and semi-parametric approaches that adapt well known techniques such as regression methods to the analysis of right censored data, e.g. right censored survival data, became popular in the statistics literature. However, the problem of choosing the best model (predictor) among a set of proposed models (predictors) in the right censored data setting have not gained much attention. In this paper, we develop a new cross-validation based model selection method to select among predictors of right censored outcomes such as survival times. The proposed method considers the risk of a given predictor based on the …

Tree-Based Multivariate Regression And Density Estimation With Right-Censored Data , Annette M. Molinaro, Sandrine Dudoit, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

We propose a unified strategy for estimator construction, selection, and performance assessment in the presence of censoring. This approach is entirely driven by the choice of a loss function for the full (uncensored) data structure and can be stated in terms of the following three main steps. (1) Define the parameter of interest as the minimizer of the expected loss, or risk, for a full data loss function chosen to represent the desired measure of performance. Map the full data loss function into an observed (censored) data loss function having the same expected value and leading to an efficient estimator …

Ranked Set Sampling Based On Binary Water Quality Data With Covariates, Paul Kvam

Department of Math & Statistics Faculty Publications

A ranked set sample (RSS) is composed of independent order statistics, formed by collecting and ordering independent subsamples, then measuring only one item from each subsample. If the cost of sampling is dominated by data measurement rather than collection or ranking, the RSS technique is known to be superior to ordinary sampling. Experiments based on binary data are not designed to exploit the advantages of ranked set sampling because categorical data typically are as easily measured as ranked, making RSS methods impractical. However, in some environmental and biological studies, the success probability of a bivariate outcome is related to one …

The Psychological Characteristics Of Obese Children, Thomas Taylor Mitchell

Loma Linda University Electronic Theses, Dissertations & Projects

Several studies along with the surgeon general of the United States have identified obesity as a national health issue. Research indicates that adult obesity treatment has relatively poor long-term effects. Consequently increased attention has been given to addressing obesity in children. This focus is particularly vital due to the fact that children are one of the populations showing the most rapid increase in overweight and obesity. The purpose of this study was to identify the some of the unique psychological characteristics of obese children by analyzing the intake data of children and adolescents attending the Loma Linda University Growing Fit …

The Radiopacity Of 29 Resin Restorative Materials Compared To Enamel, Peter Lubisich Iv

Loma Linda University Electronic Theses, Dissertations & Projects

The radiodensity of a composite resin can greatly influence a clinician's ability to analyze a restoration for failure. The purpose of this study was to determine the relative radiopacities of 29 resin restorative materials compared to enamel, and to verify the comparable radiopacity between shades of 23 of the restorative composite resin materials. Composite resin discs measuring >5mm in diameter and 2 mm in thickness were made and polymerized. Radiographs were taken of the composite resin discs, which included an aluminum step wedge for calibration and a densitometer was used to measure the radiodensity of the composite resin samples. The …

Quantification Of Bone Deposition In Onlay Bone Grafting, Gregory Alan Kammeyer

Loma Linda University Electronic Theses, Dissertations & Projects

The purpose of this study was to quantitatively evaluate; 1) mandibular posterior mono-cortical block onlay bone graft success rates 2) the amount of grafted bone retained after initial resorption in radiographic projections. Seven consecutive patients with 14 mandibular posterior implant sites requiring bone augmentation for implants of a standard 3.8mm x 10mm dimension were sought. Mono-cortical block bone grafts from the ipsilateral ascending ramus were placed then observed with two standardized CT scans, dividing them into two groups: Group A had the first CT at 4 weeks after graft placement and a second taken before Prosthodontic reconstruction, 36.7 weeks later. …

Cross-Calibration Of Stroke Disability Measures: Bayesian Analysis Of Longitudinal Ordinal Categorical Data Using Negative Dependence, Giovanni Parmigiani, Heidi W. Ashih, Gregory P. Samsa, Pamela W. Duncan, Sue Min Lai, David B. Matchar

Johns Hopkins University, Dept. of Biostatistics Working Papers

It is common to assess disability of stroke patients using standardized scales, such as the Rankin Stroke Outcome Scale (RS) and the Barthel Index (BI). The Rankin Scale, which was designed for applications to stroke, is based on assessing directly the global conditions of a patient. The Barthel Index, which was designed for general applications, is based on a series of questions about the patient’s ability to carry out 10 basis activities of daily living. As both scales are commonly used, but few studies use both, translating between scales is important in gaining an overall understanding of the efficacy of …

Randomization Tests For Small Samples: An Application For Genetic Expression Data, Gary L. Gadbury, Grier P. Page, Moonseong Heo, John D. Mountz, David B. Allison

Mathematics and Statistics Faculty Research & Creative Works

An advantage of randomization tests for small samples is that an exact P-value can be computed under an additive model. a disadvantage with very small sample sizes is that the resulting discrete distribution for P-values can make it mathematically impossible for a P-value to attain a particular degree of significance. We investigate a distribution of P-values that arises when several thousand randomization tests are conducted simultaneously using small samples, a situation that arises with microarray gene expression data. We show that the distribution yields valuable information regarding groups of genes that are differentially expressed between two groups: A treatment group …

Modern Statistical Methods For Handling Missing Repeated Measurements In Obesity Trial Data: Beyond Locf, Gary L. Gadbury, C. S. Coffey, D. B. Allison

Mathematics and Statistics Faculty Research & Creative Works

This paper brings together some modern statistical methods to address the problem of missing data in obesity trials with repeated measurements. Such missing data occur when subjects miss one or more follow-up visits or drop out early from an obesity trial. a common approach to dealing with missing data because of dropout is 'last observation carried forward' (LOCF). This method, although intuitively appealing, requires restrictive assumptions to produce valid statistical conclusions. We review the need for obesity trials, the assumptions that must be made regarding missing data in such trials, and some modern statistical methods for analyzing data containing missing …

An Extended General Location Model For Causal Inference From Data Subject To Noncompliance And Missing Values, Yahong Peng, Rod Little, Trivellore E. Raghuanthan

The University of Michigan Department of Biostatistics Working Paper Series

Noncompliance is a common problem in experiments involving randomized assignment of treatments, and standard analyses based on intention-to treat or treatment received have limitations. An attractive alternative is to estimate the Complier-Average Causal Effect (CACE), which is the average treatment effect for the subpopulation of subjects who would comply under either treatment (Angrist, Imbens and Rubin, 1996, henceforth AIR). We propose an Extended General Location Model to estimate the CACE from data with non-compliance and missing data in the outcome and in baseline covariates. Models for both continuous and categorical outcomes and ignorable and latent ignorable (Frangakis and Rubin, 1999) …

Inference For The Population Total From Probability-Proportional-To-Size Samples Based On Predictions From A Penalized Spline Nonparametric Model, Hui Zheng, Rod Little

The University of Michigan Department of Biostatistics Working Paper Series

Inference about the finite population total from probability-proportional-to-size (PPS) samples is considered. In previous work (Zheng and Little, 2003), penalized spline (p-spline) nonparametric model-based estimators were shown to generally outperform the Horvitz-Thompson (HT) and generalized regression (GR) estimators in terms of the root mean squared error. In this article we develop model-based, jackknife and balanced repeated replicate variance estimation methods for the p-spline based estimators. Asymptotic properties of the jackknife method are discussed. Simulations show that p-spline point estimators and their jackknife standard errors lead to inferences that are superior to HT or GR based inferences. This suggests that nonparametric …

On The Formation Of Weighting Adjustment Cells For Unit Nonresponse, Sonya Vartivarian, Rod Little

The University of Michigan Department of Biostatistics Working Paper Series

A method is proposed for weighting adjustments for unit nonresponse based on a crossclassification by the estimated propensity to respond and by the predicted mean of a survey outcome. Simulations to assess the performance of the method are described.

Spearman Rank Regression, Jason C. Parcon

Dissertations

The main purpose of this dissertation is to obtain an estimate of the slope parameter in a regression model that is robust to outlying values in both the x - and Y-spaces. The least squares method, though known to be optimal for normal errors, can yield estimates with infinitely large MSE's if the error distribution is thick-tailed. Regular rank-based methods like the Wilcoxon method are known to be robust to outlying values in the Y -space, but it is still grossly affected by outlying values in x -space.

This dissertation derives an estimate of the slope from an estimating …

Locally Efficient Estimation Of Nonparametric Causal Effects On Mean Outcomes In Longitudinal Studies, Romain Neugebauer, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

Marginal Structural Models (MSM) have been introduced by Robins (1998a) as a powerful tool for causal inference as they directly model causal curves of interest, i.e. mean treatment-specific outcomes possibly adjusted for baseline covariates. Two estimators of the corresponding MSM parameters of interest have been proposed, see van der Laan and Robins (2002): the Inverse Probability of Treatment Weighted (IPTW) and the Double Robust (DR) estimators. A parametric MSM approach to causal inference has been favored since the introduction of MSM. It relies on correct specification of a parametric MSM to consistently estimate the parameter of interest using the IPTW …

Adjusting For Non-Ignorable Verification Bias In Clinical Studies For Alzheimer’S Disease, Xiao-Hua Zhou, Pete Castelluccio

UW Biostatistics Working Paper Series

A common problem for comparing the relative accuracy of two screening tests for Alzheimer’s disease (D) in a two-stage design study is verification bias. If the verification bias can be assumed to be ignorable, Zhou and Higgs (2000) have proposed a maximum likelihood approach to compare the relative accuracy of screening tests in a two-stage design study. However, if the verification mechanism also depends on the unobserved disease status, the ignorable assumption does not hold. In this paper, we discuss how to use a profile likelihood approach to compare the relative accuracy of two screening tests for AD without assuming …

Behavioral Evaluation Of The Psychological Welfare And Environmental Requirements Of Agricultural Research Animals: Theory, Measurement, Ethics, And Practical Implications, Lesley A. King

Experimentation Collection

The welfare of agricultural research animals relies not only on measures of good health but also on the presence of positive emotional states and the absence of aversive or unpleasant subjective states such as fear, frustration, or association with pain. Although subjective states are not inherently observable, their interaction with motivational states can be measured through assessment of motivated behavior, which indicates the priority animals place on obtaining or avoiding specific environmental stimuli and thus allows conclusions regarding the impact of housing, husbandry, and experimental procedures on animal welfare. Preference tests and consumer demand models demonstrate that animal choices are …

On Adaptive Estimation In Orthogonal Saturated Designs, Weizhen Wang, Daniel T. Voss

Mathematics and Statistics Faculty Publications

A simple method is provided to construct a general class of individual and simultaneous confidence intervals for the effects in orthogonal saturated designs. These intervals use the data adaptively, maintain the confidence levels sharply at 1 - α at the least favorable parameter configuration, work effectively under effect sparsity, and include the intervals by Wang and Voss (2001) as a special case.

Robust Regression With High Coverage, David J. Olive, Douglas M. Hawkins

Articles and Preprints

An important parameter for several high breakdown regression algorithm estimators is the number of cases given weight one, called the coverage of the estimator. Increasing the coverage is believed to result in a more stable estimator, but the price paid for this stability is greatly decreased resistance to outliers. A simple modification of the algorithm can greatly increase the coverage and hence its statistical performance while maintaining high outlier resistance.

What Is A Reasonable Attorney Fee? An Empirical Study Of Class Action Settlements, Theodore Eisenberg, Geoffrey P. Miller

Cornell Law Faculty Publications

Determining an appropriate fee is a difficult task facing trial court judges in class action litigation. But courts rarely rely on empirical research to assess a fee’s reasonableness, due, at least in part, to the relative paucity of available information. Existing empirical studies of attorney fees in class action cases are limited in scope, and generally do not control for important variables. To help fill this gap, we analyzed data from all state and federal class actions with reported fee decisions from 1993 to 2002 in which the fee and class recovery could be determined with reasonable confidence.

We find …

On The Stability Of The Positive Radial Steady States For A Semilinear Cauchy Problem, Yinbin Deng, Yi Li, Yi Liu

Mathematics and Statistics Faculty Publications

No abstract provided.

Geographic Variation In The Morphology Of Crotalus Horridus (Serpentes: Viperidae), John Robert Allsteadt

Biological Sciences Theses & Dissertations

The Timber Rattlesnake (Crotalus horridus) occurs in discontinuous populations throughout the eastern and central United States. The species exhibits high levels of polymorphism in morphological traits, especially in coloration and pattern. Previous studies recognized either distinct northern and southern subspecies or three regional morphs (northern, southern, and western), but conflicting data sets and limited geographic sampling of previous studies have left the relationships among those regional variants unclear. In this study, univariate and multivariate statistics, together with a geographic information system, were used to analyze geographic variation in 36 morphological characters recorded from 2,420 specimens of C. horridus …

Estimation Of Parameters In Replicated Time Series Regression Models, Genming Shi

Mathematics & Statistics Theses & Dissertations

The time series regression model was widely studied in the literature by several authors. However, statistical analysis of replicated time series regression models has received little attention. In this thesis, we study the application of quasi-least squares, a relatively new method, to estimate the parameters in replicated time series models with general ARMA( p, q) correlation structure. We also study several established methods for estimating the parameters in those models, including the maximum likelihood, method of moments, and the GEE method. Asymptotic comparisons of the methods are made bV fixing the number of repeated measurements in each series, and …

Analysis Of Multivariate Data Using Kotz Type Distribution, Kusaya Plungpongpun

Mathematics & Statistics Theses & Dissertations

Most of the inferential statistical methods for multivariate data are developed under the fundamental assumption that the data are from a multivariate normal distribution. Unfortunately, one can never be sure a set of data is really from a multivariate normal distribution. There are numerous methods for checking (testing) multivariate normality, but based on many published and our own simulation studies, provided in the first chapter of this dissertation, we observe that these tests are generally not very powerful, especially for smaller sample sizes. Hence it is always beneficial to have alternative multivariate distributions available along with the methodology for using …

Resampling-Based Multiple Testing: Asymptotic Control Of Type I Error And Applications To Gene Expression Data, Katherine S. Pollard, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

We define a general statistical framework for multiple hypothesis testing and show that the correct null distribution for the test statistics is obtained by projecting the true distribution of the test statistics onto the space of mean zero distributions. For common choices of test statistics (based on an asymptotically linear parameter estimator), this distribution is asymptotically multivariate normal with mean zero and the covariance of the vector influence curve for the parameter estimator. This test statistic null distribution can be estimated by applying the non-parametric or parametric bootstrap to correctly centered test statistics. We prove that this bootstrap estimated null …

Statistical Implications Of Pooling Rna Samples For Microarray Experiments, Xuejun Peng, Constance L. Wood, Eric M. Blalock, Kuey Chu Chen, Philip W. Landfield, Arnold J. Stromberg

Statistics Faculty Publications

BACKGROUND: Microarray technology has become a very important tool for studying gene expression profiles under various conditions. Biologists often pool RNA samples extracted from different subjects onto a single microarray chip to help defray the cost of microarray experiments as well as to correct for the technical difficulty in getting sufficient RNA from a single subject. However, the statistical, technical and financial implications of pooling have not been explicitly investigated.

RESULTS: Modeling the resulting gene expression from sample pooling as a mixture of individual responses, we derived expressions for the experimental error and provided both upper and lower bounds for …

The Consequences Of Race-Blindness: Revisiting Prediction Models With Current Law School Data, Linda F. Wightman

Journal of Legal Education

No abstract provided.

Maximization By Parts In Likelihood Inference, Peter Xuekun Song, Yanqin Fan, Jack Kalbfleisch

The University of Michigan Department of Biostatistics Working Paper Series

This paper presents and examines a new algorithm for solving a score equation for the maximum likelyhood estimate in certain problems of practical interest. The method circumvents the need to compute second order derivaties of the full likelihood function. It exploits the structure of certain models that yield a natural decomposition of a very complicated likelihood function. In this decomposition, the first part is a log likelihood from a simply analyzed model and the second part is used to update estimates from the first. Convergence properties of this fixed point algorithm are examined and asymptotics are derived for estimators obtained …

Design Considerations For Efficient And Effective Microarray Studies, M. Kathleen Kerr

UW Biostatistics Working Paper Series

This paper describes the theoretical and practical issues in experimental design for gene expression microarrays. Specifically, this paper (1) discusses the basic principles of design (randomization, replication, and blocking) as they pertain to microarrays, and (2) provides some general guidelines for statisticians designing microarray studies.

Cluster Stability Scores For Microarray Data In Cancer Studies, Mark Smolkin, Debashis Ghosh

The University of Michigan Department of Biostatistics Working Paper Series

A potential benefit of profiling of tissue samples using microarrays is the generation of molecular fingerprints that will define subtypes of disease. Hierarchical clustering has been the primary analytical tool used to define disease subtypes from microarray experiments in cancer settings. Assessing cluster reliability poses a major complication in analyzing output from these procedures. While much work has been done on assessing the global question of number of clusters in a dataset, relatively little research exists on assessing stability of individual clusters. A potential benefit of profiling of tissue samples using microarrays is the generation of molecular fingerprints that will …

Double Robust Estimation In Longitudinal Marginal Structural Models, Zhuo Yu, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

Consider estimation of causal parameters in a marginal structural model for the discrete intensity of the treatment specific counting process (e.g. hazard of a treatment specific survival time) based on longitudinal observational data on treatment, covariates and survival. We assume the sequential randomization assumption (SRA) on the treatment assignment mechanism and the so called experimental treatment assignment assumption which is needed to identify the causal parameters from the observed data distribution. Under SRA, the likelihood of the observed data structure factorizes in the auxiliary treatment mechanism and the partial likelihood consisting of the product over time of conditional distributions of …