Digital Commons Network^™

Properties Of The Gar(1) Model For Time Series Of Counts, Vasiliki Karioti, Chrys Caroni

Journal of Modern Applied Statistical Methods

Models for time series count data include several proposed by Zeger and Qaqish (1988), subsequently generalized into the GARMA family. The GAR(1) model is examined in detail. The maximum likelihood estimation of the parameters will be discussed and the properties of Pearson and randomized residuals will be examined.

Variance Estimation And Construction Of Confidence Intervals For Gee Estimator, Shenghai Zhang, Mary E. Thompson

Journal of Modern Applied Statistical Methods

The sandwich estimator, also known as the robust covariance matrix estimator, has achieved increasing use in the statistical literature as well as with the growing popularity of generalized estimating equations (GEE). A modified sandwich variance estimator is proposed, and its consistency and efficiency are studied. It is compared with other variance estimators, such as a model based estimator, the sandwich estimator and a corrected sandwich estimator. Confidence intervals for regression parameters based on these estimators are discussed. Simulation studies using clustered data to compare the performance of variance estimators are reported.

Jmasm22: A Convenient Way Of Generating Normal Random Variables Using Generalized Exponential Distribution, Debasis Kundu, Anubhav Manglick

Journal of Modern Applied Statistical Methods

A convenient method to generate normal random variable using a generalized exponential distribution is proposed. The new method is compared with the other existing methods and it is observed that the proposed method is quite competitive with most of the existing methods in terms of the K − S distances and the corresponding p-values.

A Combined Individuals And Moving Range Control Chart, Michael B. C. Khoo, S. H. Quah, C. K. Ch'ng

Journal of Modern Applied Statistical Methods

An individuals control chart is usually used to monitor shifts in the process mean when it is not possible to form subgroups. The moving range of two successive process measures is used as the basis for estimating the process variability. Similar to the case of the X − R and X − S charts, the individualsmoving range (I-MR) charts are used simultaneously in the monitoring of the process mean and variance respectively for individual observations, requiring maintaining two different charts. In this article, a new approach is suggested where the measurements of both the process mean and variance are plotted …

A Combined Standard Deviation Based Data Clustering Algorithm, Kuttiannan Thangavel, Durairaj Ashok Kumar

Journal of Modern Applied Statistical Methods

The clustering problem has been widely studied because it arises in many knowledge management oriented applications. It aims at identifying the distribution of patterns and intrinsic correlations in data sets by partitioning the data points into similarity clusters. Traditional clustering algorithms use distance functions to measure similarity centroid, which subside the influences of data points. Hence, in this article a novel non-distance based clustering algorithm is proposed which uses Combined Standard Deviation (CSD) as measure of similarity. The performance of CSD based K-means approach, called K-CSD clustering algorithm, is tested on synthetic data sets. It compared favorably to widely used …

The Use Of Hierarchical Ancova In Curriculum Studies, Show-Mann Liou, Chao-Ying Joanne Peng

Journal of Modern Applied Statistical Methods

Many educational studies are carried out in intact settings, such as classrooms or groups in which individual data were collected before and after a treatment. Researchers advocate either the use of individual scores as the unit of analysis or class means. Both approaches suffer from conceptual and methodological limitations. In this article, the use of hierarchical ANCOVA for analyzing quasiexperimental data including baseline measures is designed and promoted. It is illustrated with a realworld data set collected from a curriculum study. Results showed that the hierarchical ANCOVA is a conceptually and methodologically sound approach, and is better than ANCOVA based …

Programs To Compute Distribution Functions And Critical Values For Extreme Value Ratios For Outlier Detection, George Mcbane

Peer Reviewed Articles

A set of FORTRAN subprograms is presented to compute density and cumulative distribution functions and critical values for the range ratio statistics of Dixon (1951, The Annals of Mathematical Statistics) These statistics are useful for detection of outliers in small samples

Profile Likelihood Estimation Of Partially Linear Panel Data Models With Fixed Effects, Liangjun Su, Aman Ullah

Research Collection School Of Economics

We consider consistent estimation of partially linear panel data models with fixed effects. We propose profile-likelihood-based estimators for both the parametric and nonparametric components in the models and establish convergence rates and asymptotic normality for both estimators.

Semiparametric Latent Variable Regression Models For Spatio-Temporal Modeling Of Mobile Source Particles In The Greater Boston Area, Alexandros Gryparis, Brent A. Coull, Joel Schwartz, Helen H. Suh

Harvard University Biostatistics Working Paper Series

Traffic particle concentrations show considerable spatial variability within a metropolitan area. We consider latent variable semiparametric regression models for modeling the spatial and temporal variability of black carbon and elemental carbon concentrations in the greater Boston area. Measurements of these pollutants, which are markers of traffic particles, were obtained from several individual exposure studies conducted at specific household locations as well as 15 ambient monitoring sites in the city. The models allow for both flexible, nonlinear effects of covariates and for unexplained spatial and temporal variability in exposure. In addition, the different individual exposure studies recorded different surrogates of traffic …

Estimating The Integrated Likelihood Via Posterior Simulation Using The Harmonic Mean Identity, Adrian E. Raftery, Michael A. Newton, Jaya M. Satagopan, Pavel N. Krivitsky

Memorial Sloan-Kettering Cancer Center, Dept. of Epidemiology & Biostatistics Working Paper Series

The integrated likelihood (also called the marginal likelihood or the normalizing constant) is a central quantity in Bayesian model selection and model averaging. It is defined as the integral over the parameter space of the likelihood times the prior density. The Bayes factor for model comparison and Bayesian testing is a ratio of integrated likelihoods, and the model weights in Bayesian model averaging are proportional to the integrated likelihoods. We consider the estimation of the integrated likelihood from posterior simulation output, aiming at a generic method that uses only the likelihoods from the posterior simulation iterations. The key is the …

Super Learning: An Application To Prediction Of Hiv-1 Drug Susceptibility, Sandra E. Sinisi, Maya L. Petersen, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

Many statistical methods exist that can be used to learn a predictor based on observed data. Examples include decision trees, neural networks, support vector regression, least angle regression, Logic Regression, and the Deletion/Substitution/Addition algorithm. The optimal algorithm for prediction will vary depending on the underlying data-generating distribution. In this article, we introduce a "super learner," a prediction algorithm that applies any set of candidate learners and uses cross-validation to select among them. Theory shows that asymptotically the super learner performs essentially as well or better than any of the candidate learners. We briefly present the theory behind the super learner, …

Recurrent Event Models In The Presence Of A Terminal Event: Comparison, Inference And Data Analysis, Xianghua Luo, Mei-Cheng Wang

Johns Hopkins University, Dept. of Biostatistics Working Papers

This article focuses on statistical implications of proportional rate models for recurrent event data in the presence of a terminal event. In such circumstances, various definitions of the recurrent rate function have been adopted in the proportional rate models. Although these rate functions have quite different interpretations, recognition of the differences has been lacking theoretically and practically. We compare three types of rate functions from both conceptual and quantitative perspectives; conclude that the inappropriate choice of a rate function may lead to misleading scientific conclusions. Simulations are conducted for comparisons of the focused models. Analysis of data from an AIDS …

Empirical Bayes Approach To Controlling Familywise Error: An Application To Hiv Resistance Data, Rhoderick N. Machekano, Alan E. Hubbard

U.C. Berkeley Division of Biostatistics Working Paper Series

Statistical challenges arise in identifying meaningful patterns and structures from high dimensional genomic data sets. Relating HIV genotype (sequence of amino acids) to phenotypic resistance presents a typical problem. When the HIV virus is under antiretroviral drug pressure, unfavorable mutations of the target genes often lead to greatly increased resistance of the virus to drugs, including drugs the virus has not been exposed to. Identification of mutation combinations and their correlation to drug resistance is critical in guiding efficient prescription of HIV drugs. The identification of a subset of codons associated with drug resistance from a set of several hundreds …

Survival Analysis With Change Point Hazard Functions, Melody S. Goodman, Yi Li, Ram C. Tiwari

Harvard University Biostatistics Working Paper Series

No abstract provided.

Modeling And Simulation Of Value -At -Risk In The Financial Market Area, Xiangyin Zheng

Doctoral Dissertations

Value-at-Risk (VaR) is a statistical approach to measure market risk. It is widely used by banks, securities firms, commodity and energy merchants, and other trading organizations. The main focus of this research is measuring and analyzing market risk by modeling and simulation of Value-at-Risk for portfolios in the financial market area. The objectives are (1) predicting possible future loss for a financial portfolio from VaR measurement, and (2) identifying how the distributions of the risk factors affect the distribution of the portfolio. Results from (1) and (2) provide valuable information for portfolio optimization and risk management.

The model systems chosen …

Measuring Inequality: Statistical Inference Theory With Applications, Mihaela Paun

Doctoral Dissertations

In this dissertation we develop statistical inference for the Atkinson index, one of the measures of inequality used in studying economic inequality.

Specifically, we construct empirical estimators for the Atkinson index, both in the parametric and nonparametric case, and derive formulas for the asymptotic variances for the estimators. These statistics are used for testing hypothesis and constructing confidence intervals for the Atkinson index. We test the validity and the robustness of the asymptotic theory, by simulations (using R, a language and environment for statistical computing and graphics), in the case of one and two populations. In addition to proving asymptotic …

New Estimators Of A Circular Median, Sauwanit Ratanaruamkarn

Dissertations

The specific properties of probability distributions on a circle require different definitions of several statistical concepts. For example, a median that can always be found for linear data not always exists on the circle. Several estimators of a circular median were proposed, lately by Otieno (2002), Otieno and Anderson - Cook (2003). Their work is, however, focused on the "preferred direction" which coincides with the median, mean and mode in the case of symmetric, unimodal distributions. This dissertation is focused on the estimators of a population median in a wider range of population distributions on a circle including distributions with …

Causal Effect Models For Intention To Treat And Realistic Individualized Treatment Rules, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

An important class of models in causal inference are the so-called marginal structural models which model the comparison between counterfactual outcome distributions corresponding with a static treatment intervention, conditional on user supplied baseline covariates, based on observing a longitudinal data structure on a sample of n independent and identically distributed experimental units. Identification of a static treatment regimen specific outcome distribution based on observational data requires beyond the so-called sequential randomization assumption that each experimental unit has positive probability of following the static treatment regimen. The latter assumption is called the experimental treatment assignment assumption (ETA) (which is parameter specific). …

Reliability, Effect Size, And Responsiveness And Intraclass Correlation Of Health Status Measures Used In Randomized And Cluster-Randomized Trials, Paula Diehr, Lu Chen, Donald L. Patrick, Ziding Feng, Yutaka Yasui

UW Biostatistics Working Paper Series

Background: New health status instruments are described by psychometric properties, such as Reliability, Effect Size, and Responsiveness. For cluster-randomized trials, another important statistic is the Intraclass Correlation for the instrument within clusters. Studies using better instruments can be performed with smaller sample sizes, but better instruments may be more expensive in terms of dollars, lost opportunities, or poorer data quality due to the response burden of longer instruments. Investigators often need to estimate the psychometric properties of a new instrument, or of an established instrument in a new setting. Optimal sample sizes for estimating these properties have not been studied …

Detecting Pulsatile Hormone Secretion Events: A Bayesian Approach, Tim Johnson

The University of Michigan Department of Biostatistics Working Paper Series

Many challenges arise in the analysis of pulsatile, or episodic, hormone concentration time series data. Among these challenges is the determination of the number and location of pulsatile events and the discrimination of events from noise. Analyses of these data are typically performed in two stages. In the first stage, the number and approximate location of the pulses are determined. In the second stage, a model (typically a deconvolution model) is fit to the data conditional on the number of pulses. Any error made in the first stage is carried over to the second stage. Furthermore, current methods, except two, …

Semiparametric Analysis For Correlated Recurrent And Terminal Events, Yining Ye, Jack Kalbfleisch, Doug E. Schaubel

The University of Michigan Department of Biostatistics Working Paper Series

In clinical and observational studies, recurrent event data (e.g. hospitalization) with a terminal event (e.g. death) are often encountered. In many instances, the terminal event is strongly correlated with the recurrent event process. In this article, we propose a semiparametric method to jointly model the recurrent and terminal event processes. The dependence is modeled by a shared gamma frailty that is included in both the recurrent event rate and terminal event hazard function. Marginal models are used to estimate the regression effects on the terminal and recurrent event processes and a Poisson model is used to estimate the dispersion of …

Adjusting For Covariate Effects On Classification Accuracy Using The Covariate-Adjusted Roc Curve, Holly Janes, Margaret S. Pepe

UW Biostatistics Working Paper Series

Recent scientific and technological innovations have produced an abundance of potential markers which are being investigated for their use in disease screen- ing and diagnosis. In evaluating these markers, it is often necessary to account for covariates which are associated with the marker of interest. These covariates may include subject characteristics, expertise of the test operator, test proce- dures, or aspects of specimen handling. In this paper, we propose the AROC, a covariate-adjusted measure of the classification accuracy. The AROC is the common covariate-specific ROC curve, when the covariate does not affect dis- crimination, and a weighted average of covariate-specific …

Censored Data Regression In High-Dimension And Low-Sample Size Settings For Genomic Applications, Hongzhe Li

UPenn Biostatistics Working Papers

New high-throughput technologies are generating various types of high-dimensional genomic and proteomic data and meta-data (e.g., networks and pathways) in order to obtain a systems-level understanding of various complex diseases such as human cancers and cardiovascular diseases. As the amount and complexity of the data increase and as the questions being addressed become more sophisticated, we face the great challenge of how to model such data in order to draw valid statistical and biological conclusions. One important problem in genomic research is to relate these high-throughput genomic data to various clinical outcomes, including possibly censored survival outcomes such as age …

A General Framework For Statistical Performance Comparison Of Evolutionary Computation Algorithms, David Shilane, Jarno Martikainen, Sandrine Dudoit, Seppo Ovaska

U.C. Berkeley Division of Biostatistics Working Paper Series

This paper proposes a statistical methodology for comparing the performance of evolutionary computation algorithms. A two-fold sampling scheme for collecting performance data is introduced, and these data are analyzed using bootstrap-based multiple hypothesis testing procedures. The proposed method is sufficiently flexible to allow the researcher to choose how performance is measured, does not rely upon distributional assumptions, and can be extended to analyze many other randomized numeric optimization routines. As a result, this approach offers a convenient, flexible, and reliable technique for comparing algorithms in a wide variety of applications.

The Two-Sample Problem For Failure Rates Depending On A Continuous Mark: An Application To Vaccine Efficacy, Peter B. Gilbert, Ian W. Mckeague, Yanqing Sun

UW Biostatistics Working Paper Series

The efficacy of an HIV vaccine to prevent infection is likely to depend on the genetic variation of the exposing virus. This paper addresses the problem of using data on the HIV sequences that infect vaccine efficacy trial participants to 1) test for vaccine efficacy more powerfully than procedures that ignore the sequence data; and 2) evaluate the dependence of vaccine efficacy on the divergence of infecting HIV strains from the HIV strain that is contained in the vaccine. Because hundreds of amino acid sites in each HIV genome are sequenced, it is natural to treat the divergence (defined in …

Genome Scanning Methods For Comparing Sequences Between Groups, With Application To Hiv Vaccine Trials, Peter B. Gilbert, Chunyuan Wu, David V. Jobes

UW Biostatistics Working Paper Series

Consider a placebo-controlled preventive HIV vaccine efficacy trial. An HIV amino acid sequence is measured from each volunteer who acquires HIV, and these sequences are aligned together with the reference HIV sequence represented in the vaccine. We develop genome scanning methods to identify HIV positions at which the amino acids in sequences from infected vaccine recipients tend to be more divergent from the corresponding reference amino acid than the amino acids in sequences from infected placebo recipients. We consider five two-sample test statistics, based on Euclidean, Mahalanobis, and Kullback-Leibler divergence measures. Weights are incorporated to reflect biological information contained in …

On The Equivalence Of Case-Crossover And Time Series Methods In Environmental Epidemiology, Yun Lu, Scott L. Zeger

Johns Hopkins University, Dept. of Biostatistics Working Papers

Time series and case-crossover methods are often viewed as competing alternatives in environmental epidemiologic studies. Several recent studies have compared the time series and case-crossover methods. In this paper, we show that case-crossover using conditional logistic regression is a special case of time series analysis when there is a common exposure such as in air pollution studies. This equivalence provides computational convenience for case-crossover analyses and a better understanding of time series models. Time series log-linear regression accounts for over-dispersion of the Poisson variance, while case-crossover analyses typically do not. This equivalence also permits model checking for case-crossover data using …

Evaluating Prediction Rules For T-Year Survivors With Censored Regression Models, Hajime Uno, Tianxi Cai, Lu Tian, L.J. Wei

Harvard University Biostatistics Working Paper Series

Suppose that we are interested in establishing simple, but reliable rules for predicting future t-year survivors via censored regression models. In this article, we present inference procedures for evaluating such binary classification rules based on various prediction precision measures quantified by the overall misclassification rate, sensitivity and specificity, and positive and negative predictive values. Specifically, under various working models we derive consistent estimators for the above measures via substitution and cross validation estimation procedures. Furthermore, we provide large sample approximations to the distributions of these nonsmooth estimators without assuming that the working model is correctly specified. Confidence intervals, for example, …

The Longitudinal Effect Of Self-Monitoring And Locus Of Control On Social Network Position In Friendship Networks, Gary J. Moore

Theses and Dissertations

The purpose of this research was to identify how enduring personality characteristics predict a person's location in a network, locations which in turn affect outcomes such as performance. Specifically, this thesis examines how self-monitoring and locus of control influence an individual's location in a friendship social network over time. Hierarchical Linear Modeling (HLM) was used to analyze 28 groups of students and instructors at a military training course over six and one half weeks. Self-monitoring predicted betweenness centrality in five of six time periods while locus of control predicted betweenness centrality in three of six time periods. The moderation of …

Modeling The Performance Of A Baseball Player's Offensive Production, Michael Ross Smith

Theses and Dissertations

This project addresses the problem of comparing the offensive abilities of players from different eras in Major League Baseball (MLB). We will study players from the perspective of an overall offensive summary statistic that is highly linked with scoring runs, or the Berry Value. We will build an additive model to estimate the innate ability of the player, the effect of the relative level of competition of each season, and the effect of age on performance using piecewise age curves. Using Hierarchical Bayes methodology with Gibbs sampling, we model each of these effects for each individual. The results of the …