Constant width simultaneous confidence bands in multiple linear regression with predictor variables constrained in intervals

In the last fifty years, a great deal of research effort has been made on the construction of simultaneous confidence bands for a linear regression function. Two most frequently quoted confidence bands in the statistics literature are the Scheffé type and constant width bands over a given rectangular region of the predictor variables. For the constant width bands, a method is given by Gafarian [Gafarian, A.V., 1964, Confidence bands in straight line regression. Journal of the American Statistical Association, 59, 182–213.] for the calculation of critical constants only for the special case of one predictor variable. In this article, a method is proposed to construct constant width bands when there are any number of predictor variables. A new criterion for assessing a confidence band is also proposed; it is the probability that a confidence band excludes a false regression function and can be viewed as the power function of a test associated, naturally, with a confidence band. Under this criterion, a numerical comparison between the Scheffé type and constant width bands is then carried out. It emerges from this comparison that the constant width bands can be better than the Scheffé type bands for certain designs.     

Familles of confidence bands for the survival function under the general random censorship model and the Koziol-Green model

Familles of asymptotic 100(1 – α)% level confidence bands for the survival function under the general random right-censorship (GRC) model and the proportional-hazards model of random right-censorship, also known as the Koziol-Green (KG) model, are developed. The family of bands under the GRC model is based on the well-known product-limit estimator (PLE), and this family is rich in that it contains as special cases the bands of Hall and Wellner (1980) and Gillespie and Fisher (1979), and more generally, the GF-type and HW-type bands of Csörg? and Horváth (1986), as well as new bands not previously studied. The familles of bands under the KG model are based on the maximum-likelihood estimator of F under this particular model. We compare the PLE-based bands and the MLE-based bands under the KG model. This enables us to study the loss in efficiency of the former bands when used in a setting where they are not optimal. The notion of asymptotic relative width efficiency (ARWE), defined to be the limiting ratio of the sample sizes needed by the bands to achieve the same asymptotic widths, is employed to compare two bands. Through this efficiency measure it is shown that if the censoring parameter β is known, then the PLE-based bands are highly inefficient relative to the MLE-based bands when β is large. When β is not known, the MLE-based bands are asymptotically conservative. Despite their conservatism, they still dominate the PLE-based bands when β is not too small or equivalently when the degree of censoring is not too light. We also compare the various PLE-based bands under the GRC model. The resulting information is valuable for evaluating competing PLE-based bands. We illustrate the confidence bands by utilizing the well-known Channing House data.     

Likelihood Ratio Confidence Bands in Non–parametric Regression with Censored Data

Let ( X , Y ) be a random vector, where Y denotes the variable of interest possibly subject to random right censoring, and X is a covariate. We construct confidence intervals and bands for the conditional survival and quantile function of Y given X using a non-parametric likelihood ratio approach. This approach was introduced by Thomas & Grunkemeier (1975 ), who estimated confidence intervals of survival probabilities based on right censored data. The method is appealing for several reasons: it always produces intervals inside [0, 1], it does not involve variance estimation, and can produce asymmetric intervals. Asymptotic results for the confidence intervals and bands are obtained, as well as simulation results, in which the performance of the likelihood ratio intervals and bands is compared with that of the normal approximation method. We also propose a bandwidth selection procedure based on the bootstrap and apply the technique on a real data set.     

Simultaneous Confidence Bands for Linear Regression with Covariates Constrained in Intervals

Abstract. The focus of this article is on simultaneous confidence bands over a rectangular covariate region for a linear regression model with k>1 covariates, for which only conservative or approximate confidence bands are available in the statistical literature stretching back to Working & Hotelling (J. Amer. Statist. Assoc. 24 , 1929; 73–85). Formulas of simultaneous confidence levels of the hyperbolic and constant width bands are provided. These involve only a k‐dimensional integral; it is unlikely that the simultaneous confidence levels can be expressed as an integral of less than k‐dimension. These formulas allow the construction for the first time of exact hyperbolic and constant width confidence bands for at least a small k(>1) by using numerical quadrature. Comparison between the hyperbolic and constant width bands is then addressed under both the average width and minimum volume confidence set criteria. It is observed that the constant width band can be drastically less efficient than the hyperbolic band when k>1. Finally it is pointed out how the methods given in this article can be applied to more general regression models such as fixed‐effect or random‐effect generalized linear regression models.     

Exact Simultaneous Confidence Bands for Quadratic and Cubic Polynomial Regression with Applications in Dose Response Study

Exact simultaneous confidence bands (SCBs) for a polynomial regression model are available only in some special situations. In this paper, simultaneous confidence levels for both hyperbolic and constant width bands for a polynomial function over a given interval are expressed as multidimensional integrals. The dimension of these integrals is equal to the degree of the polynomial. Hence the values can be calculated quickly and accurately via numerical quadrature provided that the degree of the polynomial is small (e.g. 2 or 3). This allows the construction of exact SCBs for quadratic and cubic regression functions over any given interval and for any given design matrix. Quadratic and cubic regressions are frequently used to characterise dose response relationships in addition to many other applications. Comparison between the hyperbolic and constant width bands under both the average width and minimum volume confidence set criteria shows that the constant width band can be much less efficient than the hyperbolic band. For hyperbolic bands, comparison between the exact critical constant and conservative or approximate critical constants indicates that the exact critical constant can be substantially smaller than the conservative or approximate critical constants. Numerical examples from a dose response study are used to illustrate the methods.     

An alternative to confidence intervals constructed after a Hausman pretest in panel data

Consider panel data modelled by a linear random intercept model that includes a time‐varying covariate. Suppose that our aim is to construct a confidence interval for the slope parameter. Commonly, a Hausman pretest is used to decide whether this confidence interval is constructed using the random effects model or the fixed effects model. This post‐model‐selection confidence interval has the attractive features that it (a) is relatively short when the random effects model is correct and (b) reduces to the confidence interval based on the fixed effects model when the data and the random effects model are highly discordant. However, this confidence interval has the drawbacks that (i) its endpoints are discontinuous functions of the data and (ii) its minimum coverage can be far below its nominal coverage probability. We construct a new confidence interval that possesses these attractive features, but does not suffer from these drawbacks. This new confidence interval provides an intermediate between the post‐model‐selection confidence interval and the confidence interval obtained by always using the fixed effects model. The endpoints of the new confidence interval are smooth functions of the Hausman test statistic, whereas the endpoints of the post‐model‐selection confidence interval are discontinuous functions of this statistic.     

Interval estimation of the largest reliability of k weibull populations

A procedure is given for obtaining a random width confidence interval for the largest reliability of k Weibull populations. The procedure does not identify the populations for which the reliability would be a maximum. The maximum likelihood estimators or the simplified linear estimators of the reliability based on type II censored data are used. The cases considered include unknown shape parameters being equal or unequal. Simultaneous confidence intervals for the k reliabilities are also obtained. Tables for the lower and upper limits in selected cases are constructed using Monte Carlo methods.     

New confidence interval estimator of the signal-to-noise ratio based on asymptotic sampling distribution

In this paper, the asymptotic distribution of the signal-to-noise ratio (SNR) is derived and a new confidence interval for the SNR is introduced. An evaluation of the performance of the new interval compared to Sharma and Krishna (S–K) (1994) confidence interval for the SNR using Monte Carlo simulations is conducted. Data were randomly generated from normal, log-normal, χ², Gamma, and Weibull distributions. Simulations revealed that the performance of S–K interval is totally dependent on the amount of noise introduced and that it has a constant width for a given sample size. The S–K interval performs poorly in four of the distributions unless the SNR is around one. It is recommended against using the S–K interval for data from log-normal distribution even with SNR = 1. Unlike the S–K interval which does not account for skewness and kurtosis of the distribution, the new confidence interval for the SNR outperforms S–K for all five distributions discussed, especially when SNR?? 2. The proposed ranked set sampling (RSS) instead of simple random sampling (SRS) has improved the performance of both intervals as measured by coverage probability.     

Families of smooth confidence bands for the survival function under the general random censorship model

Randomly right censored data often arise in industrial life testing and clinical trials. Several authors have proposed asymptotic confidence bands for the survival function when data are randomly censored on the right. All of these bands are based on the empirical estimator of the survival function. In this paper, families of asymptotic (1-)100% level confidence bands are developed from the smoothed estimate of the survival function under the general random censorship model. The new bands are compared to empirical bands, and it is shown that for small sample sizes, the smooth bands have a higher coverage probability than the empirical counterparts.     

Sample size determination for the confidence interval of mean comparison adjusted by multiple covariates

The current method of determining sample size for confidence intervals does not accommodate multiple covariate adjustment. Under the normality assumption, the effect of multiple covariate adjustment on the standard error of the mean comparison is related to a Hotelling T ² statistic. Sample size can be calculated to obtain a desired probability of achieving a predetermined width in the confidence interval of the mean comparison with multiple covariate adjustment, given that the confidence interval includes the population parameter.     

Record-breaking data: a parametric comparison of the inverse-sampling and the random-sampling schemes

In many industrial quality control experiments and destructive stress testing, the only available data are successive minima (or maxima)i.e., record-breaking data. There are two sampling schemes used to collect record-breaking data: random sampling and inverse sampling. For random sampling, the total sample size is predetermined and the number of records is a random variable while in inverse-sampling the number of records to be observed is predetermined; thus the sample size is a random variable. The purpose of this papper is to determinevia simulations, which of the two schemes, if any, is more efficient. Since the two schemes are equivalent asymptotically, the simulations were carried out for small to moderate sized record-breaking samples. Simulated biases and mean square errors of the maximum likelihood estimators of the parameters using the two sampling schemes were compared. In general, it was found that if the estimators were well behaved, then there was no significant difference between the mean square errors of the estimates for the two schemes. However, for certain distributions described by both a shape and a scale parameter, random sampling led to estimators that were inconsistent. On the other hand, the estimated obtained from inverse sampling were always consistent. Moreover, for moderated sized record-breaking samples, the total sample size that needs to be observed is smaller for inverse sampling than for random sampling.     

Inferences on the intraclass correlation coefficients in the unbalanced two-way random effects model with interaction

In this paper, the hypothesis testing and interval estimation for the intraclass correlation coefficients are considered in a two-way random effects model with interaction. Two particular intraclass correlation coefficients are described in a reliability study. The tests and confidence intervals for the intraclass correlation coefficients are developed when the data are unbalanced. One approach is based on the generalized p-value and generalized confidence interval, the other is based on the modified large-sample idea. These two approaches simplify to the ones in Gilder et al. [2007. Confidence intervals on intraclass correlation coefficients in a balanced two-factor random design. J. Statist. Plann. Inference 137, 1199–1212] when the data are balanced. Furthermore, some statistical properties of the generalized confidence intervals are investigated. Finally, some simulation results to compare the performance of the modified large-sample approach with that of the generalized approach are reported. The simulation results indicate that the modified large-sample approach performs better than the generalized approach in the coverage probability and expected length of the confidence interval.     

Estimation on dependent right censoring scheme in an ordinary bivariate geometric distribution

Discrete lifetime data are very common in engineering and medical researches. In many cases the lifetime is censored at a random or predetermined time and we do not know the complete survival time. There are many situations that the lifetime variable could be dependent on the time of censoring. In this paper we propose the dependent right censoring scheme in discrete setup when the lifetime and censoring variables have a bivariate geometric distribution. We obtain the maximum likelihood estimators of the unknown parameters with their risks in closed forms. The Bayes estimators as well as the constrained Bayes estimates of the unknown parameters under the squared error loss function are also obtained. We considered an extension to the case where covariates are present along with the data. Finally we provided a simulation study and an illustrative example with a real data.     

Bayesian confidence intervals for means and variances of lognormal and bivariate lognormal distributions

The lognormal distribution is currently used extensively to describe the distribution of positive random variables. This is especially the case with data pertaining to occupational health and other biological data. One particular application of the data is statistical inference with regards to the mean of the data. Other authors, namely Zou et al. (2009), have proposed procedures involving the so-called “method of variance estimates recovery” (MOVER), while an alternative approach based on simulation is the so-called generalized confidence interval, discussed by Krishnamoorthy and Mathew (2003). In this paper we compare the performance of the MOVER-based confidence interval estimates and the generalized confidence interval procedure to coverage of credibility intervals obtained using Bayesian methodology using a variety of different prior distributions to estimate the appropriateness of each. An extensive simulation study is conducted to evaluate the coverage accuracy and interval width of the proposed methods. For the Bayesian approach both the equal-tail and highest posterior density (HPD) credibility intervals are presented. Various prior distributions (Independence Jeffreys' prior, Jeffreys'-Rule prior, namely, the square root of the determinant of the Fisher Information matrix, reference and probability-matching priors) are evaluated and compared to determine which give the best coverage with the most efficient interval width. The simulation studies show that the constructed Bayesian confidence intervals have satisfying coverage probabilities and in some cases outperform the MOVER and generalized confidence interval results. The Bayesian inference procedures (hypothesis tests and confidence intervals) are also extended to the difference between two lognormal means as well as to the case of zero-valued observations and confidence intervals for the lognormal variance. In the last section of this paper the bivariate lognormal distribution is discussed and Bayesian confidence intervals are obtained for the difference between two correlated lognormal means as well as for the ratio of lognormal variances, using nine different priors.     

Confidence bands for the laplace distribution

Based on a random sample from the Laplace population with unknown shape and scale parameters, one- and two-sided confidence bands on the entire cumulative distribution function and simultaneous confidence intervals for the interval probabilities under the distribution are constructed using Kolmogorov–Smirnov type statistics. Small sample and asymptotic percentiles of the relevant statistics are provided.     

A penalized spline estimator for fixed effects panel data models

Estimating nonlinear effects of continuous covariates by penalized splines is well established for regressions with cross-sectional data as well as for panel data regressions with random effects. Penalized splines are particularly advantageous since they enable both the estimation of unknown nonlinear covariate effects and inferential statements about these effects. The latter are based, for example, on simultaneous confidence bands that provide a simultaneous uncertainty assessment for the whole estimated functions. In this paper, we consider fixed effects panel data models instead of random effects specifications and develop a first-difference approach for the inclusion of penalized splines in this case. We take the resulting dependence structure into account and adapt the construction of simultaneous confidence bands accordingly. In addition, the penalized spline estimates as well as the confidence bands are also made available for derivatives of the estimated effects which are of considerable interest in many application areas. As an empirical illustration, we analyze the dynamics of life satisfaction over the life span based on data from the German Socio-Economic Panel. An open-source software implementation of our methods is available in the R package pamfe.     

Local regression when the responses are Interval-Censored

Conditional expectation imputation and local-likelihood methods are contrasted with a midpoint imputation method for bivariate regression involving interval-censored responses. Although the methods can be extended in principle to higher order polynomials, our focus is on the local constant case. Comparisons are based on simulations of data scattered about three target functions with normally distributed errors. Two censoring mechanisms are considered: the first is analogous to current-status data in which monitoring times occur according to a homogeneous Poisson process; the second is analogous to a coarsening mechanism such as would arise when the response values are binned. We find that, according to a pointwise MSE criterion, no method dominates any other when interval sizes are fixed, but when the intervals have a variable width, the local-likelihood method often performs better than the other methods, and midpoint imputation performs the worst. Several illustrative examples are presented.     

A Simultaneous Confidence Band for Dense Longitudinal Regression

We present a method of using local linear smoothing to construct simultaneous confidence bands for the mean function of densely spaced functional data. Our approach works well under mild conditions. In addition, the local linear estimator and its accompanying confidence band enjoy semiparametric efficiency in the sense that they are asymptotically equivalent to the counterparts obtained from the random trajectories entirely observed without errors. We illustrate the performance of the proposed confidence band through a simulation study. Furthermore, an application in food science is presented.     

基于Gibbs抽样算法的面板数据分位回归方法

 文章讨论了含有随机效应的面板数据模型,利用非对称Laplace分布与分位回归之间的关系,文章建立了一种贝叶斯分层分位回归模型。通过对非对称Laplace分布的分解,文章给出了Gibbs抽样算法下模型参数的点估计及区间估计,模拟结果显示,在处理含随机效应的面板数据模型中,特别是在误差非正态的情况下,本文的方法优于传统的均值模型方法。文章最后利用新方法对我国各地区经济与就业面板数据进行了实证研究,得到了有利于宏观调控的有用信息。     

Interval estimation for misclassification rate parameters in a complementary Poisson model

We investigate three interval estimators for binomial misclassification rates in a complementary Poisson model where the data are possibly misclassified: a Wald-based interval, a score-based interval, and an interval based on the profile log-likelihood statistic. We investigate the coverage and average width properties of these intervals via a simulation study. For small Poisson counts and small misclassification rates, the intervals can perform poorly in terms of coverage. The profile log-likelihood confidence interval (CI) is often proved to outperform the other intervals with good coverage and width properties. Lastly, we apply the CIs to a real data set involving traffic accident data that contain misclassified counts.