Identification of parametric Rasch-type models

In modern Item Response Theory, the Rasch model is viewed as a Generalized Linear Mixed Model, where the item parameters correspond to the fixed-effects, whereas the person specific parameters are the random-effects. The statistical model, bearing on the observable variables only, is obtained after integrating out the random-effects. Although it is widely accepted that the parameters of this model are identified, it is hard to find a correct justification. Furthermore, the meaning of the parameters of the Rasch model – as well as of its extensions – is typically based on the fixed-effects specification of the model, that is, when the person specific parameters are also treated as fixed-effects. The contribution of this paper is to provide an explicit proof of the identification of the random-effects Rasch model. The proof is valid for a large class of Rasch-type models. It is also shown that such a proof can be applied to analyze the identification of Explanatory Rasch Models. Finally, the meaning of the parameters of interest with respect to the different data generating process is discussed.     

An Approximate Closed Form Solution for a Quadratic Dose–Response Model*

The quadratic dose–response model is often used in radiobiology studies. Since the model is nonlinear, the least squares estimators of the parameters of the model are determined by an iterative procedure. We give a simple closed form approximation for estimating the parameters of the model.     

Approximate conditional inference and the linear functional model

Approximate conditional inference is developed for the slope parameter of the linear functional model with two variables. It is shown that the model can be transformed so that the slope parameter becomes an angle and nuisance parameters are radial distances. If the nuisance parameters are known an exact confidence interval based on a location-type conditional distribution is available for the angle. More gen¬erally, confidence distributions are used to average the conditional distribution over the nuisance parameters yielding an approximate conditional confidence interval that reflects the precision indicated by the data. An example is analyzed.     

Option pricing under model and parameter uncertainty using predictive densities

The theoretical price of a financial option is given by the expectation of its discounted expiry time payoff. The computation of this expectation depends on the density of the value of the underlying instrument at expiry time. This density depends on both the parametric model assumed for the behaviour of the underlying, and the values of parameters within the model, such as volatility. However neither the model, nor the parameter values are known. Common practice when pricing options is to assume a specific model, such as geometric Brownian Motion, and to use point estimates of the model parameters, thereby precisely defining a density function.We explicitly acknowledge the uncertainty of model and parameters by constructing the predictive density of the underlying as an average of model predictive densities, weighted by each model's posterior probability. A model's predictive density is constructed by integrating its transition density function by the posterior distribution of its parameters. This is an extension to Bayesian model averaging. Sampling importance-resampling and Monte Carlo algorithms implement the computation. The advantage of this method is that rather than falsely assuming the model and parameter values are known, inherent ignorance is acknowledged and dealt with in a mathematically logical manner, which utilises all information from past and current observations to generate and update option prices. Moreover point estimates for parameters are unnecessary. We use this method to price a European Call option on a share index.     

Bayesian analysis of spatial point processes in the neighbourhood of Voronoi networks

A model for an inhomogeneous Poisson process with high intensity near the edges of a Voronoi tessellation in 2D or 3D is proposed. The model is analysed in a Bayesian setting with priors on nuclei of the Voronoi tessellation and other model parameters. An MCMC algorithm is constructed to sample from the posterior, which contains information about the unobserved Voronoi tessellation and the model parameters. A major element of the MCMC algorithm is the reconstruction of the Voronoi tessellation after a proposed local change of the tessellation. A simulation study and examples of applications from biology (animal territories) and material science (alumina grain structure) are presented.     

Optimal designs for the power logistic model

We investigate non-sequential designs for estimating model parameters in a power logistic model when the power is assumed to be approximately known and only the ranges for the other two parameters are available. The sensitivity of these designs to nominal values of all the three parameters are studied and our proposed optimal designs are shown to be reasonably robust under moderate deviation from the assumed model. An application to a toxicity experiment involving adult beetles is discussed, including the benefits of using an optimal design.     

Mixed model for analyzing geographic variability in mortality rates

"A mixed model is proposed for the analysis of geographic variability in mortality rates. In addition to demographic parameters and random geographic parameters, the model includes additional random-effects parameters to adjust for extra-Poisson variability. The model uses a gamma-Poisson distribution with a random scale parameter having an inverse gamma prior. An empirical Bayes approach is used to estimate relative risks for geographic regions and annual rates for demographic groups within each region. Lung cancer in Missouri is used to motivate and illustrate the procedure."     

Weibull accelerated life tests when there are competing causes of failure

Accelerated life testing of a product under more severe than normal conditions is commonly used to reduce test time and costs. Data collected at such accelerated conditions are used to obtain estimates of the parameters of a stress translation function. This function is then used to make inference about the product's life under normal operating conditions. We consider the problem of accelerated life tests when the product of interest is a p component series system. Each of the components is assumed to have an independent Weibull time to failure distribution with different shape parameters and different scale parameters which are increasing functions stress. A general model i s used for the scale parameter includes the standard engineering models as special This model also has an appealing biological interpretation     

A proportional hazards regression model with change-points in the baseline function

In this article, we consider a new regression model for counting processes under a proportional hazards assumption. This model is motivated by the need of understanding the evolution of the booking process of a railway company. The main novelty of the approach consists in assuming that the baseline hazard function is piecewise constant, with unknown times of jump (these times of jump are estimated from the data as model parameters). Hence, the parameters of the model can be separated into two different types: parameters that measure the influence of the covariates, and parameters from a multiple change-point model for the baseline. Cox??s semiparametric regression can be seen as a limit case of our model. We develop an iterative procedure to estimate the different parameters, and a test procedure that allows to perform change-point detection in the baseline. Our technique is supported by simulation studies and a real data analysis, which show that our model can be a reasonable alternative to Cox??s regression model, particularly in the presence of tied event times.     

Parametric estimation of location and scale parameters in ranked set sampling

This paper develops parametric inference for the parameters of location-scale family of distributions based on a ranked set sample. Likelihood function incorporates within-set ranking errors into the model through a missing data mechanism. The maximum likelihood estimators of the location-scale and missing data model parameters are constructed and an EM-algorithm is provided. It is shown that the proposed estimator is robust against imperfect ranking error and provides higher efficiency over its competitors.     

Application of Dirichlet mixture of normals in growth curve models

In this paper, we present growth curve models with an auxiliary variable which contains an uncertain data distribution based on mixtures of standard components, such as normal distributions. The multimodality of the auxiliary random variable motivates and necessitates the use of mixtures of normal distributions in our model. We have observed that Dirichlet process priors, composed of discrete and continuous components, are appropriate in addressing the two problems of determining the number of components and estimating the parameters simultaneously and are especially useful in the aforementioned multimodal scenario. A model for the application of Dirichlet mixture of normals (DMN) in growth curve models under Bayesian formulation is presented and algorithms for computing the number of components, as well as estimating the parameters are also rendered. The simulation results show that our model gives improved goodness of fit statistics over models without DMN and the estimates for the number of components and for parameters are reasonably accurate.     

Laplace-normal mixtures fitted to wind shear data

A mixture model with Laplace and normal components is fitted to wind shear data available in grouped form. A set of equations is presented for iteratively estimating the parameters of the model using an application of the EM algorithm. Twenty-four sets of data are examined with this technique, and the model is found to give a good fit to the data. Some hypotheses about the parameters in the model are discussed in light of the estimates obtained.     

Design issues for population growth models

We briefly review and discuss design issues for population growth and decline models. We then use a flexible growth and decline model as an illustrative example and apply optimal design theory to find optimal sampling times for estimating model parameters, specific parameters and interesting functions of the model parameters for the model with two real applications. Robustness properties of the optimal designs are investigated when nominal values or the model is mis-specified, and also under a different optimality criterion. To facilitate use of optimal design ideas in practice, we also introduce a website for generating a variety of optimal designs for popular models from different disciplines.     

Bayesian updating of atmospheric dispersion after a nuclear accident

Summary.  We consider a Bayesian forecasting system to predict the dispersal of contamination on a large scale grid in the event of an accidental release of radioactivity. The statistical model is built on a physical model for atmospheric dispersion and transport called MATCH. Our spatiotemporal model is a dynamic linear model where the state parameters are the (essentially, deterministic) predictions of MATCH; the distributions of these are updated sequentially in the light of monitoring data. One of the distinguishing features of the model is that the number of these parameters is very large (typically several hundreds of thousands) and we discuss practical issues arising in its implementation as a realtime model. Our procedures have been checked against a variational approach which is used widely in the atmospheric sciences. The results of the model are applied to test data from a tracer experiment.     

A recursive approach to parameter estimation in regression and time series models

In this paper we discuss the recursive (or on line) estimation in (i) regression and (ii) autoregressive integrated moving average (ARIMA) time series models. The adopted approach uses Kalman filtering techniques to calculate estimates recursively. This approach is used for the estimation of constant as well as time varying parameters. In the first section of the paper we consider the linear regression model. We discuss recursive estimation both for constant and time varying parameters. For constant parameters, Kalman filtering specializes to recursive least squares. In general, we allow the parameters to vary according to an autoregressive integrated moving average process and update the parameter estimates recursively. Since the stochastic model for the parameter changes will "be rarely known, simplifying assumptions have to be made. In particular we assume a random walk model for the time varying parameters and show how to determine whether the parameters are changing over time. This is illustrated with an example.     

Bayesian parameter estimation via variational methods

We consider a logistic regression model with a Gaussian prior distribution over the parameters. We show that an accurate variational transformation can be used to obtain a closed form approximation to the posterior distribution of the parameters thereby yielding an approximate posterior predictive model. This approach is readily extended to binary graphical model with complete observations. For graphical models with incomplete observations we utilize an additional variational transformation and again obtain a closed form approximation to the posterior. Finally, we show that the dual of the regression problem gives a latent variable density model, the variational formulation of which leads to exactly solvable EM updates.     

Optimality of quasi-score in the multivariate mean–variance model with an application to the zero-inflated Poisson model with measurement errors

In a multivariate mean–variance model, the class of linear score (LS) estimators based on an unbiased linear estimating function is introduced. A special member of this class is the (extended) quasi-score (QS) estimator. It is ‘extended’ in the sense that it comprises the parameters describing the distribution of the regressor variables. It is shown that QS is (asymptotically) most efficient within the class of LS estimators. An application is the multivariate measurement error model, where the parameters describing the regressor distribution are nuisance parameters. A special case is the zero-inflated Poisson model with measurement errors, which can be treated within this framework.     

Economic-statistical design of X ¥ charts for non-normal data by considering quality loss

When the X ¥ control chart is used to monitor a process, three parameters should be determined: the sample size, the sampling interval between successive samples, and the control limits of the chart. Duncan presented a cost model to determine the three parameters for an X ¥ chart. Alexander et al. combined Duncan's cost model with the Taguchi loss function to present a loss model for determining the three parameters. In this paper, the Burr distribution is employed to conduct the economic-statistical design of X ¥ charts for non-normal data. Alexander's loss model is used as the objective function, and the cumulative function of the Burr distribution is applied to derive the statistical constraints of the design. An example is presented to illustrate the solution procedure. From the results of the sensitivity analyses, we find that small values of the skewness coefficient have no significant effect on the optimal design; however, a larger value of skewness coefficient leads to a slightly larger sample size and sampling interval, as well as wider control limits. Meanwhile, an increase on the kurtosis coefficient results in an increase on the sample size and wider control limits.     

Analysis of a jump-diffusion option pricing model with serially correlated jump sizes

This paper extends the classical jump-diffusion option pricing model to incorporate serially correlated jump sizes which have been documented in recent empirical studies. We model the series of jump sizes by an autoregressive process and provide an analysis on the underlying stock return process. Based on this analysis, the European option price and the hedging parameters under the extended model are derived analytically. Through numerical examples, we investigate how the autocorrelation of jump sizes influences stock returns, option prices and hedging parameters, and demonstrate its effects on hedging portfolios and implied volatility smiles. A calibration example based on real market data is provided to show the advantage of incorporating the autocorrelation of jump sizes.     

多指标综合评价模型鲁棒性度量方法研究

综合评价模型的鲁棒性是指模型参数（或环境参数）改变时,评价排序结果保持相对稳定的能力。通过构建综合评价模型鲁棒性度量指标,从无量纲化方法和样本容量两个方面分析度量综合评价模型的鲁棒性。实证分析结论表明,当模型参数（或环境参数）改变时,各种多指标综合评价模型存在鲁棒性差异。该研究为综合评价模型的优化选择提供了一种新的途径。