Change Points with Linear Trend Followed by Abrupt Change for the Exponential Distribution

We discuss testing procedures to detect if a random sequence of exponentially distributed random variables has been subjected to a linear trend change followed by an abrupt change. We propose three statistics and explore their distribution theories. As an illustration, we applied these tests to Stanford heart transplant data and airport inter arrival data.     

The quotient of two correlated normal variables with applications

The purpose of this paper is to present new results on the distribution of the quotient of two correlated normal variables and its use to solve many outstanding problems in all fields of science, economics, industry, commerce, education, and medicine. The results are relatively simple and may encourage one to use the exact distribution of the quotient. Examples in medicine, economics and industry will illustrate the applications. Further, these results open a new vista for the solution not only of problems involving bivariate distributions but also multivariate distributions.     

On testing hypotheses with divergence statistics

Using divergence measures based on entropy functions, a procedure to test statistical hypotheses is proposed. Replacing the parameters by suitable estimators in the expresion of the divergence measure, the test statistics are obtained. Asymptotic distributions for these statistics are given in several cases when maximum likelihood estimators are considered, so they can be used to construct confidence intervals and to test statistical hypotheses based on one or more samples. These results can also be applied to multinomial populations. Tests of goodness of fit and tests of homogeneity can be constructed.     

Description of a subfamily of the discrete pearson system as generalized-binomial distributions

In this paper we prove that a subfamily of distributions of the discrete Pearson, system, containing the Pólya distribution without replacement and hence the hypergeometric distribution, can be described as generalized-binomial distributions, i.e., the distribution of the number of successes which occur in independent trials. It is also shown that the probability of success will necessarily be different in each trial, with the exception of deterministic ones. As a consequence, all the properties of the generalized-binomial distribution will apply to this subfamily. Thus, applications to hypothesis testing and confidence intervals in the Pólya distribution are considered.     

A review of model selection procedures relevant to the weibull distribution

A number of goodness-of-fit and model selection procedures related to the Weibull distribution are reviewed. These procedures include probability plotting, correlation type goodness-of-fit tests, and chi-square goodness-of-fit tests. Also the Kolmogorow-Smirniv, Kuiper, and Cramer-Von Mises test statistics for completely specified hypothesis based on censored data are reviewed, and these test statistics based on complete samples for the unspecified parameters case are considered. Goodness-of-fit tests based on sample spacings, and a goodness-of-fit test for the Weibull process, is also discussed.

Model selection procedures for selecting between a Weibull and gamma model, a Weibull and lognormal model, and for selecting from among all three models are considered. Also tests of exponential versus Weibull and Weibull versus generalized gamma are mentioned.     

Improved Prediction Intervals and Distribution Functions

Abstract.  The plug-in solution is usually not entirely adequate for computing prediction intervals, as their coverage probability may differ substantially from the nominal value. Prediction intervals with improved coverage probability can be defined by adjusting the plug-in ones, using rather complicated asymptotic procedures or suitable simulation techniques. Other approaches are based on the concept of predictive likelihood for a future random variable. The contribution of this paper is the definition of a relatively simple predictive distribution function giving improved prediction intervals. This distribution function is specified as a first-order unbiased modification of the plug-in predictive distribution function based on the constrained maximum likelihood estimator. Applications of the results to the Gaussian and the generalized extreme-value distributions are presented.     

Conditional interval estimation of the mean following rejection of a two sided test

The effect of rejecting a two-sided preliminary test of significance for the mean of a normal distribution upon subsequent interval estimation of the mean is examined. For the case where the variance is known, conditional confidence intervals may be shorter than unconditional intervals, in contrast to the one-sided preliminary test case examined by Meeks and D’Agostino (1983, The American Statistician, 7, 134-136) . For the case where the variance is unknown and must be estimated by the sample variance, it is shown that customary intervals do not offer uniformly greater or lesser coverage than the nominal level.     

近单位根过程脉冲响应函数的置信区间

本文提出一个构造近单位根自回归过程脉冲响应函数的置信区间的新方法。新方法首先修正自回归系数估计的偏误,然后利用标准自举方法构造脉冲响应函数的置信区间。蒙特卡罗模拟结果表明在小样本时新方法的表现要明显优于已有的方法。新方法的实际覆盖率能够一致地达到或超过名义置信水平。     

Confidence intervals based on the deviance statistic for the hyperparameters in state space models

The main objective of this work is to evaluate the performance of confidence intervals, built using the deviance statistic, for the hyperparameters of state space models. The first procedure is a marginal approximation to confidence regions, based on the likelihood test, and the second one is based on the signed root deviance profile. Those methods are computationally efficient and are not affected by problems such as intervals with limits outside the parameter space, which can be the case when the focus is on the variances of the errors. The procedures are compared to the usual approaches existing in the literature, which includes the method based on the asymptotic distribution of the maximum likelihood estimator, as well as bootstrap confidence intervals. The comparison is performed via a Monte Carlo study, in order to establish empirically the advantages and disadvantages of each method. The results show that the methods based on the deviance statistic possess a better coverage rate than the asymptotic and bootstrap procedures.     

Characterization of the geometric distribution using properties of order statistics

Using certain properties of order statistics, the geometric distribution has been characterized when the components are independent and identically distributed. When the components are independent, the geometric distribution has been characterized in the class of either IFR or DFR discrete distributions. In particular, Ferguson's (1967) characterization theorem for independent components in a sample of size two has been extended in several directions.     

Small area housing stress estimation in Australia: Calculating confidence intervals for a spatial microsimulation model

This study provides small area housing stress estimates by tenure type in Australia with a way of calculating confidence intervals for a spatial microsimulation model. Findings reveal that prevalence of housing stress for private-renter, buyer, public-renter and owner households are 59.6%, 33.2%, 6.9%, and 0.3%, respectively. Almost two-thirds of these households are located in statistical local areas (SLAs) in eight capital cities, and a large number of them are in Sydney and Melbourne. Estimates for private renters and buyers are significantly high in some capitals and southeast coastal regions. About 95.7% of SLAs show accurate estimates with narrow confidence intervals.     

Antipodally symmetric distributions for orientation statistics

The conventional antipodally symmetric Bingham matrix distribution on the Stiefel manifold is generalized. Large sample maximum likelihood estimation and uniformity tests are discussed, and a parametric model for axial orientations (X-shapes) is suggested. A generalization of the Khatri-Mardia matrix distribution is developed to provide a model suitable for hybrids (T-shapes).     

Exact permutation inference for two sample repeated measures data

Experiments in which very few units are measured many times sometimes present particular difficulties. Interest often centers on simple location shifts between two treatment groups, but appropriate modeling of the error distribution can be challenging. For example, normality may be difficult to verify, or a single transformation stabilizing variance or improving normality for all units and all measurements may not exist. We propose an analysis of two sample repeated measures data based on the permutation distribution of units. This provides a distribution free alternative to standard analyses. The analysis includes testing, estimation and confidence intervals. By assuming a certain structure in the location shift model, the dimension of the problem is reduced by analyzing linear combinations of the marginal statistics. Recently proposed algorithms for computation of two sample permutation distributions, require only a few seconds for experiments having as many as 100 units and any number of repeated measures. The test has high asymptotic efficiency and good power with respect to tests based on the normal distribution. Since the computational burden is minimal, approximation of the permutation distribution is unnecessary.     

A comparison of estimators of the precision of two instruments

Three different estimators of the variance in errors of measurement, or the “imprecisions” of instruments, have been proposed (Hahn and Nelson, 1970; Maxwell, 1974; Hanumara: 1975) for the case where two measurements by one instrument and one measurement by a second instrument are available for each item. The estimators are compared by studying the probabilities of negative values, the biases, the variances, and the mean square errors. Based on these results, Hahn and Nelson's estimators for the variances in errors of measurement are recommended. Simultaneous confidence intervals of the imprecisions of instruments are also given.     

Bayesian prediction of order statistics based on k-record values from exponential distribution

Bayesian prediction of order statistics as well as the mean of a future sample based on observed record values from an exponential distribution are discussed. Several Bayesian prediction intervals and point predictors are derived. Finally, some numerical computations are presented for illustrating all the proposed inferential procedures.     

Estimation and prediction for the generalized inverted exponential distribution based on progressively first-failure-censored data with application

In this paper, the estimation of parameters for a generalized inverted exponential distribution based on the progressively first-failure type-II right-censored sample is studied. An expectation–maximization (EM) algorithm is developed to obtain maximum likelihood estimates of unknown parameters as well as reliability and hazard functions. Using the missing value principle, the Fisher information matrix has been obtained for constructing asymptotic confidence intervals. An exact interval and an exact confidence region for the parameters are also constructed. Bayesian procedures based on Markov Chain Monte Carlo methods have been developed to approximate the posterior distribution of the parameters of interest and in addition to deduce the corresponding credible intervals. The performances of the maximum likelihood and Bayes estimators are compared in terms of their mean-squared errors through the simulation study. Furthermore, Bayes two-sample point and interval predictors are obtained when the future sample is ordinary order statistics. The squared error, linear-exponential and general entropy loss functions have been considered for obtaining the Bayes estimators and predictors. To illustrate the discussed procedures, a set of real data is analyzed.     

Simultaneous confidence intervals by iteratively adjusted alpha for relative effects in the one-way layout

A bootstrap based method to construct 1−α simultaneous confidence intervals for relative effects in the one-way layout is presented. This procedure takes the stochastic correlation between the test statistics into account and results in narrower simultaneous confidence intervals than the application of the Bonferroni correction. Instead of using the bootstrap distribution of a maximum statistic, the coverage of the confidence intervals for the individual comparisons are adjusted iteratively until the overall confidence level is reached. Empirical coverage and power estimates of the introduced procedure for many-to-one comparisons are presented and compared with asymptotic procedures based on the multivariate normal distribution.     

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.     

Effects of some design factors on the distribution of similarity indices in cluster analysis

This article investigates the effects of number of clusters, cluster size, and correction for chance agreement on the distribution of two similarity indices, namely, Jaccard and Rand indices. Skewness and kurtosis are calculated for the two indices and their corrected forms then compared with those of the normal distribution. Three clustering algorithms are implemented: complete linkage, Ward, and K-means. Data were randomly generated from bivariate normal distributions with specified means and variance covariance matrices. Three-way ANOVA is performed to assess the significance of the design factors using skewness and kurtosis of the indices as responses. Test statistics for testing skewness and kurtosis and observed power are calculated. Simulation results showed that independent of the clustering algorithms or the similarity indices used, the interaction effect cluster size x number of clusters and the main effects of cluster size and number of clusters were found always significant for skewness and kurtosis. The three way interaction of cluster size x correction x number of clusters was significant for skewness of Rand and Jaccard indices using all clustering algorithms, but was not significant using Ward's method for both Rand and Jaccard indices, while significant for Jaccard only using complete linkage and K-means algorithms. The correction for chance agreement was significant for skewness and kurtosis using Rand and Jaccard indices when complete linkage method is used. Hence, such design factors must be taken into consideration when studying distribution of such indices.     

The folded t distribution

Measurements are frequently recorded without their algebraicsign. As a consequence the underlying distribution of measurements is replaced by a distribution of absolute measurements. When the underlying distribution is t the resulting distribution is called the “folded-t distribution”. Here we study this distribution, we find the relationship between the folded-t distribution and a special case of the folded normal distribution and we derive relationships of the folded-t distribution to other distributions pertaining to computer generation. Also tables are presented which give areas of the folded-t distribution.