A Multivariate Model for Repeated Failure Time Measurements

A parametric multivariate failure time distribution is derived from a frailty-type model with a particular frailty distribution. It covers as special cases certain distributions which have been used for multivariate survival data in recent years. Some properties of the distribution are derived: its marginal and conditional distributions lie within the parametric family, and association between the component variates can be positive or, to a limited extent, negative. The simple closed form of the survivor function is useful for right-censored data, as occur commonly in survival analysis, and for calculating uniform residuals. Also featured is the distribution of ratios of paired failure times. The model is applied to data from the literature     

Curved Chaotic Map Time Series Models and Their Stochastic Reversals

This paper considers two types of chaotic map time series models, including the well-known tent, logistic and binary-shift maps as special cases; these are called curved tent and curved binary families. Deterministic behaviour is investigated by invariant distributions, Lyapunov exponents, and by serial dependency. Stochastic time reversal of the families is shown to produce models which have a broader range of stochastic and chaotic properties than their deterministic counterparts. The marginal distributions may have concentrations and restricted supports and are shown to be a non-standard class of invariant distribution. Dependenc y is generally weaker with the reversed stochastic models. The work gives a broad statistical account of deterministic and stochastically reversed map models, such as are emerging in random number generation, communica tion systems and cryptography     

De minimis and equity in risk

Indices and orderings are developed for evaluating alternative strategies in the management of risk. They reflect the goals of reducing individual and collective risks, of increasing equity, and of assigning priority to the reduction and to the equity of high risks. Individual risk is defined as the (random or non-random) level of exposure to a danger. In particular the role of a lower negligibility level is investigated. A class of indices is proposed which involves two parameters, a negligibility level and a parameter of inequality aversion, and several interpretations of the indices are discussed. We provide a set of eight axioms which are necessary and sufficient for this class of indices, and we present an approach to deal with partial information on the parameters.     

Inference for quantile measures of kurtosis,peakedness, and tail weight

Insight into measures of peakedness, heavy-tailedness, and kurtosis can be gained by studying Ruppert’s ratios of interquantile ranges. They are not only monotone in Horn’s measure of peakedness when applied to the central portion of the population, but also monotone in the practical tail-index of Morgenthaler and Tukey, when applied to the tails. Non-parametric confidence intervals are found for Ruppert’s ratios, and sample sizes required to obtain such intervals for a pre-specified relative width and level are provided. In addition, the empirical power of distribution-free tests for peakedness and bimodality are found for some symmetric distributions.     

Multivariate log-skewed distributions with normal kernel and their applications

We introduce two classes of multivariate log-skewed distributions with normal kernel: the log canonical fundamental skew-normal (log-CFUSN) and the log unified skew-normal. We also discuss some properties of the log-CFUSN family of distributions. These new classes of log-skewed distributions include the log-normal and multivariate log-skew normal families as particular cases. We discuss some issues related to Bayesian inference in the log-CFUSN family of distributions, mainly we focus on how to model the prior uncertainty about the skewing parameter. Based on the stochastic representation of the log-CFUSN family, we propose a data augmentation strategy for sampling from the posterior distributions. This proposed family is used to analyse the US national monthly precipitation data. We conclude that a high-dimensional skewing function lead to a better model fit.     

Some Results Based on a Version of the Generalized Dynamic Entropies

We propose new dynamic measures of uncertainty based on the notion of generalized dynamic entropy introduced in Di Crescenzo and Longobardi (2006). These can uniquely determine distribution functions in continuous and discrete cases, and the characterizations of some well-known distributions are provided. We also define some orderings and aging notions based on the generalized dynamic measures, and prove some of their properties, obtaining as corollaries results that have recently appeared in the literature.     

Percentage points of the statistics for testing hypotheses on mean vectors of multivariate normal distributions with missing observations

A problem of testing of hypotheses on the mean vector of a multivariate normal distribution with unknown and positive definite covariance matrix is considered when a sample with a special, though not unusual, pattern of missing observations from that population is available. The approximate percentage points of the test statistic are obtained and their accuracy has been checked by comparing them with some exact percentage points which are calculated for complete samples and some special incomplete samples. The approximate percentage points are in good agreement with exact percentage points. The above work is extended to the problem of testing the hypothesis of equality of two mean vectors of two multivariate normal distributions with the same, unknown covariance matrix     

Incomplete categorical data analysis: a bayesian perspective

In this paper the Bayesian analysis of incomplete categorical data under informative general censoring proposed by Paulino and Pereira (1995) is revisited. That analysis is based on Dirichlet priors and can be applied to any missing data pattern. However, the known properties of the posterior distributions are scarce and therefore severe limitations to the posterior computations remain. Here is shown how a Monte Carlo simulation approach based on an alternative parameterisation can be used to overcome the former computational difficulties. The proposed simulation approach makes available the approximate estimation of general parametric functions and can be implemented in a very straightforward way.     

Implementing particle swarm optimization algorithm to estimate the mixture of two Weibull parameters with censored data

We present the maximum likelihood estimation (MLE) via particle swarm optimization (PSO) algorithm to estimate the mixture of two Weibull parameters with complete and multiple censored data. A simulation study is conducted to assess the performance of the MLE via PSO algorithm, quasi-Newton method and expectation-maximization (EM) algorithm for different parameter settings and sample sizes in both uncensored and censored cases. The simulation results showed that the PSO algorithm outperforms the quasi-Newton method and the EM algorithm in most cases regarding bias and root mean square errors. Two numerical examples are used to demonstrate the performance of our proposed method.