A modification of the partial sequential probability ratio test

For the problem of discriminating between two simple hypoth¬eses concerning a Koopman - Darmois parameter, a modification of the partial sequential probability ratio test is proposed where instead of drawing only one fixed sample, two fixed samples are drawn and then Wald's SPRT is started. The OC and the ASN func¬tions are derived. Numerical comparisons are made with Wald's and Read's procedures for testing the normal mean with known variance. For some parameter values, the test procedure has a lower ASN than that of Read's procedure.     

Minimax sequential tests of composite hypotheses on the drift of a Wiener process

A Wiener process with unknown drift parameter μ is, beginning at O, observed continuously and one has to decide between the hypotheses μ≤0 and μ>0. For loss functions of the form sμ^r and linear cost functions one wants to determine a minimax sequential test. Generalizing the results of DeGroot (1960) a minimax test in the class of all symmetrical SPRT’s is given in an explicit form. On the other hand it is shown that this SPRT is, in general, no longer minimax in the class of all sequential tests.     

A sequential sampling plan for the inverse gaussian mean

The inverse gaussian distribution is a flexible model which has been extensively applied in the theory of generalized linear models and accelerated life testing where early failure times predominate. More recently it has received attention in areas such as quality control, and as an underlying model that provides an alternative to the analysis of variance. In reliability testing and acceptance sampling data acquisition is often in the face of scarce resources and may be both costly and time-consuming. In such settings it is desirable to reach a statistically sound decision as quickly as possible. Based on sequential probability ratio tests (SPRT), sequential sampling plans provide one method of arriving at a timely, statistically based decision. A sequential sampling plan for the inverse gaussian process mean when the value of the shape parameter of the density is known is presented in this paper.     

On the linear combinations of the means of inverse gaussian distributions

Consider the two parameter Inverse Gaussian distribution with mean μ and scale parameter λ.

Suppose one is interested in testing a problem on a linear combination for the means of Inverse Gaussian distributions. For this problem a test and confidence intervals are proposed when: (1) λ’s are known and; (2) λ’s are unknown.

Finally an application of the procedures is illustrated with a data set of failure times of high-speed turbine bearings.     

An R package for testing goodness of fit: goft

This R package implements three types of goodness-of-fit tests for some widely used probability distributions where there are unknown parameters, namely tests based on data transformations, on the ratio of two estimators of a dispersion parameter, and correlation tests. Most of the considered tests have been proved to be powerful against a wide range of alternatives and some new ones are proposed here. The package's functionality is illustrated with several examples by using some data sets from the areas of environmental studies, biology and finance, among others.     

Tests for symmetry of two-piece normal distribution

In this paper, classical optimum tests for symmetry of two-piece normal distribution is derived. Uniformly most powerful one-sided test for the skewness parameter is obtained when the location and scale parameters are known and is compared with sequential probability ratio test. An ad-hoc test for symmetry and likelihood ratio test for symmetry for large samples, can be found in literature for this distribution. But in this paper, we derive exact likelihood ratio test for symmetry, when location parameter is known. The exact power of the test is evaluated for different sample sizes.     

Modified kolmogorov-smirnov tests of goodness of fit

The Kolmogorov-Smirnov (K–S) one-sided and two-sided tests of goodness of fit based on the test statistics D⁺ _n D^? _n and D_n are equivalent to tests based on taking the cumulative probability of the i–th order statistic of a sample of size n to be (i–.5)/n. Modified test statistics C⁺ _n, C^? _n and C_n are obtained by taking the cumulative probability to be i/(n+l). More generally, the cumula-tive probability may be taken to be (i?δ)/(n+l?2δ), as suggested by Blom (1958), where 0 less than or equal δ less than or equal .5. Critical values of the test statis-tics can be found by interpolating inversely in tables of the proba-bility integrals obtained by setting a=l/(n+l?2δ) in an expression given by Pyke (1959). Critical values for the D's (corresponding to δ=.5) have been tabulated to 5DP by Miller (1956) for n=1(1)100. The authors have made analogous tabulations for the C's (corresponding to δ=0) [previously tabulated by Durbin (1969) for n=1(1)60(2)100] and for the test statistics E⁺ _n, E^? _n and E_n corresponding to δ f.3. They have also made a Monte Carlo comparison of the power of the modified tests with that of the K–S test for several hypothetical distributions. In a number of cases, the power of the modified tests is greater than that of the K–S test, especially when the standard deviation is greater under the alternative than under the null hypo-thesis.     

A Gaussian alternative to using improper confidence intervals

The problem posed by exact confidence intervals (CIs) which can be either all-inclusive or empty for a nonnegligible set of sample points is known to have no solution within CI theory. Confidence belts causing improper CIs can be modified by using margins of error from the renewed theory of errors initiated by J. W. Tukey—briefly described in the article—for which an extended Fraser's frequency interpretation is given. This approach is consistent with Kolmogorov's axiomatization of probability, in which a probability and an error measure obey the same axioms, although the connotation of the two words is different. An algorithm capable of producing a margin of error for any parameter derived from the five parameters of the bivariate normal distribution is provided. Margins of error correcting Fieller's CIs for a ratio of means are obtained, as are margins of error replacing Jolicoeur's CIs for the slope of the major axis. Margins of error using Dempster's conditioning that can correct optimal, but improper, CIs for the noncentrality parameter of a noncentral chi-square distribution are also given.     

A generalization of anderson's modified sequential probability ratio test

A generalization of Anderson's sequential probability ratio test procedure is proposed in which the continuation region is bounded by a pair of converging lines up to a certain stage of the experiment and later by another pair of converging lines until the procedure is truncated at a predetermined stage of the experiment. The OC and the ASN functions have been derived. For certain parameter values the proposed procedure attains lower average sample numbers than that attainable by any other known procedure.     

Bayesian Variable Selection in Markov Mixture Models

Monte Carlo simulation is used to evaluate the actual confidence levels of five different approximations for confidence intervals for the probability of success in Markov dependent trials. The approximations involve the conditional probability of success as a nuisance parameter, and the effects of substituting Klotz's (1973), Price's (1976), and a new estimator are also evaluated. The new estimator is less biased and tends to increase the confidence level. A program for calculating the estimator and the confidence interval approximations is available.     

Statistical inference for the extended generalized inverse Gaussian model

This paper deals with the extended generalized inverse Gaussian (EGIG) distribution which has more than one turning point of the failure rate for certain values of the parameters. The EGIG model is a versatile model for analysing lifetime data and has one additional parameter, δ, than the GIG model's three parameters [B. Jorgensen, Statistical Properties of the Generalized Inverse Gaussian Distribution, Springer-Verlag, New York, 1982]. For the EGIG model, the maximum-likelihood estimation of the four parameters is discussed and a score test is developed for testing the importance of the additional parameter, δ. A non-central chi-square approximation to the power of the score test is provided. Simulation studies are carried out to examine the performance of the score test and the Wald confidence intervals. Finally, an example discussed by Jorgensen [5] is provided to illustrate that the EGIG model fits the data better than the GIG of Jorgensen [5]. Three other examples are presented and the power comparisons are displayed for each.     

Two stochastic models for the geeta distribution

The first stochastic model is based upon two urns A and B, where A contains a fixed number of white and black balls and B is empty. The player selects an integer β ≥, 2 and draws the balls one by one (with replacement) from urn A and balls of the same colour are put in urn B. The process is continued as long as the number of white balls in B exceeds (β-1) times the number of black balls in B. The player stops after drawing β(x-1) balls and is declared to be a winner if urn B has (x-1) black balls. It is shown that x has the Geeta distribution.

Assuming that the mean μ is a function of two parameters θ and β it has been shown that for small changes inthe value of θ there exists a difference-differential equation which leads to the Geeta distribution.     

Moments of the decisive sample number and of the average time to termination of sequential tests

The moments of the decisive sample number and of the average time to termination of sequential tests are derived. All moments depend only on the probability of continuation of the sequential test at each step, or at each instant of time for time processes. Two theorems are proved. Although the second theorem may be proved by quoting a result used in Barlow and Proschan (1965), p. 33, or both theorems from a lemma in Govindarajulu (1975), pp. 518-519, the author’s direct proofs illuminate the results.     

Stepwise likelihood ratio statistics in sequential studies

Summary.  It is well known that in a sequential study the probability that the likelihood ratio for a simple alternative hypothesis H ₁ versus a simple null hypothesis H ₀ will ever be greater than a positive constant c will not exceed 1/ c under H ₀. However, for a composite alternative hypothesis, this bound of 1/ c will no longer hold when a generalized likelihood ratio statistic is used. We consider a stepwise likelihood ratio statistic which, for each new observation, is updated by cumulatively multiplying the ratio of the conditional likelihoods for the composite alternative hypothesis evaluated at an estimate of the parameter obtained from the preceding observations versus the simple null hypothesis. We show that, under the null hypothesis, the probability that this stepwise likelihood ratio will ever be greater than c will not exceed 1/ c . In contrast, under the composite alternative hypothesis, this ratio will generally converge in probability to ∞. These results suggest that a stepwise likelihood ratio statistic can be useful in a sequential study for testing a composite alternative versus a simple null hypothesis. For illustration, we conduct two simulation studies, one for a normal response and one for an exponential response, to compare the performance of a sequential test based on a stepwise likelihood ratio statistic with a constant boundary versus some existing approaches.     

Goodness-of-Fit of the Distribution of Time-to-First-Occurrence in Recurrent Event Models

Imperfect repair models are a class of stochastic models that deal with recurrent phenomena. This article focuses on the Block, Borges, and Savits (1985) age-dependent minimal repair model (the BBS model) in which a system that fails at time t undergoes one of two types of repair: with probability p(t), a perfect repair is performed, or with probability 1-p(t), a minimal repair is performed. The goodness-of-fit problem of interest concerns the initial distribution of the failure ages. In particular, interest is on testing the null hypothesis that the hazard rate function of the time-to-first-event-occurrence, λ(·), is equal to a prespecified hazard rate function λ₀(·). This paper extends the class of hazard-based smooth goodness-of-fit tests introduced in Peña (1998a) to the case where data accrual is from a BBS model. The goodness-of-fit tests are score tests derived by reformulating Neyman's idea of smooth tests in terms of hazard functions. Omnibus as well as directional tests are developed and simulation results are presented to illustrate the sensitivities of the proposed tests for certain types of alternatives.     

A two-stage sequential design and analysis for comparing binomial rates with possibly random dropouts

Many clinical trials involve two-stage sequential designs with one interim analysis (Elashoff and Reedy, 1984, Biometrics 41, 791-795.) In this paper we present a situation where events are counted only at two fixed calendar time points and some patients may dropout during the time intervals. In the two-stage case, naive application of Tsiatis’s (1984, JASA 77, 855-861) logrank and Wilcoxon tests, which are for continuous survival time, is shown to lead to conservative type-I error rates and lower power. The two-stage sequential boundaries can also be calculated directly, rather than by simulation as was done by DeMets and Gail (1985, Biometrics 41, 1039-1044) with the assumption of some survival models, and are shown to be more flexible than the Pocock (1977, Biometrika 64, 191-199) and O’Brien-Fleming (1983, Biometrics 35, 549-556) boundaries since the former do not require an assumption on the correlation of the test statistics for the two stages. Repeated confidence intervals are also discussed. The design and approach are motivated by clinical trials studying treatment effects on vertebral fracture rates in elderly osteoporotic women. An example (Tilyard, et al. New England Journal of Medicine, 1992) is given to illustrate the method.     

Two-sided tests for change in level for correlated data

The classical problem of change point is considered when the data are assumed to be correlated. The nuisance parameters in the model are the initial level μ and the common variance σ². The four cases, based on none, one, and both of the parameters are known are considered. Likelihood ratio tests are obtained for testing hypotheses regarding the change in level, δ, in each case. Following Henderson (1986), a Bayesian test is obtained for the two sided alternative. Under the Bayesian set up, a locally most powerful unbiased test is derived for the case μ=0 and σ²=1. The exact null distribution function of the Bayesian test statistic is given an integral representation. Methods to obtain exact and approximate critical values are indicated.     

On Some Aspects of Unbiased Estimation of Parameters in Quasi-Binomial Distributions

In this article, an attempt has been made to settle the question of existence of unbiased estimator of the key parameter p of the quasi-binomial distributions of Type I (QBD I) and of Type II (QBD II), with/without any knowledge of the other parameter φ appearing in the expressions for probability functions of the QBD's. This is studied with reference to a single observation, a random sample of finite size m as also with samples drawn by suitably defined sequential sampling rules.     

Sequential procedures based on M-estimators with discontinuous score functions

Bounded-width sequential confidence intervals and sequential tests for regression parameter based on M-estimators are extended to the case where the score-functions generating the M-estimators have jump-discontinuities. In the context of the asymptotic normality of the stopping variable, for the confidence interval problem, it is observed that the jump-discontinuities induce a slower rate of convergence. The proofs of the main theorems rest on the weak convergence of some related processes and this is also studied.     

Simultaneous confidence sets and confidence intervals for multiple ratios

This paper deals with the problem of simultaneously estimating multiple ratios. In the simplest case of only one ratio parameter, Fieller's theorem (J. Roy. Statist. Soc. Ser. B 16 (1954) 175) provides a confidence interval for the single ratio. For multiple ratios, there is no method available to construct simultaneous confidence intervals that exactly satisfy a given familywise confidence level. Many of the methods in use are conservative since they are based on probability inequalities. In this paper, first we consider exact simultaneous confidence sets based on the multivariate t-distribution. Two approaches of determining the exact simultaneous confidence sets are outlined. Second, approximate simultaneous confidence intervals based on the multivariate t-distribution with estimated correlation matrix and a resampling approach are discussed. The methods are applied to ratios of linear combinations of the means in the one-way layout and ratios of parameter combinations in the general linear model. Extensive Monte Carlo simulation is carried out to compare the performance of the various methods with respect to the stability of the estimated critical points and of the coverage probabilities.