Estimation and diagnostic for skew-normal partially linear models

Partially linear models (PLMs) are an important tool in modelling economic and biometric data and are considered as a flexible generalization of the linear model by including a nonparametric component of some covariate into the linear predictor. Usually, the error component is assumed to follow a normal distribution. However, the theory and application (through simulation or experimentation) often generate a great amount of data sets that are skewed. The objective of this paper is to extend the PLMs allowing the errors to follow a skew-normal distribution [A. Azzalini, A class of distributions which includes the normal ones, Scand. J. Statist. 12 (1985), pp. 171–178], increasing the flexibility of the model. In particular, we develop the expectation-maximization (EM) algorithm for linear regression models and diagnostic analysis via local influence as well as generalized leverage, following [H. Zhu and S. Lee, Local influence for incomplete-data models, J. R. Stat. Soc. Ser. B 63 (2001), pp. 111–126]. A simulation study is also conducted to evaluate the efficiency of the EM algorithm. Finally, a suitable transformation is applied in a data set on ragweed pollen concentration in order to fit PLMs under asymmetric distributions. An illustrative comparison is performed between normal and skew-normal errors.     

Comparison of dichotomized and distributional approaches in rare event clinical trial design: a fixed Bayesian design

This research was motivated by our goal to design an efficient clinical trial to compare two doses of docosahexaenoic acid supplementation for reducing the rate of earliest preterm births (ePTB) and/or preterm births (PTB). Dichotomizing continuous gestational age (GA) data using a classic binomial distribution will result in a loss of information and reduced power. A distributional approach is an improved strategy to retain statistical power from the continuous distribution. However, appropriate distributions that fit the data properly, particularly in the tails, must be chosen, especially when the data are skewed. A recent study proposed a skew-normal method. We propose a three-component normal mixture model and introduce separate treatment effects at different components of GA. We evaluate operating characteristics of mixture model, beta-binomial model, and skew-normal model through simulation. We also apply these three methods to data from two completed clinical trials from the USA and Australia. Finite mixture models are shown to have favorable properties in PTB analysis but minimal benefit for ePTB analysis. Normal models on log-transformed data have the largest bias. Therefore we recommend finite mixture model for PTB study. Either finite mixture model or beta-binomial model is acceptable for ePTB study.     

Tests for Structural Changes in Time Series of Counts

We propose methods for detecting structural changes in time series with discrete‐valued observations. The detector statistics come in familiar L2‐type formulations incorporating the empirical probability generating function. Special emphasis is given to the popular models of integer autoregression and Poisson autoregression. For both models, we study mainly structural changes due to a change in distribution, but we also comment for the classical problem of parameter change. The asymptotic properties of the proposed test statistics are studied under the null hypothesis as well as under alternatives. A Monte Carlo power study on bootstrap versions of the new methods is also included along with a real data example.     

Model selection criteria in beta regression with varying dispersion

We address the issue of model selection in beta regressions with varying dispersion. The model consists of two submodels, namely: for the mean and for the dispersion. Our focus is on the selection of the covariates for each submodel. Our Monte Carlo evidence reveals that the joint selection of covariates for the two submodels is not accurate in finite samples. We introduce two new model selection criteria that explicitly account for varying dispersion and propose a fast two step model selection scheme which is considerably more accurate and is computationally less costly than usual joint model selection. Monte Carlo evidence is presented and discussed. We also present the results of an empirical application.     

A comparison of L-, LQ-, TL-moment and maximum likelihood high quantile estimates of the GPD and GEV distribution

In this article, L-moments, LQ-moments and TL-moments of the generalized Pareto and generalized extreme-value distributions are derived up to the fourth order. The first three L-, LQ- and TL-moments are used to obtain estimators of their parameters. Performing a simulation study, high-quantile estimates based on L-, LQ-, and TL-moments are compared to the maximum likelihood estimate with respect to their sample mean squared error. This consists of identifying an optimal combination of parameters α and p both considered in the range [0, 0.5] for estimating quantiles by LQ-moments. The results show L-moment and maximum likelihood methods outperform other methods.     

Copulas related to piecewise monotone functions of the interval and associated processes

In this work, we derive the copulas related to vectors obtained from the so-called chaotic stochastic processes. These are defined by the iteration of certain piecewise monotone functions of the interval [0, 1] to some initial random variable. We study some of its properties and present some examples. Since often these types of copulas do not have closed formulas, we provide a general approximation method which converges uniformly to the true copula. Our results cover a wide class of processes, including the so-called Manneville–Pomeau processes. The general theory is applied to the parametric estimation in certain chaotic processes. A Monte Carlo simulation study is also presented.     

On fiducial generators

Fiducial inference has been gaining presence recently and it is the intention of the present article to look at the notion of fiducial generators; meaning procedures to simulate parameter values that in some sense correspond to simulations from some implicit fiducial distribution. It is well known that when the distribution has group structure, stemming from the natural pivotal associated, a fiducial may be obtained. It is in the non group distributions that there appears to be still room for finding a fiducial distribution. Recently some general procedures have been proposed for dealing with generalized fiducials, but these depend on certain choices for a structural equation or a fiducial equation, as in Hannig (2009 Hannig, J. (2009). On generalized fiducial inference. Stat. Sin. 19:491–544.[Web of Science ®] , [Google Scholar]) or Taraldsen and Lindqvist (2013 Taraldsen, G., Lindqvist, B.H. (2013). Fiducial theory and optimal inference. Ann. Stat. 41(1):323–341.[Crossref], [Web of Science ®] , [Google Scholar]), respectively. A brief presentation is made of an earlier approach to fiducial inference for multivariate parameters, as in Brillinger (1962 Brillinger, D.R. (1962). Examples bearing on the definition of fiducial probability with a bibliography. Ann. Math. Stat. 33(4):1349–1355.[Crossref] , [Google Scholar]), and the implied fiducial generator introduced in Engen and Lillegård (1997 Engen, S., Lillegård, M. (1997). Stochastic simulation conditioned on sufficient statistics. Biometrika 84(1):235–240.[Crossref], [Web of Science ®] , [Google Scholar]), trying to connect them. Three interesting non group distributions are seen; two of them, the truncated exponential and the two-parameter gamma, already reported in literature. A third non group distribution is analyzed; the inverse Gaussian, connecting the fiducial that results following Brillinger (1962 Brillinger, D.R. (1962). Examples bearing on the definition of fiducial probability with a bibliography. Ann. Math. Stat. 33(4):1349–1355.[Crossref] , [Google Scholar]), with a result pertaining confidence limits for the shape parameter in Hsieh (1990 Hsieh, H.K. (1990). Inferences on the coefficient of variation of an inverse-Gaussian distribution. Commun. Stat. - Theory Methods 19(5):1589–1605.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). In the three cases, comparisons are made with the Bayesian posteriors that have been known to be close numerically. Some discussion is made on the issue of singularities of the fiducial density and its connection with densities that do not integrate to unity. As to the case of discrete observables, some comments are made for the Bernoulli distribution, only.     

Generalized Kerridge’s inaccuracy measure for conditionally specified models

Recently, conditional Renyi’s divergence of order α and Kerridge’s inaccuracy measures are studied by Navarro et al. (2014 Navarro, J., Sunoj, S.M., Linu, M.N. (2014). Characterizations of bivariate models using some dynamic conditional information divergence measures. Commun. Stat. Theory Methods 43:1939–1948.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). In the present article, a generalized dynamic conditional Kerridge’s inaccuracy measure is introduced, which can be represented as the sum of conditional Renyi’s divergence and Renyi’s entropy. Some useful bounds are obtained using the concept of likelihood ratio order. The results are extended to weighted distributions. Sufficient conditions are obtained for the monotonicity of the proposed measure. Characterizations for bivariate exponential conditional distribution are presented based on the proposed measure.     

The exponentiated reduced Kies distribution: Properties and applications

Here we consider an exponentiated version of the reduced Kies distribution and discuss some of its properties. The parameters of the distribution are estimated by the method of maximum likelihood and illustrated with the help of certain real-life data sets. Asymptotic behavior of the maximum likelihood estimators of the parameters of the distribution is also studied by using certain simulated data sets.