Modelling the duration of spells of withdrawal from a female dominated labour force

It is of particular importance to nurse manpower planners to determine the likelihood and time-scale of return to service of nurses who have temporarily withdrawn from the labour force. Such behaviour is common in female dominated professions and is referred to in the nursing literature as being in ‘limbo’. In this paper we develop non-parametric, parametric and network models to describe the distribution of stay in ‘limbo’ and estimate the probability of return to service. These models are fitted to data from the Northern Ireland nursing service. It is inherent in the problem that the data are statistically incomplete. This must be taken into account in the statistical methods used to estimate the various model parameters.     

Nonparametric maximum likelihood estimation for artificially truncated absence data

In manpower planning it is cornmoniy tue case tnat employees withuraw from active service for a period of time before returning to take up post at a later date. Such periods of absence are frequently of major concern to employers who are anxious to ensure that employees return as soon as possible. The distribution of duration of such periods of absence are therefore of considerable interest as is the probability that such employees will ever return to active service. In this paper we derive a nonparametric estimator for such a lifetime distribution based on renewal data which are subject to various forms of incompleteness, namely right censoring, left and right truncation, and forward recurrence. Artificial truncation is used to ensure that the data are time homogeneous. A nonparametric maximum likelihood estimator for the lifetime.     

A two phases queueing system with Bernoulli vacation schedule under multiple vacation policy

This paper deals with a single server Poisson arrival queue with two phases of heterogeneous service along with a Bernoulli schedule vacation model, where after two successive phases service the server either goes for a vacation with probability p (0≤p≤1) or may continue to serve the next unit, if any, with probability q(=1−p). Further the concept of multiple vacation policy is also introduced here. We obtained the queue size distributions at a departure epoch and at a random epoch, Laplace Stieltjes Transform of the waiting time distribution and busy period distribution along with some mean performance measures. Finally we discuss some statistical inference related issues.     

Time Series Modeling of Degradation Due to Outdoor Weathering

Environmental variables have an important effect on the reliability of many products such as coatings and polymeric composites. Long-term prediction of the performance or service life of such products must take into account the probabilistic/stochastic nature of the outdoor weather. In this article, we propose a time series modeling procedure to model the time series data of daily accumulated degradation. Daily accumulated degradation is the total amount of degradation accrued within one day and can be obtained by using a degradation rate model for the product and the weather data. The fitted model of the time series can then be used to estimate the future distribution of cumulative degradation over a period of time, and to compute reliability measures such as the probability of failure. The modeling technique and estimation method are illustrated using the degradation of a solar reflector material. We also provide a method to construct approximate confidence intervals for the probability of failure.     

Data reconciliation of nonnormal observations with nonlinear constraints

This paper presents a new method for the reconciliation of data described by arbitrary continuous probability distributions, with the focus on nonlinear constraints. The main idea, already applied to linear constraints in a previous paper, is to restrict the joint prior probability distribution of the observed variables with model constraints to get a joint posterior probability distribution. Because in general the posterior probability density function cannot be calculated analytically, it is shown that it has decisive advantages to sample from the posterior distribution by a Markov chain Monte Carlo (MCMC) method. From the resulting sample of observed and unobserved variables various characteristics of the posterior distribution can be estimated, such as the mean, the full covariance matrix, marginal posterior densities, as well as marginal moments, quantiles, and HPD intervals. The procedure is illustrated by examples from material flow analysis and chemical engineering.     

The discrete Lindley distribution: properties and applications

Modelling count data is one of the most important issues in statistical research. In this paper, a new probability mass function is introduced by discretizing the continuous failure model of the Lindley distribution. The model obtained is over-dispersed and competitive with the Poisson distribution to fit automobile claim frequency data. After revising some of its properties a compound discrete Lindley distribution is obtained in closed form. This model is suitable to be applied in the collective risk model when both number of claims and size of a single claim are implemented into the model. The new compound distribution fades away to zero much more slowly than the classical compound Poisson distribution, being therefore suitable for modelling extreme data.     

The negative contagion reflection of the polya-eggenberger distribution

The Polya-Eggenberger distribution Involves drawing a ball from an urn containing black and white balls and, after each drawing, returning the ball together with s balls of the same color, The model represents positive contagion since the added balls are the same color as the one drawn, See Johnson and Kotz, (1977),

This paper derives and examines the probability distribution which results from the Polya-Eggenberger model with only one change namely, the s additional balls added after each drawing are of the opposite color, producing a negative contagion model.

Formulas in closed form are presented for the probability distribution function, the mean and variance, all binomial moments and, where s is greater than or equal to the number of balls in the urn at start, the mode, A formula for the mode is conjectured where s is less than the number of balls in the urn at start.

Finally, the probability of obtaining k black balls in n drawings is shown in certain instances to be equal to A_nk/n!

where A_nk are the Eulerian numbers.     

An application of lomax distributions in receiver operating characteristic(roc)curve analysis

Receiver operating characteristic(ROC)curves are useful for studying the performance of diagnostic tests. ROC curves occur in many fields of applications including psychophysics, quality control and medical diagnostics. In practical situations, often the responses to a diagnostic test are classified into a number of ordered categories. Such data are referred to as ratings data. It is typically assumed that the underlying model is based on a continuous probability distribution. The ROC curve is then constructed from such data using this probability model. Properties of the ROC curve are inherited from the model. Therefore, understanding the role of different probability distributions in ROC modeling is an interesting and important area of research. In this paper the Lomax distribution is considered as a model for ratings data and the corresponding ROC curve is derived. The maximum likelihood estimation procedure for the related parameters is discussed. This procedure is then illustrated in the analysis of a neurological data example.     

Computational analysis of the queue with working breakdowns and delaying repair under a Bernoulli-schedule-controlled policy

This paper considers a single server queueing system with working breakdowns and delaying repair under a Bernoulli-schedule-controlled policy. At a breakdown instant, the system either goes to repair period immediately with probability p, or continues to provide auxiliary service for the current customers with probability q = 1 ? p. While the system resides in the auxiliary service period, it may go to repair period if there is no customer at the epoch of service completion or the occurrence of breakdown. By using the matrix analytic method and the spectral expansion method, we respectively obtain the steady state distribution to make the straightforward computation of performance measures and the Laplace-Stieltjes transform of the stationary sojourn time of an arbitrary customer. In addition, some numerical examples are presented to show the impact of parameters on the performance measures.     

Combining Inverse Probability Weighting and Multiple Imputation to Improve Robustness of Estimation

Inverse probability weighting (IPW) and multiple imputation are two widely adopted approaches dealing with missing data. The former models the selection probability, and the latter models data distribution. Consistent estimation requires correct specification of corresponding models. Although the augmented IPW method provides an extra layer of protection on consistency, it is usually not sufficient in practice as the true data‐generating process is unknown. This paper proposes a method combining the two approaches in the same spirit of calibration in sampling survey literature. Multiple models for both the selection probability and data distribution can be simultaneously accounted for, and the resulting estimator is consistent if any model is correctly specified. The proposed method is within the framework of estimating equations and is general enough to cover regression analysis with missing outcomes and/or missing covariates. Results on both theoretical and numerical investigation are provided.     

On some probability distributions associated with random walks

The probability distribution of the total number of games to ruin in a gambler's ruin random walk with initial position n, the probability distribution of the total size of an epidemic starting with n cases and the probability distribution of the number of customers served during a busy period M/M/1 when the service starts with n waiting customers are identical. All these can be easily obtained by using Lagrangian expansions instead of long combinatorial methods. The binomial, trinomial, quadrinomial and polynomial random walks of a particle have been considered with an absorbing barrier at 0 when the particle starts its walks from a point n, and the pgfs. and the probability distributions of the total number of jumps (trials) before absorption at 0 have been obtained. The values for the mean and variance of such walks have also been given.     

Estimation of the Parameters of a Two-phase Tandem Queue with a Second Optional Service

In this article, we consider a two-phase tandem queueing model with a second optional service. In this model, the service is done by two phases. The first phase of service is essential for all customers and after the completion of the first phase of service, any customer receives the second phase of service with probability α, or leaves the system with probability 1 ? α. Also, there are two heterogeneous servers which work independently, one of them providing the first phase of service and the other a second phase of service. In this model, our main purpose is to estimate the parameters of the model, traffic intensity, and mean system size, in the steady state, via maximum likelihood and Bayesian methods. Furthermore, we find asymptotic confidence intervals for mean system size. Finally, by a simulation study, we compute the confidence levels and mean length for asymptotic confidence intervals of mean system size with a nominal level 0.95.     

Analysis of a two incongruent arrivals and services M/G/1 queue with random feedback

A queuing system with two incongruent arrivals and services is considered. Two kinds of customers enter the system by Poisson process and the service times are assumed to have general distribution. After first kind service completion, it may feedback to repeat the first service, leave the system or go to give second service. The same policy is applied for the other kind of customer. All stochastic processes involved in this system are independent. We derive the probability generating function for each kind and for the system that yield the performance measures. Some numerical approaches examined the validity of the results.     

Analyzing dependence in incidence of diabetes and heart problem using generalized bivariate geometric models with covariates

For analyzing incidence data on diabetes and health problems, the bivariate geometric probability distribution is a natural choice but remained unexplored largely due to lack of models linking covariates with the probabilities of bivariate incidence of correlated outcomes. In this paper, bivariate geometric models are proposed for two correlated incidence outcomes. The extended generalized linear models are developed to take into account covariate dependence of the bivariate probabilities of correlated incidence outcomes for diabetes and heart diseases for the elderly population. The estimation and test procedures are illustrated using the Health and Retirement Study data. Two models are shown in this paper, one based on conditional-marginal approach and the other one based on the joint probability distribution with an association parameter. The joint model with association parameter appears to be a very good choice for analyzing the covariate dependence of the joint incidence of diabetes and heart diseases. Bootstrapping is performed to measure the accuracy of estimates and the results indicate very small bias.     

A tampered Brownian motion process model for partial step-stress accelerated life testing

In partial step-stress accelerated life testing, models extrapolating data obtained under more severe conditions to infer the lifetime distribution under normal use conditions are needed. Bhattacharyya (Invited paper for 46th session of the ISI, 1987) proposed a tampered Brownian motion process model and later derived the probability distribution from a decay process perspective without linear assumption. In this paper, the model is described and the features of the failure time distribution are discussed. The maximum likelihood estimates of the parameters in the model and their asymptotic properties are presented. An application of models for step-stress accelerated life test to fields other than engineering is described and illustrated by applying the tampered Brownian motion process model to data taken from a clinical trial.     

中国证券市场股价指数VaR研究

This paper compares the models in Predicting the VaR of stock price indexes in China. Our resultsindicate that the normal model usually underestimate the VaR when given probability is 0.01 or 0.02, andthe weighted normal model usually overestimate the VaR when given probability is 0.04 or 0.05. Thehistorical simulating model and Logistic distribution model are superior to normal model and to weightednomal model in predicting the VaR.     

On predictive distributions and Bayesian networks

In this paper we are interested in discrete prediction problems for a decision-theoretic setting, where the task is to compute the predictive distribution for a finite set of possible alternatives. This question is first addressed in a general Bayesian framework, where we consider a set of probability distributions defined by some parametric model class. Given a prior distribution on the model parameters and a set of sample data, one possible approach for determining a predictive distribution is to fix the parameters to the instantiation with the maximum a posteriori probability. A more accurate predictive distribution can be obtained by computing the evidence (marginal likelihood), i.e., the integral over all the individual parameter instantiations. As an alternative to these two approaches, we demonstrate how to use Rissanen's new definition of stochastic complexity for determining predictive distributions, and show how the evidence predictive distribution with Jeffrey's prior approaches the new stochastic complexity predictive distribution in the limit with increasing amount of sample data. To compare the alternative approaches in practice, each of the predictive distributions discussed is instantiated in the Bayesian network model family case. In particular, to determine Jeffrey's prior for this model family, we show how to compute the (expected) Fisher information matrix for a fixed but arbitrary Bayesian network structure. In the empirical part of the paper the predictive distributions are compared by using the simple tree-structured Naive Bayes model, which is used in the experiments for computational reasons. The experimentation with several public domain classification datasets suggest that the evidence approach produces the most accurate predictions in the log-score sense. The evidence-based methods are also quite robust in the sense that they predict surprisingly well even when only a small fraction of the full training set is used.     

ESTIMATION OF PARAMETERS IN THE MARGINAL FISHER DISTRIBUTION1

The Fisher distribution is frequently used as a model for the probability distribution of directional data, which may be specified either in terms of unit vectors or angular co-ordinates (co-latitude and azimuth). If, in practical situations, only the co-latitudes can be observed, the available data must be regarded as a sample from the corresponding marginal distribution. This paper discusses the estimation by Maximum Likelihood (ML) and the Method of Moments of the two parameters of this marginal Fisher distribution. The moment estimators are generally simpler to compute than the ML estimators, and have high asymptotic efficiency.     

Marshall–Olkin q-Weibull distribution and max–min processes

In this paper, we introduce a new probability model known as Marshall–Olkin q-Weibull distribution. Various properties of the distribution and hazard rate functions are considered. The distribution is applied to model a biostatistical data. The corresponding time series models are developed to illustrate its application in times series modeling. We also develop different types of autoregressive processes with minification structure and max–min structure which can be applied to a rich variety of contexts in real life. Sample path properties are examined and generalization to higher orders are also made. The model is applied to a time series data on daily discharge of Neyyar river in Kerala, India.     

Higher-order and non-stationary properties of lampard's stochastic reversible counter system

Summary. In this paper we discuss higher-order and non-stationary properties of LAM-PAED'S (1968) stochastic reversible counter model whose output is a stochastic point process with MABKOV dependent time intervals. Statistical properties of the primed counter system, are also considered. A characterization of the probability density function of the interval lengths and the probability distribution of the counts by means of a mixed POXSSOH process is then derived. The distribution of counts of the primed counter system is shown to be mathematically equivalent to the birth-and-death process with immigration. A generalization of the problem when initial counts and incremental counts have independent distributions is also discussed