On the complex dynamics of functional-coefficients nonlinear autoregressive time series models

In this article, a new chaotic functional-coefficient nonlinear autoregressive time series model is formulated. Asymptotic stability of the equilibria of the skeleton is studied. Complex dynamics of the proposed model are investigated by means of the bifurcations, time series diagrams, and phase portraits. The effects of noise intensity on its dynamics and the intermittency phenomenon are also discussed via simulation. Two chaotic indicators, namely, the fractal dimension and the Lyapunov exponent methods are investigated for the model.     

A stratified unrelated question randomized response model using Neyman allocation

This paper applies stratified random sampling using Neyman allocation to Mangat et al. (1992 Mangat, N.S., Singh, R., Singh, S. (1992). An improved unrelated question randomized response strategy. Cal. Stat. Assoc. Bull. 42:277–281.[Crossref] , [Google Scholar]) unrelated question randomized response (RR) strategy for both completely truthful reporting and less than completely truthful reporting. It is shown that, for the prior information given, our new model is more efficient in terms of variance (in the case of completely truthful reporting) and mean square error (in terms of less than completely truthful reporting) than Kim and Elam's (2007 Kim, J.M., Elam, M.E. (2007). A stratified unrelated question randomized response model. Stat. Papers 48:215–233.[Crossref], [Web of Science ®] , [Google Scholar]) model. Numerical illustrations and graphs are also given in support of the present study.     

Bayesian prediction of order statistics with fixed and random sample sizes based on k-record values from Pareto distribution

In this paper, the two-parameter Pareto distribution is considered and the problem of prediction of order statistics from a future sample and that of its geometric mean are discussed. The Bayesian approach is applied to construct predictors based on observed k-record values for the cases when the future sample size is fixed and when it is random. Several Bayesian prediction intervals are derived. Finally, the results of a simulation study and a numerical example are presented for illustrating all the inferential procedures developed here.     

An optional randomized response model for estimating a rare sensitive attribute using Poisson distribution

The crux of this article is to estimate the mean of the number of persons possessing a rare sensitive attribute based on the Mangat (1991 Mangat, N.S. (1991). An optional randomized response sampling technique using non–stigmatized attribute. Statistica. 51(4):595–602. [Google Scholar]) randomization device by utilizing the Poisson distribution in simple random sampling and stratified sampling. Properties of the proposed randomized response (RR) model have been studied along with recommendations. It is also shown that the proposed model is more efficient than that of Land et al. (2011 Land, M., Singh, S., Sedory, S.A. (2011). Estimation of a rare attribute using Poisson distribution. Statistics. 46(3):351–360.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) in simple random sampling and that of Lee et al. (2013 Lee, G.S., Uhm, D., Kim, J.M. (2013). Estimation of a rare sensitive attribute in stratified sampling using Poisson distribution. Statistics. 47(3):575–589.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) in stratified random sampling when the proportion of persons possessing a rare unrelated attribute is known. Numerical illustrations are also given in support of the present study.     

Slashed generalized Rayleigh distribution

We introduce a new family of distributions suitable for fitting positive data sets with high kurtosis which is called the Slashed Generalized Rayleigh Distribution. This distribution arises as the quotient of two independent random variables, one being a generalized Rayleigh distribution in the numerator and the other a power of the uniform distribution in the denominator. We present properties and carry out estimation of the model parameters by moment and maximum likelihood (ML) methods. Finally, we conduct a small simulation study to evaluate the performance of ML estimators and analyze real data sets to illustrate the usefulness of the new model.     

Gamma-Maxwell distribution

A new two-parameter distribution, the gamma-Maxwell distribution, isproposed and studied. We generate the new distribution using the gamma-G generator of distributions. The proposal distribution can be seen as an extension of the Maxwell distribution with more flexibility in terms of the distribution asymmetry and kurtosis. We study some probability properties, discuss maximum-likelihood estimation and present a real data application indicating that the new distribution can improve the ordinary Maxwell distribution in fitting real data.     

Bayesian quantile regression for skew-normal linear mixed models

Linear mixed models have been widely used to analyze repeated measures data which arise in many studies. In most applications, it is assumed that both the random effects and the within-subjects errors are normally distributed. This can be extremely restrictive, obscuring important features of within-and among-subject variations. Here, quantile regression in the Bayesian framework for the linear mixed models is described to carry out the robust inferences. We also relax the normality assumption for the random effects by using a multivariate skew-normal distribution, which includes the normal ones as a special case and provides robust estimation in the linear mixed models. For posterior inference, we propose a Gibbs sampling algorithm based on a mixture representation of the asymmetric Laplace distribution and multivariate skew-normal distribution. The procedures are demonstrated by both simulated and real data examples.     

The performance of model averaged tail area confidence intervals

We investigate the exact coverage and expected length properties of the model averaged tail area (MATA) confidence interval proposed by Turek and Fletcher, CSDA, 2012, in the context of two nested, normal linear regression models. The simpler model is obtained by applying a single linear constraint on the regression parameter vector of the full model. For given length of response vector and nominal coverage of the MATA confidence interval, we consider all possible models of this type and all possible true parameter values, together with a wide class of design matrices and parameters of interest. Our results show that, while not ideal, MATA confidence intervals perform surprisingly well in our regression scenario, provided that we use the minimum weight within the class of weights that we consider on the simpler model.     

Parallel and interacting stochastic approximation annealing algorithms for global optimisation

We present the parallel and interacting stochastic approximation annealing (PISAA) algorithm, a stochastic simulation procedure for global optimisation, that extends and improves the stochastic approximation annealing (SAA) by using population Monte Carlo ideas. The efficiency of standard SAA algorithm crucially depends on its self-adjusting mechanism which presents stability issues in high dimensional or rugged optimisation problems. The proposed algorithm involves simulating a population of SAA chains that interact each other in a manner that significantly improves the stability of the self-adjusting mechanism and the search for the global optimum in the sampling space, as well as it inherits SAA desired convergence properties when a square-root cooling schedule is used. It can be implemented in parallel computing environments in order to mitigate the computational overhead. As a result, PISAA can address complex optimisation problems that it would be difficult for SAA to satisfactory address. We demonstrate the good performance of the proposed algorithm on challenging applications including Bayesian network learning and protein folding. Our numerical comparisons suggest that PISAA outperforms the simulated annealing, stochastic approximation annealing, and annealing evolutionary stochastic approximation Monte Carlo.