Correlation between the numbers of two types of children in a family with the mpsd for the family size

The family size (sibship size) N is regarded as an integer-valued random variable having a Modified Power Series distribution (MPSD) of Gupta (1974). The family produces two types of children, with probabilities p and q (p+q =1) . It is proved that the correlation between the numbers B and C of these children is positive or negative according as the function log f(θ) is convex or concave with respect to the function g(θ), (see Section 2). This condition is a simple and a natural extension of the one given by Rao et al (1973). Several examples are discussed to illustrate the result.     

On an application of the diffusion approximation of the generalized birth and death process for stochastic population projection

Assuming that both birth and death rates are density and time dependent, a diffusion approximation of the generalized birth and death process has been considered in this paper to obtain a suitable stochastic population model describing the population size and its moments. A simple method of estimating the parameters of the model Is discussed. The predictions of the expected size of the population, and the variance are made and compared with the corresponding census figures as well as with another deterministic projection series made for the corresponding period.     

A Discrete Distribution Arising as a Solution of a Linear Difference Equation: Extension of the Non Central Negative Binomial Distribution

This article considers a discrete distribution that arises as the dominant solution of a linear difference equation. Basic properties and various chance mechanisms that lead to this distribution are given. In particular, its formulation as a weighted distribution and a mixed Poisson process are proposed. Parameter estimation by (a) using a combination of observed frequencies and moments and (b) maximum likelihood are examined. An example of goodness of fit is considered.     

Nonparametric estimation of the sojourn time distributions for a multipath model

Summary. We use a multipath (multistate) model to describe data with multiple end points. Statistical inference based on the intermediate end point is challenging because of the problems of nonidentifiability and dependent censoring. We study nonparametric estimation for the path probability and the sojourn time distributions between the states. The methodology proposed can be applied to analyse cure models which account for the competing risk of death. Asymptotic properties of the estimators proposed are derived. Simulation shows that the methods proposed have good finite sample performance. The methodology is applied to two data sets.     

Asymptotic properties of the maximum-likelihood estimator for a class of birth-and-death processes admitting a unique stationary distribution

The asymptotic properties of the maximum-likelihood estimator of the parameter vector for a class of birth-and-death processes admitting a unique stationary distribution are studied. Also, it is shown that identifiability of the parameter vector with respect to the likelihood implies that the Fisher information matrix is of full rank. Two special cases of biological interest are presented. One of these, the exponential birth-and-death process, is proposed as a more appropriate model of density dependence than the logistic process.     

Peculiar bias properties of the ols estimator when applied to a dynamic model with autocorrelated disturbances

This study reveals that contrary to the conventional wisdom among econometricians, the bias of the OLS estimator can be quite small when the estimator is applied to a geometrically distributed lag model, y_t<ce:glyph name="dbnd6"/> α + βx _t+ λy _t-₁. + u_t, with autocorrelated disturbances, be they AR(1), MA(1), MA(2), AR(2), and ARMA(1,1). This happens when λ is large and x_tis smoothly trended (e.g., a real GNP series). In fact, the bias of the OLS estimator becomes zero at one parameter combination, and the OLS estimator performs well over a wide range around this parameter combination. By decomposing the disturbance term into two parts, the paper also explains why OLS shows such an unexpected property. These findings have both pedagogical and practical significance.     

Estimation of the mean quality-adjusted survival using a multistate model for the sojourn times

In clinical trials, it may be of interest taking into account physical and emotional well-being in addition to survival when comparing treatments. Quality-adjusted survival time has the advantage of incorporating information about both survival time and quality-of-life. In this paper, we discuss the estimation of the expected value of the quality-adjusted survival, based on multistate models for the sojourn times in health states. Semiparametric and parametric (with exponential distribution) approaches are considered. A simulation study is presented to evaluate the performance of the proposed estimator and the jackknife resampling method is used to compute bias and variance of the estimator.     

A competing risks derivation of a mixture model for the analysis of survival data

Boag (1949) and Berkson and Gage (1952) proposed a mixture model for the analysis of survival time data when aproportion of treated patients are cured. This paper presents a derivation of the Boag/Berkson-Gage mixture model as well as a eneralization of the model based on the theory of competing risks. The assumptions underlying the model are stated and discussed and a general likelihood function is obtained. Use of the model is illustrated ith data from the Stanford Heart Transplant Program.     

Construction of a semi-Markov model for the performance of a football team in the presence of missing data

Using play-by-play data from the very beginning of the professional football league in Turkey, a semi-Markov model is presented for describing the performance of football teams. The official match results of the selected teams during 55 football seasons are used and winning, drawing and losing are considered as Markov states. The semi-Markov model is constructed with transition rates inferred from the official match results. The duration between the last match of a season and the very first match of the following season is much longer than any other duration during the season. Therefore these values are considered as missing values and estimated by using expectation–maximization algorithm. The effect of the sojourn time in a state to the performance of a team is discussed as well as mean sojourn times after losing/winning are estimated. The limiting probabilities of winning, drawing and losing are calculated. Some insights about the performance of the selected teams are presented.     

A formal test for the stationarity of the incidence rate using data from a prevalent cohort study with follow-up

In a prevalent cohort study with follow-up, the incidence process is not directly observed since only the onset times of prevalent cases can be ascertained. Assessing the “stationarity” of the underlying incidence process can be important for at least three reasons, including an improvement in efficiency when estimating the survivor function. We propose, for the first time, a formal test for stationarity using data from a prevalent cohort study with follow-up. The test makes use of a characterization of stationarity, an extension of this characterization developed in this paper, and of a test for matched pairs of right censored data. We report the results from a power study assuming varying degrees of departure from the null hypothesis of stationarity. The test is also applied to data obtained as part of the Canadian Study of Health and Aging (CSHA) to verify whether the incidence rate of dementia amongst the elderly in Canada has remained constant.     

Statistical inference for the effect of magnetic brain stimulation on a motoneurone

A stochastic model is developed for the possible excitatory and inhibitory effects of a stimulus to the brain on the activity of single human motoneurones. The model consists of a Wiener process for the build-up of underlying potential, a deterministic effect due to the stimulus and a random lag from brain to muscle. Direct likelihood inference for its parameters seems impossible, and we study the use of simulation to estimate the log-likelihood for the parameters of substantive interest. Monte Carlo methods yield point and confidence interval estimates of the membrane excitability underlying excitatory and inhibitory effects. The main qualitative conclusion is that both excitatory and inhibitory effects are unambiguously present. The contribution of statistical analysis to this problem is to provide accurate and apparently reliable inference for the quantities of neurophysiological interest. More generally, our methodology has the potential to make accurate likelihood-based inferences in challenging problems, but the computational burden can be large, particularly if the model is not fully adequate for the data, as in our application.     

The asymptotic covariance structure of log-linear model estimated parameters for the multiple recapture census

This paper presents a matrix formulation for log-linear model analysis of the incomplete contingency table which arises from multiple recapture census data. Explicit matrix product expressions are given for the asymptotic covariance structure of the maximum likelihood estimators of both the log-linear model parameter vector and the predicted value vector for the observed and missing cells. These results are illustrated for data pertaining to a population of children possessing a common congenital anomaly.     

Testing for a Change of the Innovation Distribution in Nonparametric Autoregression: The Sequential Empirical Process Approach

We consider a nonparametric autoregression model under conditional heteroscedasticity with the aim to test whether the innovation distribution changes in time. To this end, we develop an asymptotic expansion for the sequential empirical process of nonparametrically estimated innovations (residuals). We suggest a Kolmogorov–Smirnov statistic based on the difference of the estimated innovation distributions built from the first ?ns?and the last n ? ?ns? residuals, respectively (0 ≤ s ≤ 1). Weak convergence of the underlying stochastic process to a Gaussian process is proved under the null hypothesis of no change point. The result implies that the test is asymptotically distribution‐free. Consistency against fixed alternatives is shown. The small sample performance of the proposed test is investigated in a simulation study and the test is applied to a data example.     

Accounting for informative non‐compliance with a bivariate exponential model in the design of endpoint trials

Failure to adjust for informative non‐compliance, a common phenomenon in endpoint trials, can lead to a considerably underpowered study. However, standard methods for sample size calculation assume that non‐compliance is non‐informative. One existing method to account for informative non‐compliance, based on a two‐subpopulation model, is limited with respect to the degree of association between the risk of non‐compliance and the risk of a study endpoint that can be modelled, and with respect to the maximum allowable rates of non‐compliance and endpoints. In this paper, we introduce a new method that largely overcomes these limitations. This method is based on a model in which time to non‐compliance and time to endpoint are assumed to follow a bivariate exponential distribution. Parameters of the distribution are obtained by equating them with the study design parameters. The impact of informative non‐compliance is investigated across a wide range of conditions, and the method is illustrated by recalculating the sample size of a published clinical trial. Copyright © 2005 John Wiley & Sons, Ltd.     

The exact hypothesis test for the shape parameter of a new two-parameter distribution with the bathtub shape or increasing failure rate function under progressive censoring with random removals

This study considers the exact hypothesis test for the shape parameter of a new two-parameter distribution with the shape of a bathtub or increasing failure rate function under type II progressive censoring with random removals, where the number of units removed at each failure time follows a binomial or a uniform distribution. Several test statistics are proposed and one numerical example is provided to illustrate the proposed hypothesis test for the shape parameter. Finally, a simulation study is performed to compare the power performances of all proposed test statistics. We concluded that the test statistic w ₁ is more attractive than other methods as it has better performance than other test statistics for most cases based on the criteria of maximum power.