Left truncation in linked data: A practical guide to understanding left truncation and applying it using SAS and R

Time-to-event data such as time to death are broadly used in medical research and drug development to understand the efficacy of a therapeutic. For time-to-event data, right censoring (data only observed up to a certain point of time) is common and easy to recognize. Methods that use right censored data, such as the Kaplan–Meier estimator and the Cox proportional hazard model, are well established. Time-to-event data can also be left truncated, which arises when patients are excluded from the sample because their events occur before a specific milestone, potentially resulting in an immortal time bias. For example, in a study evaluating the association between biomarker status and overall survival, patients who did not live long enough to receive a genomic test were not observed in the study. Left truncation causes selection bias and often leads to an overestimate of survival time. In this tutorial, we used a nationwide electronic health record-derived de-identified database to demonstrate how to analyze left truncated and right censored data without bias using example code from SAS and R.     

Stochastic hydrology: An introduction to wet statistics for dry statisticians

The paper is intended to give non-initiates some idea of the nature of stochastic hydrology. After a brief historical review it concentrates on three recent examples of stochastic modelling procedures that have aroused interest among hydrologists, namely the Hurst Effect, Short-Term Runoff Models, and Stochastic Reservoir Theory.     

Computation of exact confidence intervals from discrete data using studentized test statistics

A new area of research interest is the computation of exact confidence limits or intervals for a scalar parameter of interest from discrete data by inverting a hypothesis test based on a studentized test statistic. See, for example, Chan and Zhang (1999), Agresti and Min (2001) and Agresti (2003) who deal with a difference of binomial probabilities and Agresti and Min (2002) who deal with an odds ratio. However, neither (1) a detailed analysis of the computational issues involved nor (2) a reliable method of computation that deals effectively with these issues is currently available. In this paper we solve these two problems for a very broad class of discrete data models. We suppose that the distribution of the data is determined by (,) where is a nuisance parameter vector. We also consider six different studentized test statistics. Our contributions to (1) are as follows. We show that the P-value resulting from the hypothesis test, considered as a function of the null-hypothesized value of , has both jump and drop discontinuities. Numerical examples are used to demonstrate that these discontinuities lead to the failure of simple-minded approaches to the computation of the confidence limit or interval. We also provide a new method for efficiently computing the set of all possible locations of these discontinuities. Our contribution to (2) is to provide a new and reliable method of computing the confidence limit or interval, based on the knowledge of this set.     

Analysis of means: a generalized approach using R

Papers on the analysis of means (ANOM) have been circulating in the quality control literature for decades, routinely describing it as a statistical stand-alone concept. Therefore, we clarify that ANOM should rather be regarded as a special case of a much more universal approach known as multiple contrast tests (MCTs). Perceiving ANOM as a grand-mean-type MCT paves the way for implementing it in the open-source software R. We give a brief tutorial on how to exploit R's versatility and introduce the R package ANOM for drawing the familiar decision charts. Beyond that, we illustrate two practical aspects of data analysis with ANOM: firstly, we compare merits and drawbacks of ANOM-type MCTs and ANOVA F-test and assess their respective statistical powers, and secondly, we show that the benefit of using critical values from multivariate t-distributions for ANOM instead of simple Bonferroni quantiles is oftentimes negligible.     

The estimation of R 2 and adjusted R 2 in incomplete data sets using multiple imputation

The coefficient of determination, known also as the R ², is a common measure in regression analysis. Many scientists use the R ² and the adjusted R ² on a regular basis. In most cases, the researchers treat the coefficient of determination as an index of ‘usefulness’ or ‘goodness of fit,’ and in some cases, they even treat it as a model selection tool. In cases in which the data is incomplete, most researchers and common statistical software will use complete case analysis in order to estimate the R ², a procedure that might lead to biased results. In this paper, I introduce the use of multiple imputation for the estimation of R ² and adjusted R ² in incomplete data sets. I illustrate my methodology using a biomedical example.     

Some results on estimating a change-point using non-parametric type statistics

A Mann-Whitney type statistic is used to estimate a change-point when a change, at an unknown point in a sequence of random variables, has taken place. This estimate is compared, using Monte Carlo techniques, with the normal theory maximum likelihood estimate, when a location change has occurred, for different underlying distributions ranging from the normal to the long tailed “normal over uniform” distribution. The distribution of the Mann-Whitney type estimate remains fairly constant over the various distributions. Two generalisations of the statistic are considered and investigated.     

On the role of data,statistics and decisions in a pandemic

AStA Advances in Statistical Analysis - A pandemic poses particular challenges to decision-making because of the need to continuously adapt decisions to rapidly changing evidence and available...     

Analyzing panel count data with a dependent observation process and a terminal event

Panel count data occur in many fields and a number of approaches have been developed. However, most of these approaches are for situations where there is no terminal event and the observation process is independent of the underlying recurrent event process unconditionally or conditional on the covariates. In this paper, we discuss a more general situation where the observation process is informative and there exists a terminal event which precludes further occurrence of the recurrent events of interest. For the analysis, a semiparametric transformation model is presented for the mean function of the underlying recurrent event process among survivors. To estimate the regression parameters, an estimating equation approach is proposed in which an inverse survival probability weighting technique is used. The asymptotic distribution of the proposed estimates is provided. Simulation studies are conducted and suggest that the proposed approach works well for practical situations. An illustrative example is provided. The Canadian Journal of Statistics 41: 174–191; 2013 © 2012 Statistical Society of Canada     

On using percentiles to fit data by a johnson distribution

The well-known Johnson system of distributions was developed by N. L. Johnson (1949). Slifker and Shapiro (1980) presented a criterion for choosing a member from the three distributional classes (S_B,S_L, and S_v) in the Johnson system to fit a set of data. The criterion is based on the value of a quantile ratio which depends on a specified positive z value and the parameters of the distribution. In this paper, we present some properties of the quantile ratio for various distributions and for some selected z values. Some comments are made on using the criterion for selecting a Johnson distribution to fit empirical data.     

Fitting the generalized Pareto distribution to data based on transformations of order statistics

Generalized Pareto distribution (GPD) has been widely used to model exceedances over thresholds. In this article we propose a new method called weighted nonlinear least squares (WNLS) to estimate the parameters of the GPD. The WNLS estimators always exist and are simple to compute. Some asymptotic results of the proposed method are provided. The simulation results indicate that the proposed method performs well compared to existing methods in terms of mean squared error and bias. Its advantages are further illustrated through the analysis of two real data sets.     

Analyzing multivariate longitudinal binary data: A generalized estimating equations approach

The authors consider regression analysis for binary data collected repeatedly over time on members of numerous small clusters of individuals sharing a common random effect that induces dependence among them. They propose a mixed model that can accommodate both these structural and longitudinal dependencies. They estimate the parameters of the model consistently and efficiently using generalized estimating equations. They show through simulations that their approach yields significant gains in mean squared error when estimating the random effects variance and the longitudinal correlations, while providing estimates of the fixed effects that are just as precise as under a generalized penalized quasi‐likelihood approach. Their method is illustrated using smoking prevention data.     

Scan statistics for detecting a local change in variance for two-dimensional normal data

In this article scan statistics for detecting a local change in variance for two-dimensional normal data are discussed. When the precise size of the rectangular window, where a local change in variance has occurred, is unknown, multiple and variable window scan statistics are proposed. A simulation study is presented to evaluate the performance of the scan statistics investigated in this article via comparison of power. A method for estimating the rectangular region, where a change in variance has occurred, and the size of the change in variance is also discussed.     

Analyzing survival data with highly negatively skewed distribution: The Gompertz-sinh family

In this article, we explore a new two-parameter family of distribution, which is derived by suitably replacing the exponential term in the Gompertz distribution with a hyperbolic sine term. The resulting new family of distribution is referred to as the Gompertz-sinh distribution, and it possesses a thicker and longer lower tail than the Gompertz family, which is often used to model highly negatively skewed data. Moreover, we introduce a useful generalization of this model by adding a second shape parameter to accommodate a variety of density shapes as well as nondecreasing hazard shapes. The flexibility and better fitness of the new family, as well as its generalization, is demonstrated by providing well-known examples that involve complete, group, and censored data.     

Discussion on On the role of data,statistics and decisions in a pandemic

AStA Advances in Statistical Analysis - The authors make an important contribution presenting a comprehensive and thoughtful overview about the many different aspects of data, statistics and data...