Local Box–Cox transformation on time-varying parametric models for smoothing estimation of conditional CDF with longitudinal data

Nonparametric estimation and inferences of conditional distribution functions with longitudinal data have important applications in biomedical studies, such as epidemiological studies and longitudinal clinical trials. Estimation approaches without any structural assumptions may lead to inadequate and numerically unstable estimators in practice. We propose in this paper a nonparametric approach based on time-varying parametric models for estimating the conditional distribution functions with a longitudinal sample. Our model assumes that the conditional distribution of the outcome variable at each given time point can be approximated by a parametric model after local Box–Cox transformation. Our estimation is based on a two-step smoothing method, in which we first obtain the raw estimators of the conditional distribution functions at a set of disjoint time points, and then compute the final estimators at any time by smoothing the raw estimators. Applications of our two-step estimation method have been demonstrated through a large epidemiological study of childhood growth and blood pressure. Finite sample properties of our procedures are investigated through a simulation study. Application and simulation results show that smoothing estimation from time-variant parametric models outperforms the existing kernel smoothing estimator by producing narrower pointwise bootstrap confidence band and smaller root mean squared error.     

Effect of individual observations on the Box–Cox transformation

In this paper, we consider the influence of individual observations on inferences about the Box–Cox power transformation parameter from a Bayesian point of view. We compare Bayesian diagnostic measures with the ‘forward’ method of analysis due to Riani and Atkinson. In particular, we look at the effect of omitting observations on the inference by comparing particular choices of transformation using the conditional predictive ordinate and the k _d measure of Pettit and Young. We illustrate the methods using a designed experiment. We show that a group of masked outliers can be detected using these single deletion diagnostics. Also, we show that Bayesian diagnostic measures are simpler to use to investigate the effect of observations on transformations than the forward search method.     

An algorithm for estimating Box–Cox transformation parameter in ANOVA

In this study, we construct a feasible region, in which we maximize the likelihood function, by using Shapiro–Wilk and Bartlett's test statistics to obtain Box–Cox power transformation parameter for solving the issues of non-normality and/or heterogeneity of variances in analysis of variance (ANOVA). Simulation studies illustrate that the proposed approach is more successful in attaining normality and variance stabilization, and is at least as good as the usual maximum likelihood estimation (MLE) in estimating the transformation parameter for different conditions. Our proposed method is illustrated on two real-life datasets. Moreover, the proposed algorithm is released under R package AID under the name of “boxcoxfr” for implementation.     

Variable selection in joint mean and variance models of Box–Cox transformation

In many applications, a single Box–Cox transformation cannot necessarily produce the normality, constancy of variance and linearity of systematic effects. In this paper, by establishing a heterogeneous linear regression model for the Box–Cox transformed response, we propose a hybrid strategy, in which variable selection is employed to reduce the dimension of the explanatory variables in joint mean and variance models, and Box–Cox transformation is made to remedy the response. We propose a unified procedure which can simultaneously select significant variables in the joint mean and variance models of Box–Cox transformation which provide a useful extension of the ordinary normal linear regression models. With appropriate choice of the tuning parameters, we establish the consistency of this procedure and the oracle property of the obtained estimators. Moreover, we also consider the maximum profile likelihood estimator of the Box–Cox transformation parameter. Simulation studies and a real example are used to illustrate the application of the proposed methods.     

Non parametric confidence intervals for the extinction probability of a Galton–Watson branching process

It is demonstrated that the confidence intervals (CIs) for the probability of eventual extinction and other parameters of a Galton–Watson branching process based upon the maximum likelihood estimators can often have substantially lower coverage when compared to the desired nominal confidence coefficient, especially in small, more realistic sample sizes. The same conclusion holds for the traditional bootstrap CIs. We propose several adjustments to these CIs, which greatly improves coverage in most cases. We also make a correction in an asymptotic variance formula given in Stigler (1971 Stigler, S.M. (1971). The estimation of the probability of extinction and other parameters associated with branching processes. Biometrika 58(3):499–508.[Crossref], [Web of Science ®] , [Google Scholar]). The focus here is on implementation of the CIs which have good coverage, in a wide variety of cases. We also consider expected CI lengths. Some recommendations are made.     

Box–Cox transformation for spatial linear models: a study on lattice data

In this article, we extended the classic Box–Cox transformation to spatial linear models. For a comparative study, the proposed models were applied to a real data set of Chinese population growth and economic development with three different structures: no spatial correction, conditional autoregressive and simultaneous autoregressive. Maximal likelihood method was used to estimate the Box–Cox parameter λ and other parameters in the models. The residuals of the models were analyzed through Moran’s I and Geary’s c.     

Prediction Techniques for Box–Cox Regression Models

This article reviews several techniques useful for forming point and interval predictions in regression models with Box-Cox transformed variables. The techniques reviewed—plug-in, mean squared error analysis, predictive likelihood, and stochastic simulation—take account of nonnormality and parameter uncertainty in varying degrees. A Monte Carlo study examining their small-sample accuracy indicates that uncertainty about the Box–Cox transformation parameter may be relatively unimportant. For certain parameters, deterministic point predictions are biased, and plug-in prediction intervals are also biased. Stochastic simulation, as usually carried out, leads to badly biased predictions. A modification of the usual approach renders stochastic simulation predictions largely unbiased.     

Effects of the generalized Box–Cox transformation on Type I error rate and power of Hotelling's T 2

Most multivariate statistical techniques rely on the assumption of multivariate normality. The effects of nonnormality on multivariate tests are assumed to be negligible when variance–covariance matrices and sample sizes are equal. Therefore, in practice, investigators usually do not attempt to assess multivariate normality. In this simulation study, the effects of skewed and leptokurtic multivariate data on the Type I error and power of Hotelling's T ² were examined by manipulating distribution, sample size, and variance–covariance matrix. The empirical Type I error rate and power of Hotelling's T ² were calculated before and after the application of generalized Box–Cox transformation. The findings demonstrated that even when variance–covariance matrices and sample sizes are equal, small to moderate changes in power still can be observed.     

Improving TAIEX forecasting using fuzzy time series with Box–Cox power transformation

Box–Cox together with our newly proposed transformation were implemented in three different real world empirical problems to alleviate noisy and the volatility effect of them. Consequently, a new domain was constructed. Subsequently, universe of discourse for transformed data was established and an approach for calculating effective length of the intervals was then proposed. Considering the steps above, the initial forecasts were performed using frequently used fuzzy time series (FTS) methods on transformed data. Final forecasts were retrieved from initial forecasted values by proper inverse operation. Comparisons of the results demonstrate that the proposed method produced more accurate forecasts compared with existing FTS on original data.     

A Realized Stochastic Volatility Model With Box–Cox Transformation

This article presents a new class of realized stochastic volatility model based on realized volatilities and returns jointly. We generalize the traditionally used logarithm transformation of realized volatility to the Box–Cox transformation, a more flexible parametric family of transformations. A two-step maximum likelihood estimation procedure is introduced to estimate this model on the basis of Koopman and Scharth (2013 Koopman, S.J., Scharth, M. (2013), The Analysis of Stochastic Volatility in the Presence of Daily Realised Measures, Journal of Financial Econometrics, 11, 76–115.[Crossref], [Web of Science ®] , [Google Scholar]). Simulation results show that the two-step estimator performs well, and the misspecified log transformation may lead to inaccurate parameter estimation and certain excessive skewness and kurtosis. Finally, an empirical investigation on realized volatility measures and daily returns is carried out for several stock indices.     

A modification of the Box–Meyer method for finding the active factors in screening experiments

Screening experiments are conducted to identify a few active factors among a large number of factors. For the objective of identifying active factors, Box and Meyer provided an innovative approach, the Box–Meyer method (BMM). With the use of means models, we propose a modification of the BMM in this paper. Compared with the original BMM, the modified BMM (MBMM) can circumvent the problem that the original BMM runs into, namely that it may fail to identify some active factors due to the ignorance of higher order interactions. Furthermore, the number of explanatory variables in the MBMM is smaller. Therefore, the computational complexity is reduced. Finally, three examples with different types of designs are used to demonstrate the wide applicability of the MBMM.     

A two-stage estimation in the Clayton–Oakes model with marginal linear transformation models for multivariate failure time data

This paper considers the analysis of multivariate survival data where the marginal distributions are specified by semiparametric transformation models, a general class including the Cox model and the proportional odds model as special cases. First, consideration is given to the situation where the joint distribution of all failure times within the same cluster is specified by the Clayton–Oakes model (Clayton, Biometrika 65:141–151, l978; Oakes, J R Stat Soc B 44:412–422, 1982). A two-stage estimation procedure is adopted by first estimating the marginal parameters under the independence working assumption, and then the association parameter is estimated from the maximization of the full likelihood function with the estimators of the marginal parameters plugged in. The asymptotic properties of all estimators in the semiparametric model are derived. For the second situation, the third and higher order dependency structures are left unspecified, and interest focuses on the pairwise correlation between any two failure times. Thus, the pairwise association estimate can be obtained in the second stage by maximizing the pairwise likelihood function. Large sample properties for the pairwise association are also derived. Simulation studies show that the proposed approach is appropriate for practical use. To illustrate, a subset of the data from the Diabetic Retinopathy Study is used.     

An improved method of estimation for the parameters of the Birnbaum–Saunders distribution

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution, frequently used for analysing lifetime data. In this paper, we propose a simple method of estimation for the parameters of the two-parameter BS distribution by making use of some key properties of the distribution. Compared with the maximum likelihood estimators and the modified moment estimators, the proposed method has smaller bias, but having the same mean square errors as these two estimators. We also discuss some methods of construction of confidence intervals. The performance of the estimators is then assessed by means of Monte Carlo simulations. Finally, an example is used to illustrate the method of estimation developed here.     

Box–Cox t random intercept model for estimating usual nutrient intake distributions

The issue of estimating usual nutrient intake distributions and prevalence of inadequate nutrient intakes is of interest in nutrition studies. Box–Cox transformations coupled with the normal distribution are usually employed for modeling nutrient intake data. When the data present highly asymmetric distribution or include outliers, this approach may lead to implausible estimates. Additionally, it does not allow interpretation of the parameters in terms of characteristics of the original data and requires back transformation of the transformed data to the original scale. This paper proposes an alternative approach for estimating usual nutrient intake distributions and prevalence of inadequate nutrient intakes through a Box–Cox t model with random intercept. The proposed model is flexible enough for modeling highly asymmetric data even when outliers are present. Unlike the usual approach, the proposed model does not require a transformation of the data. A simulation study suggests that the Box–Cox t model with random intercept estimates the usual intake distribution satisfactorily, and that it should be preferable to the usual approach particularly in cases of highly asymmetric heavy-tailed data. In applications to data sets on intake of 19 micronutrients, the Box–Cox t models provided better fit than its competitors in most of the cases.     

Adjusting Cornish–Fisher expansions and confidence intervals for the effect of roundoff

Suppose we want to estimate some smooth function of two types of parameters. The first can be estimated by sample means, while the second is known exactly up to the number of decimal places recorded, that is they are subject to roundoff. We obtain the Cornish–Fisher expansions and associated nonparametric confidence intervals for such functions. These results are illustrated by a simulation study.     

A simple generalization of the Box–Muller method for obtaining a pair of correlated standard normal variables

In the simple and widely used method of Box–Muller [G. Box and M. Muller, A note on the generation of random normal deviates, Ann. Math. Statist. 29 (1958), pp. 610–611], from a pair of uniform and independent random variables in (0,1), a pair of standard and independent normal variables is obtained. In this article, we present a very simple and elegant generalization of this method to obtain a pair of correlated standard normal variables with a given coefficient of correlation. This generalized method, which is computationally very easy, is interpreted in geometric terms, considering a translation of the uniform interval (0,1) and a rotation of a defined angle, both related to the coefficient of correlation. Some numerical results are simulated and statistically analysed, proving that the generalization is extremely simple and powerful.     

Approximate tolerance intervals for the discrete Poisson–Lindley distribution

The Poisson–Lindley distribution is a compound discrete distribution that can be used as an alternative to other discrete distributions, like the negative binomial. This paper develops approximate one-sided and equal-tailed two-sided tolerance intervals for the Poisson–Lindley distribution. Practical applications of the Poisson–Lindley distribution frequently involve large samples, thus we utilize large-sample Wald confidence intervals in the construction of our tolerance intervals. A coverage study is presented to demonstrate the efficacy of the proposed tolerance intervals. The tolerance intervals are also demonstrated using two real data sets. The R code developed for our discussion is briefly highlighted and included in the tolerance package.     

Comparison of estimation methods for the Marshall–Olkin extended Lindley distribution

The aim of this paper is to compare the parameters' estimations of the Marshall–Olkin extended Lindley distribution obtained by six estimation methods: maximum likelihood, ordinary least-squares, weighted least-squares, maximum product of spacings, Cramér–von Mises and Anderson–Darling. The bias, root mean-squared error, average absolute difference between the true and estimate distributions' functions and the maximum absolute difference between the true and estimate distributions' functions are used as comparison criteria. Although the maximum product of spacings method is not widely used, the simulation study concludes that it is highly competitive with the maximum likelihood method.     

Non parametric estimation of the structural expectation of a stochastic increasing function

This article introduces a non parametric warping model for functional data. When the outcome of an experiment is a sample of curves, data can be seen as realizations of a stochastic process, which takes into account the variations between the different observed curves. The aim of this work is to define a mean pattern which represents the main behaviour of the set of all the realizations. So, we define the structural expectation of the underlying stochastic function. Then, we provide empirical estimators of this structural expectation and of each individual warping function. Consistency and asymptotic normality for such estimators are proved.     

A generalised likelihood framework for partially observed capture–recapture–recovery models

We provide a closed form likelihood expression for multi-state capture–recapture–recovery data when the state of an individual may be only partially observed. The corresponding sufficient statistics are presented in addition to a matrix formulation which facilitates an efficient calculation of the likelihood. This likelihood framework provides a consistent and unified framework with many standard models applied to capture–recapture–recovery data as special cases.