Markov Beta and Gamma Processes for Modelling Hazard Rates

This paper generalizes the discrete time independent increment beta process of Hjort (1990 ), for modelling discrete failure times, and also generalizes the independent gamma process for modelling piecewise constant hazard rates ( Walker and Mallick, 1997 ). The generalizations are from independent increment to Markov increment prior processes allowing the modelling of smoothness. We derive posterior distributions and undertake a full Bayesian analysis.     

Order statistics from non-identical power function random variables

In this paper, we derive some recurrence relations satisfied by the single and the product moments of order statistics arising from n independent and non-identically distributed power function random variables. These recurrence relations will enable one to compute all the single and the product moments of all order statistics in a simple recursive manner. The results for the multiple-outlier model are deduced as special cases. The results are further generalized to the case of truncated power function random variables.     

A flexible parametric survival model which allows a bathtub-shaped hazard rate function

A new parametric (three-parameter) survival distribution, the lognormal–power function distribution, with flexible behaviour is introduced. Its hazard rate function can be either unimodal, monotonically decreasing or can exhibit a bathtub shape. Special cases include the lognormal distribution and the power function distribution, with finite support. Regions of parameter space where the various forms of the hazard-rate function prevail are established analytically. The distribution lends itself readily to accelerated life regression modelling. Applications to five data sets taken from the literature are given. Also it is shown how the distribution can behave like a Weibull distribution (with negative aging) for certain parameter values.     

Statistical power to detect violation of the proportional hazards assumption when using the Cox regression model

The use of the Cox proportional hazards regression model is wide-spread. A key assumption of the model is that of proportional hazards. Analysts frequently test the validity of this assumption using statistical significance testing. However, the statistical power of such assessments is frequently unknown. We used Monte Carlo simulations to estimate the statistical power of two different methods for detecting violations of this assumption. When the covariate was binary, we found that a model-based method had greater power than a method based on cumulative sums of martingale residuals. Furthermore, the parametric nature of the distribution of event times had an impact on power when the covariate was binary. Statistical power to detect a strong violation of the proportional hazards assumption was low to moderate even when the number of observed events was high. In many data sets, power to detect a violation of this assumption is likely to be low to modest.     

Estimating the intensity function of a power law process at the current time: time truncated case

The power law process, a nonhomogeneous Poisson process with intensity function µ(t) = (β/θ)(t/θ) , is frequently used to model the occurence of events in time. Often, an important quantity is the value of the intensity function at the current time, that is, the time when data collection is ceased. In this article, the problem of estimating this quantity is addressed when the data are time truncated, that is, when data collection is stopped at a predetermined time T. The class of multiples of the conditional MLE is suggested, and some members are analyzed. In addition, the class of estimators formed by first performing a preliminary test of significance on the parameter β is analyzed. Expressions for the bias and MSE of these estimators are derived and evaluated for several values of the parameters     

The limiting power of the durbin-watson test

We compute the limiting povwer of the Durbin-Watson test in a general linear regression model Our treatment includes previous results due to W. Kramer and H, Zeisel as well as new results. In particular, we provide new insight under which conditions the limiting power is zero or one. Stochastic-simulations correspond to our investigations.     

A generalization for two-sided power distributions and adjusted method of moments

A new rich class of generalized two-sided power (TSP) distributions, where their density functions are expressed in terms of the Gauss hypergeometric functions, is introduced and studied. In this class, the symmetric distributions are supported by finite intervals and have normal shape densities. Our study on TSP distributions also leads us to a new class of discrete distributions on {0, 1, …, k}. In addition, a new numerical method for parameter estimation using moments is given.     

Discretely Observed Diffusions: Approximation of the Continuous-time Score Function

We discuss parameter estimation for discretely observed, ergodic diffusion processes where the diffusion coefficient does not depend on the parameter. We propose using an approximation of the continuous-time score function as an estimating function. The estimating function can be expressed in simple terms through the drift and the diffusion coefficient and is thus easy to calculate. Simulation studies show that the method performs well.     

The perils with the misuse of predictive power

In early drug development, especially when studying new mechanisms of action or in new disease areas, little is known about the targeted or anticipated treatment effect or variability estimates. Adaptive designs that allow for early stopping but also use interim data to adapt the sample size have been proposed as a practical way of dealing with these uncertainties. Predictive power and conditional power are two commonly mentioned techniques that allow predictions of what will happen at the end of the trial based on the interim data. Decisions about stopping or continuing the trial can then be based on these predictions. However, unless the user of these statistics has a deep understanding of their characteristics important pitfalls may be encountered, especially with the use of predictive power. The aim of this paper is to highlight these potential pitfalls. It is critical that statisticians understand the fundamental differences between predictive power and conditional power as they can have dramatic effects on decision making at the interim stage, especially if used to re-evaluate the sample size. The use of predictive power can lead to much larger sample sizes than either conditional power or standard sample size calculations. One crucial difference is that predictive power takes account of all uncertainty, parts of which are ignored by standard sample size calculations and conditional power. By comparing the characteristics of each of these statistics we highlight important characteristics of predictive power that experimenters need to be aware of when using this approach.     

A Bayesian nonparametric estimator of a multivariate survival function

Using reinforced processes related to beta-Stacy process and generalized Pólya urn scheme jointly with a structure assumption about dependence, a Bayesian nonparametric prior and a predictive estimator for a multivariate survival function are provided. This estimator can be computed through an easy implementation of a Gibbs sampler algorithm. Moreover consistency of the estimator is studied.     

The inverse power Lindley distribution in the presence of left-censored data

In this study, classical and Bayesian inference methods are introduced to analyze lifetime data sets in the presence of left censoring considering two generalizations of the Lindley distribution: a first generalization proposed by Ghitany et al. [Power Lindley distribution and associated inference, Comput. Statist. Data Anal. 64 (2013), pp. 20–33], denoted as a power Lindley distribution and a second generalization proposed by Sharma et al. [The inverse Lindley distribution: A stress–strength reliability model with application to head and neck cancer data, J. Ind. Prod. Eng. 32 (2015), pp. 162–173], denoted as an inverse Lindley distribution. In our approach, we have used a distribution obtained from these two generalizations denoted as an inverse power Lindley distribution. A numerical illustration is presented considering a dataset of thyroglobulin levels present in a group of individuals with differentiated cancer of thyroid.     

On the robustness of the power function of the one-sample test for the negative exponential distribution

The robustness of the power function of the standard one-sample parametric test for the mean of the negative exponential distribution is examined. The main form of departure from the exponential assumption is a mixture of negative exponential components although an alternative Gamma distribution is also examined. It is found that the test is sensitive to these departures although the effect of mixtures with short tails is less dramatic than those with long tails.     

The exact power function of an exact test of a regression model against multiple separate alternatives

A common procedure for testing a regression model against separate alternatives is to check the statistical significance of predictions from the latter appended as artificial regressors in the model under test. This paper derives the previously-unknown exact small-sample power function of such a procedure. It is then demonstrated that the procedure is biased in small samples.     

Bayesian inference and prediction analysis of the power law process based on a natural conjugate prior

Under a natural conjugate prior with four hyperparameters, the importance sampling (IS) technique is applied to the Bayesian analysis of the power law process (PLP). Samples of the parameters of the PLP are obtained from IS. Based on these samples, not only the posterior analysis of parameters and some parameter functions in the PLP are performed conveniently, but also single-sample and two-sample prediction procedures are constructed easily. Furthermore, the sensitivity of the posterior mean of the parameter functions in the PLP is studied with respect to the hyperparameters of the natural conjugate prior and it can guide the selections of the hyperparameters directly. Coupled this sensitivity with the relations between the prior moments and the hyperparameters in the natural conjugate prior, it is possible to give directions about the selections of the prior moments to a certain degree. After some numerical experiments illustrate the rationality and feasibility of the proposed methods, an engineering example demonstrates its application.     

Modelling Aggregation on the Large Scale and Regularity on the Small Scale in Spatial Point Pattern Datasets

We consider a dependent thinning of a regular point process with the aim of obtaining aggregation on the large scale and regularity on the small scale in the resulting target point process of retained points. Various parametric models for the underlying processes are suggested and the properties of the target point process are studied. Simulation and inference procedures are discussed when a realization of the target point process is observed, depending on whether the thinned points are observed or not. The paper extends previous work by Dietrich Stoyan on interrupted point processes.     

Frailty models power variance function with cure fraction and latent risk factors negative binomial

In this article, we propose a flexible cure rate model, which is an extension of Cancho et al. (2011 Cancho, V.G., Rodrigues, J., de Castro, M. (2011). A flexible model for survival data with a cure rate: A Bayesian approach. J. Appl. Stat. 38:57–70.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) model, by incorporating a power variance function (PVF) frailty term in latent risk. The model is more flexible in terms of dispersion and it also quantifies the unobservable heterogeneity. The parameter estimation is reached by maximum likelihood estimation procedure and Monte Carlo simulation studies are considered to evaluate the proposed model performance. The practical relevance of the model is illustrated in a real data set of preventing cancer recurrence.     

Empirical likelihood inference for the bivariate survival function under univariate censoring

In this article, we apply the empirical likelihood method to make inference on the bivariate survival function of paired failure times by estimating the survival function of censored time with the Kaplan–Meier estimator. Adjusted empirical likelihood (AEL) confidence intervals for the bivariate survival function are developed. We conduct a simulation study to compare the proposed AEL method with other methods. The simulation study shows the proposed AEL method has better performance than other existing methods. We illustrate the proposed method by analyzing the skin graft data.     

Bayesian Semiparametric Cure Rate Model with an Unknown Threshold

Abstract.  We propose a Bayesian semiparametric model for survival data with a cure fraction. We explicitly consider a finite cure time in the model, which allows us to separate the cured and the uncured populations. We take a mixture prior of a Markov gamma process and a point mass at zero to model the baseline hazard rate function of the entire population. We focus on estimating the cure threshold after which subjects are considered cured. We can incorporate covariates through a structure similar to the proportional hazards model and allow the cure threshold also to depend on the covariates. For illustration, we undertake simulation studies and a full Bayesian analysis of a bone marrow transplant data set.     

A characterization of the inverse autocorrelation function

The inverse autocorrelation function of a weakly stationary stochastic process X_t at lag h, γⁱ _h, is shown to equal the negative of the partial correlation between random variables X_t and X_t+h after elimination of the influence of random variables X_k, k≠t5,t+h.     

A derivation of the asymptotic distribution of the partial autocorrelation function of an autoregressive process

The asymptotic distribution of the sample partial autocorrelation function of an autoregressive process is derived from that of the sample autocorrelation function.