Asymptotic bootstrap validity for finite markov chains

Two methods of bootstrap, viz., standard, and conditional, are presented for estimating the transition probabilities of a finite state Markov chain. Asymptotic validity of the bootstrap estimates are established for both methods. An applica- tion to a bootstrapped statistic for testing independence is briefly discussed together with some simulation results.     

Likelihood ratio test for testing order of continuous time finite markov chains

Anderson and Goodman ( 1957) have obtained the likelihood ratio tests and chi-square tests for testing the hypothesis about the order of discrete time finite Markov chains, On the similar lines we have obtained likeli¬hood ratio tests and chi-square tests (asymptotic) for testing hypotheses about the order of continuous time Markov chains (MC) with finite state space.     

The group inverse of finite homogeneous QBD processes

Generalized inverses of I?P, where P is a stochastic matrix, play an important role in the theory of Markov chains. In particular, the group inverse (I?P)^# has a probabilistic interpretation and is well suited for algorithmic implementation. We determine (I?P)^# for finite homogeneous quasi-birth-and-death (QBD) processes by exploiting both the structure of the process and the probabilistic properties of the group inverse.     

Development of Specific hypothesis tests for estimated markov chains

We develop and evaluate the validity and power of two specific tests for the transition probabilities in a Markov chain estimated from aggregate frequency data. The two null hypotheses considered are (1) constancy of the diagonal elements of the one-step transition probability matrix and (2) an arbitrarily chosen transition probability’s being equal to a specific value. The formation of tests uses a general framework for statistical inference on estimated Markov processes; we also indicate how this framework can be used to form tests for a variety of other hypotheses. The validity and power performance of the two tests formed in this paper are examined in factorially designed Monte Carlo experiments. The results indicate that the proposed tests lead to type I error probabilities which are close to the desired levels and to high power against even small deviations from the null hypotheses considered.     

The Equivalence of Two Rules of Classification for Two Populations

In many applied classification problems, the populations of interest are defined in terms of ranges for the dependent variable. In these situations, it is intuitively appealing to classify individuals into the respective populations based on their estimated conditional expectation. On the other hand, based on theoretical considerations, one may wish to use the classification rule based on the posterior probabilities. This article shows that under certain conditions these two classification rules are equivalent.     

Estimation of finite population variance in successive sampling using multi-auxiliary variables

ABSTRACT

In this paper, a general class of estimators for estimating the finite population variance in successive sampling on two occasions using multi-auxiliary variables has been proposed. The expression of variance has also been derived. Further, it has been shown that the proposed general class of estimators is more efficient than the usual variance estimator and the class of variance estimators proposed by Singh et al. (2011) when we used more than one auxiliary variable. In addition, we support this with the aid of numerical illustration.     

Comparison of regression estimators using Pitman measures of nearness

Received: August 5, 1999; revised version: June 14, 2000     

Constructing an unbiased estimator of population mean in finite populations using auxiliary information

The objective of this paper is to construct an unbiased estimator (up to order 0(1/n)) of the population mean of the study variatey which is more efficient than the sample mean of the ‘n’ obsrvedy-values. In particular, the unbiased estimators are discussed for the cases of positive and negative correlations of the study variatey and the auxiliary variatex.     

Confidence limits for estimates of totals from stratified samples, with application to medicare Part B overpayment audits

Superpopulation models are proposed that should be appropriate for modelling sample-based audits of Medicare payments and other overpayment situations. Simulations are used to estimate the coverage probabilities of confidence intervals formed using the standard Stratified Expansion and Combined Ratio estimators of the total. Despite severe departures from the usual model of normal deviations, these methods have actual coverage probabilities reasonably close to the nominal level specified by the US government's sampling guidelines. An exception occurs when all claims from a single sampling unit are either completely allowed, or completely denied, and for this situation an alternative is explored. A balanced sampling design is also examined, but shown to make no improvement over ordinary stratified samples used in conjunction with ratio estimates.     

A bayesian method in determining the order of a finite state markov chain

In this paper, we use the Bayesian method in the application of hypothesis testing and model selection to determine the order of a Markov chain. The criteria used are based on Bayes factors with noninformative priors. Com¬parisons with the commonly used AIC and BIC criteria are made through an example and computer simulations. The results show that the proposed method is better than the AIC and BIC criteria, especially for Markov chains with higher orders and larger state spaces.     

Fitting finite mixture models using iterative Monte Carlo classification

Parameters of a finite mixture model are often estimated by the expectation–maximization (EM) algorithm where the observed data log-likelihood function is maximized. This paper proposes an alternative approach for fitting finite mixture models. Our method, called the iterative Monte Carlo classification (IMCC), is also an iterative fitting procedure. Within each iteration, it first estimates the membership probabilities for each data point, namely the conditional probability of a data point belonging to a particular mixing component given that the data point value is obtained, it then classifies each data point into a component distribution using the estimated conditional probabilities and the Monte Carlo method. It finally updates the parameters of each component distribution based on the classified data. Simulation studies were conducted to compare IMCC with some other algorithms for fitting mixture normal, and mixture t, densities.     

Approximation of distributions by using the Anderson Darling statistic

ABSTRACT

In practice, it is often not possible to find an appropriate family of distributions which can be used for fitting the sample distribution with high precision. In these cases, it seems to be opportune to search for the best approximation by a family of distributions instead of an exact fit. In this paper, we consider the Anderson–Darling statistic with plugged-in minimum distance estimator for the parameter vector. We prove asymptotic normality of the Anderson–Darling statistic which is used for a test of goodness of approximation. Moreover, we introduce a measure of discrepancy between the sample distribution and the model class.     

Equivariant estimation of a normal mean using a normal concomitant variable for covariance adjustment

Equivariant point estimators of one component of a bivariate normal mean vector are considered when the second component is known. Equivariant point estimators are characterized and compared in terms of their risk functions with respect to a normalized squared-error loss function. Specific point estimators that dominate the usual estimator when the squared correlation coefficient is sufficiently large are provided.     

The computation of standard normal distribution integral in any required precision based on reliability method

Based on reliability theory, the value of the standard normal distribution integral can be obtained by calculating the probability of the failure domain of the linear performance function. After the sample space is divided into some sub-sample spaces, a number of sub-failure domains are obtained. In the paper, the methods of computing the probabilities of sub-failure domains are discussed. All the formulae and the steps of computing the standard normal distribution integral which meet any required precision are given in the paper. Examples show that it is easy for the method to compute the standard normal distribution integral.     

Correcting the bias due to dependent censoring of the survival estimator by conditioning

We consider the conditional estimation of the survival function of the time T₂ to a second event as a function of the time T₁ to a first event when there is a censoring mechanism acting on their sum T₁+T₂. The problem has been motivated by a treatment interruption study aimed at improving the quality of life of HIV-infected patients. We base the analysis on the survival function of T₂ given that T₁∈I, where I represents a period of scientific interest (1 trimester, 1 year, 2 years, etc.) and propose a non-parametric estimator for the survival function of T₂ given that T₁∈I, which takes into account both the selection bias and the heterogeneity due to the dependent censoring. The proposed estimator for the survival function uses the risk group of T₂ conditioned on the categories of T₁ and corrects for the dependent censoring using weights defined by the observed values of T₁. The estimator, properly normalized, converges weakly to a zero-mean Gaussian process. We estimate the variance of the limiting process via a bootstrap methodology. Properties of the proposed estimator are illustrated by an extensive simulation study. The motivating data set is analysed by means of this new methodology.     

Estimators of contingent probabilities and means with actuarial applications

Contingent probabilities and means are ubiquitous in actuarial science. The correct interpretation of contingent probabilities and means as well as the probability theory behind them have been addressed by researchers. In this article, we explore their statistical aspect. We give non-parametric estimators of contingent probabilities and means. Then, we show that our estimators are strongly consistent. Moreover, we give the asymptotic distributions of our estimators. Finally, we provide several examples to demonstrate the applications of these estimators in actuarial science.     

On the Non-attainability of Chebyshev Bounds

In this department The American Statistician publishes articles, reviews, and notes of interest to teachers of the first mathematical statistics course and of applied statistics courses. The department includes the Accent on Teaching Materials section; suitable contents for the section are described under the section heading. Articles and notes for the department, but not intended specifically for the section, should be useful to a substantial number of teachers of the indicated types of courses or should have the potential for fundamentally affecting the way in which a course is taught.     

On the finite sample performance of exogeneity tests of revankar, revankar and hartley and wu-hausman

In this paper we compare the performance of the exogeneity tests of Revankar, Revankar and Hartley and Wu-Hausman for the cases of two and three included endogenous variables. The distribution and power functions are evaluated using the conditional distributions given in Kariya and Hodoshima. Our results indicate that the Revankar's test is the most powerful for large values of the concentration parameter and the Revankar and Hartley test is the most powerful for small values of the concentration parameter.     

Modelling accelerated life test data by using a Bayesian approach

Summary. Because of the high reliability of many modern products, accelerated life tests are becoming widely used to obtain timely information about their time-to-failure distributions. We propose a general class of accelerated life testing models which are motivated by the actual failure process of units from a limited failure population with a positive probability of not failing during the technological lifetime. We demonstrate a Bayesian approach to this problem, using a new class of models with non-monotone hazard rates, the hazard model with potential scope for use far beyond accelerated life testing. Our methods are illustrated with the modelling and analysis of a data set on lifetimes of printed circuit boards under humidity accelerated life testing.     

On finite sample properties of nonparametric discrete asymmetric kernel estimators

The discrete kernel method was developed to estimate count data distributions, distinguishing discrete associated kernels based on their asymptotic behaviour. This study investigates the class of discrete asymmetric kernels and their resulting non-consistent estimators, but this theoretical drawback of the estimators is balanced by some interesting features in small/medium samples. The role of modal probability and variance of discrete asymmetric kernels is highlighted to help better understand the performance of these kernels, in particular how the binomial kernel outperforms other asymmetric kernels. The performance of discrete asymmetric kernel estimators of probability mass functions is illustrated using simulations, in addition to applications to real data sets.