The unit-inverse Gaussian distribution: A new alternative to two-parameter distributions on the unit interval

A new two-parameter distribution over the unit interval, called the Unit-Inverse Gaussian distribution, is introduced and studied in detail. The proposed distribution shares many properties with other known distributions on the unit interval, such as Beta, Johnson S_B, Unit-Gamma, and Kumaraswamy distributions. Estimation of the parameters of the proposed distribution are obtained by transforming the data to the inverse Gaussian distribution. Unlike most distributions on the unit interval, the maximum likelihood or method of moments estimators of the parameters of the proposed distribution are expressed in simple closed forms which do not need iterative methods to compute. Application of the proposed distribution to a real data set shows better fit than many known two-parameter distributions on the unit interval.     

Estimating Transition Probabilities for Ignorable Intermittent Missing Data in a Discrete-Time Markov Chain

This article considers a discrete-time Markov chain for modeling transition probabilities when multiple successive observations are missing at random between two observed outcomes using three methods: a na\"?ve analog of complete-case analysis using the observed one-step transitions alone, a non data-augmentation method (NL) by solving nonlinear equations, and a data-augmentation method, the Expectation-Maximization (EM) algorithm. The explicit form of the conditional log-likelihood given the observed information as required by the E step is provided, and the iterative formula in the M step is expressed in a closed form. An empirical study was performed to examine the accuracy and precision of the estimates obtained in the three methods under ignorable missing mechanisms of missing completely at random and missing at random. A dataset from the mental health arena was used for illustration. It was found that both data-augmentation and nonaugmentation methods provide accurate and precise point estimation, and that the na\"?ve method resulted in estimates of the transition probabilities with similar bias but larger MSE. The NL method and the EM algorithm in general provide similar results whereas the latter provides conditional expected row margins leading to smaller standard errors.     

An iterative algorithm for solving sparse linear equations

This article introduces a new iterative technique for solving systems of linear equations of the kind Ax = b. Convergence, and with a given rate, is guaranteed with the square nonsingular matrix A being non-negative. The iterative algorithm depends on a scheme derived from Bayesian updating. The algorithm is shown to compare very favorably with the wisely used GMRES routine. With the algorithm being easy to code, it has the potential to be highly useable.     

Optimal window width choice in spectral density estimation

This paper deals with optimal window width choice in on-parametric lag or spectral window estimation of the spectral density of a stationary zero-mean process. Several approaches are reviewed: cross-validation-based methods as described by Hurvich(1985) BelträHo and Bloomfield (1987) and Hurvich and Belträo (1990); an iterative pro-cedure developed by Bühlmann (1996); and a bootstrap approach followed by Franke and Hardle (1992). These methods are compared in terms of the mean square error,the mean square percentage error, and a third measure of the istance between the true spectral density and its estimate. The comparison is based on a simulation study, the simulated processes being in he class of ARMA (5,5) processes. On the basis of simu-lation evidence we suggest to use a slightly modified version of Biihlmann's (1996)iterative method. This paper also makes a minor correction of the bootstrap criterion by Franke and Härdle (1992).     

Likelihood methods in randomized trials with noncompliance and subsequent nonresponse

This article considers likelihood methods for estimating the causal effect of treatment assignment for a two-armed randomized trial assuming all-or-none treatment noncompliance and allowing for subsequent nonresponse. We first derive the observed data likelihood function as a closed form expression of the parameter given the observed data where both response and compliance state are treated as variables with missing values. Then we describe an iterative procedure which maximizes the observed data likelihood function directly to compute a maximum likelihood estimator (MLE) of the causal effect of treatment assignment. Closed form expressions at each iterative step are provided. Finally we compare the MLE with an alternative estimator where the probability distribution of the compliance state is estimated independent of the response and its missingness mechanism. Our work indicates that direct maximum likelihood inference is straightforward for this problem. Extensive simulation studies are provided to examine the finite sample performance of the proposed methods.     

A Study on Reliability for A Two-Item Cold Standby Markov Repairable System with Neglected Failures

This article studies reliability for a Markov repairable two-item cold standby system with neglected failures. In the system, if a failed time of the system is too short (less than a given critical value) to cause the system to fail, then the failed time may be omitted from the downtime record, i.e., the failure effect could be neglected. In ion-channel modeling, this situation is called the time interval omission problem. The availability indices and the mean downtime are presented as two measures of reliability for this repairable system. Some numerical examples are shown to illustrate the results obtained in this article.     

Subexponential asymptotics of the waiting time distribution in a single-server queue with multiple Markovian arrival streams

This paper considers subexponential asymptotics of the tail distributions of waiting times in stationary work-conserving single-server queues with multiple Markovian arrival streams, where all arrival streams are modulated by the underlying Markov chain with finite states and service time distributions may differ for different arrival streams. Under the assumption that the equilibrium distribution of the overall (i.e., customer-average) service time distribution is subexponential, a subexponential asymptotic formula is first shown for the virtual waiting time distribution, using a closed formula recently found by the author. Further when customers are served on a FIFO basis, the actual waiting time and sojourn time distributions of customers from respective arrival streams are shown to have the same asymptotics as the virtual waiting time distribution.     

Approximate MLEs of the parameters of location-scale models under type II censoring

Log-location-scale distributions are widely used parametric models that have fundamental importance in both parametric and semiparametric frameworks. The likelihood equations based on a Type II censored sample from location-scale distributions do not provide explicit solutions for the para-meters. Statistical software is widely available and is based on iterative methods (such as, Newton Raphson Algorithm, EM algorithm etc.), which require starting values near the global maximum. There are also many situations that the specialized software does not handle. This paper provides a method for determining explicit estimators for the location and scale parameters by approximating the likelihood function, where the method does not require any starting values. The performance of the proposed approximate method for the Weibull distribution and Log-Logistic distributions is compared with those based on iterative methods through the use of simulation studies for a wide range of sample size and Type II censoring schemes. Here we also examine the probability coverages of the pivotal quantities based on asymptotic normality. In addition, two examples are given.     

Bayesian density estimation using bernstein polynomials

We propose a Bayesian nonparametric procedure for density estimation, for data in a closed, bounded interval, say [0,1]. To this aim, we use a prior based on Bemstein polynomials. This corresponds to expressing the density of the data as a mixture of given beta densities, with random weights and a random number of components. The density estimate is then obtained as the corresponding predictive density function. Comparison with classical and Bayesian kernel estimates is provided. The proposed procedure is illustrated in an example; an MCMC algorithm for approximating the estimate is also discussed.     

An extension of the von mises distribution

Directional distribution theory is very useful for the estimation of directional spectra needed for the analysis of time series data. A four parameter directional exponential family is discussed. Depending on the values of its parameters this distribution can be unimodal symmetric, bimodal symmetric, unimodal non-symmetric, or bimodal non-symmetric. The moments of this distribution are found, and equations leading to maximum-likelihood estimates of the parameters along with an outline on numerical procedures for solving these equations are given. FORTRAN subroutines implementing these procedures are available from the authors. Finally, some applications of the new directional density are given.     

Bayesian inference and prediction of order statistics for a Type-II censored Weibull distribution

This paper describes the Bayesian inference and prediction of the two-parameter Weibull distribution when the data are Type-II censored data. The aim of this paper is twofold. First we consider the Bayesian inference of the unknown parameters under different loss functions. The Bayes estimates cannot be obtained in closed form. We use Gibbs sampling procedure to draw Markov Chain Monte Carlo (MCMC) samples and it has been used to compute the Bayes estimates and also to construct symmetric credible intervals. Further we consider the Bayes prediction of the future order statistics based on the observed sample. We consider the posterior predictive density of the future observations and also construct a predictive interval with a given coverage probability. Monte Carlo simulations are performed to compare different methods and one data analysis is performed for illustration purposes.     

Robust non-parametric smoothing of non-stationary time series

Motivated by the need of extracting local trends and low frequency components in non-stationary time series, this paper discusses methods of robust non-parametric smoothing. Basic approach is the combination of the parametric M-estimation with kernel and local polynomial regression methods. The result is an iterative estimator that retains a linear structure, but has kernel weights also in the direction of the prediction errors. The design of smoothing coefficients is carried out with robust cross-validation criteria and rules of thumb. The method works well both to remove the influence of patches of outliers and to detect the local breaks and persistent structural change in time series.     

Regression models for analyzing clustered binary and continuous outcomes under an assumption of exchangeability

Scientific experiments commonly result in clustered discrete and continuous data. Existing methods for analyzing such data include the use of quasi-likelihood procedures and generalized estimating equations to estimate marginal mean response parameters. In applications to areas such as developmental toxicity studies, where discrete and continuous measurements are recorded on each fetus, or clinical ophthalmologic trials, where different types of observations are made on each eye, the assumption that data within cluster are exchangeable is often very reasonable. We use this assumption to formulate fully parametric regression models for clusters of bivariate data with binary and continuous components. The regression models proposed have marginal interpretations and reproducible model structures. Tractable expressions for likelihood equations are derived and iterative schemes are given for computing efficient estimates (MLEs) of the marginal mean, correlations, variances and higher moments. We demonstrate the use the ‘exchangeable’ procedure with an application to a developmental toxicity study involving fetal weight and malformation data.     

Tests for Regression Models Fitted to Survey Data

Data from complex surveys are being used increasingly to build the same sort of explanatory and predictive models as those used in the rest of statistics. Unfortunately the assumptions underlying standard statistical methods are not even approximately valid for most survey data. The problem of parameter estimation has been largely solved, at least for routine data analysis, through the use of weighted estimating equations, and software for most standard analytical procedures is now available in the major statistical packages. One notable omission from standard software is an analogue of the likelihood ratio test. An exception is the Rao–Scott test for loglinear models in contingency tables. In this paper we show how the Rao–Scott test can be extended to handle arbitrary regression models. We illustrate the process of fitting a model to survey data with an example from NHANES.     

Weibull inference using trimmed samples and prior information

Trimmed samples are commonly used in several branches of statistical methodology, especially when the presence of contaminated data is suspected. Assuming that certain proportions of the smallest and largest observations from a Weibull sample are unknown or have been eliminated, a Bayesian approach to point and interval estimation of the scale parameter, as well as hypothesis testing and prediction, is presented. In many cases, the use of substantial prior information can significantly increase the quality of the inferences and reduce the amount of testing required. Some Bayes estimators and predictors are derived in closed-forms. Highest posterior density estimators and credibility intervals can be computed using iterative methods. Bayes rules for testing one- and two-sided hypotheses are also provided. An illustrative numerical example is included.     

New concepts and applications in growth phenomena

This paper reviews earlier research on new concepts for growth curve modelling, expands on this work and explores the practical feedback from industry. The definition of growth is refined and the consequences of the new concepts are considered in greater depth. Reference is made to diffusion as a subset of growth, and the importance of the generalized Norden model is highlighted by comparison with the Bass and logistic models for dominant growth or diffusion elements. The extensive implications of the new principles are explored, including growth differentials, and open and closed temporal and spatial systems. Various examples are given, distinguishing between 'micro' (one product) and 'macro' (tens of products) case studies. These range from a study of the Spanish tractor economy to the production of a helicopter. Cashflows are shown to be a growth process. Finally, the stochastic issues are examined both for a dominant growth element and, by use of an example of the construction of a new model, at the level of stochastic differential equations.     

Non-parametric maximum likelihood estimation of interval-censored failure time data subject to misclassification

The paper considers non-parametric maximum likelihood estimation of the failure time distribution for interval-censored data subject to misclassification. Such data can arise from two types of observation scheme; either where observations continue until the first positive test result or where tests continue regardless of the test results. In the former case, the misclassification probabilities must be known, whereas in the latter case, joint estimation of the event-time distribution and misclassification probabilities is possible. The regions for which the maximum likelihood estimate can only have support are derived. Algorithms for computing the maximum likelihood estimate are investigated and it is shown that algorithms appropriate for computing non-parametric mixing distributions perform better than an iterative convex minorant algorithm in terms of time to absolute convergence. A profile likelihood approach is proposed for joint estimation. The methods are illustrated on a data set relating to the onset of cardiac allograft vasculopathy in post-heart-transplantation patients.     

Bootstrap methods for stationary functional time series

Bootstrap methods for estimating the long-run covariance of stationary functional time series are considered. We introduce a versatile bootstrap method that relies on functional principal component analysis, where principal component scores can be bootstrapped by maximum entropy. Two other bootstrap methods resample error functions, after the dependence structure being modeled linearly by a sieve method or nonlinearly by a functional kernel regression. Through a series of Monte-Carlo simulation, we evaluate and compare the finite-sample performances of these three bootstrap methods for estimating the long-run covariance in a functional time series. Using the intraday particulate matter (\(\hbox {PM}_{10}\)) dataset in Graz, the proposed bootstrap methods provide a way of constructing the distribution of estimated long-run covariance for functional time series.     

Seasonal patterns of fertility measures: theory and data

The distribution of births by month exhibits a seasonal pattern in most populations. The monthly marital fertility rate for an area of Bangladesh provides a good example of the seasonal periodicity. Seasonal patterns of measures of reproduction in a population of married women are considered. Equations are developed that predict the seasonal patterns of these alternative measures under the assumption that the fertility rate (R) follows a trigonometric curve. This is followed by an empirical analysis of the measures in a Bangladesh population that has a pronounced seasonal fertility. The investigation is intended both to validate the theoretical framework developed in the 1st part of the paper as well as to determine whether seasonal variation in actual populations is sufficiently large to affect the alternative measures significantly. 4 measures are considered: pregnancy prevalence (PP)--the proportion of married women who are pregnant at the survey date; mean open birth interval (MOI)--the time from the last live birth to the date of the survey for parous women and from the time of marriage to the date of the survey for nulliparous women; mean closed interval-birth (MCIB)--the mean interval between the last 2 live births for married women who have a birth in the period immediately preceding the survey date; and mean closed interval-woman (MCIW)--the mean interval between the last 2 live births for women who have had at least 2 children by the time of the survey. It is assumed that the seasonal pattern of the fertility rate of a population follows a cosine curve and that there is no trend in annual fertility from year to year. The lag and relative variability of the other measures are considered in comparison with the fertility rate curve. The predictions from this theoretical effort, when compared with observed patterns and trigonometric regression results for each measure in data from Bangladesh, are shown to be quite accurate. The figure and regression results show that R, PP, and MOI have definite seasonal periodicity, but MCIB and MCIW do not display any seasonal patterns. If there is a secular trend in fertility in addition to seasonality, these relationships between the seasonal patterns of the measures may no longer hold. There is a disadvantage to using closed interval measures, for they are unable to detect effects of limiting of childbearing in a population since they are based only on information from women who have births.