Search

Scale mixtures of normal distributions are useful for statistical procedures involving symmetric and heavy-tailed data. Ferreira, Lachos, and Bolfarine (2016) defined a multivariate skewed version of these distributions that offers much-needed flexibility by combining both skewness and heavy tails. In this work, we develop a linear mixed model based on skew scale mixtures of normal distributions, with emphasis on the skew Student-t normal, skew–slash and skew–contaminated normal distributions. Using the hierarchical structure of the model, we develop maximum likelihood estimation of the model parameters via the expectation-maximization (EM) algorithm. In addition, the standard errors are obtained via the approximate information matrix and the local influence analysis is explored under some perturbation schemes. To examine the performance and the usefulness of the proposed method, we present simulation studies and analyze a real dataset. [ABSTRACT FROM AUTHOR]

In this article, we introduce a new distribution for modeling positive data sets with high kurtosis, the modified slashed generalized exponential distribution. The new model can be seen as a modified version of the slashed generalized exponential distribution. It arises as a quotient of two independent random variables, one being a generalized exponential distribution in the numerator and a power of the exponential distribution in the denominator. We studied various structural properties (such as the stochastic representation, density function, hazard rate function and moments) and discuss moment and maximum likelihood estimating approaches. Two real data sets are considered in which the utility of the new model in the analysis with high kurtosis is illustrated. [ABSTRACT FROM AUTHOR]

In this paper, we consider a generalization of the modified slash distribution. We define the new family through the quotient between an elliptically distributed random variable and the power of an exponential random variable with parameter equals to 2, both independent. We use the same idea to extend the model for the multivariate case and study general important properties from the resultant family. We perform inference by the method of moments and maximum likelihood. We present a simulation study which indicates satisfactory parameter recovery by using the estimation approaches. Illustrations reveals that it has potential for doing well in real problems. [ABSTRACT FROM AUTHOR]

In this paper, the Rayleigh-Lindley (RL) distribution is introduced, obtained by compounding the Rayleigh and Lindley discrete distributions, where the compounding procedure follows an approach similar to the one previously studied by Adamidis and Loukas in some other contexts. The resulting distribution is a two-parameter model, which is competitive with other parsimonious models such as the gamma and Weibull distributions. We study some properties of this new model such as the moments and the mean residual life. The estimation was approached via EM algorithm. The behavior of these estimators was studied in finite samples through a simulation study. Finally, we report two real data illustrations in order to show the performance of the proposed model versus other common two-parameter models in the literature. The main conclusion is that the model proposed can be a valid alternative to other competing models well established in the literature. [ABSTRACT FROM AUTHOR]

In [Martínez-Flórez G, Bolfarine H, Gómez HW. Skew-normal alpha-power model. Statistics. 2014;48(6):1414-1428] the authors make an error in calculating the information matrix of the skew-normal alpha-power model under the normality hypothesis . Here we prove that the information matrix under is singular and we proposed a reparametrization that allows one to calculate the asymptotic distribution of the maximum likelihood estimator. [ABSTRACT FROM AUTHOR]

A new family of slash distributions, the modified slashed-Rayleigh distribution, is proposed and studied. This family is an extension of the ordinary Rayleigh distribution, being more flexible in terms of distributional kurtosis. It arises as a quotient of two independent random variables, one being a Rayleigh distribution in the numerator and the other a power of the exponential distribution in denominator. We present properties of the proposed family. In addition, we carry out estimation of the model parameters by moment and maximum likelihood methods. Finally, we conduct a small-scale simulation study to evaluate the performance of the maximum likelihood estimators and apply the results to a real data set, revealing its good performance. [ABSTRACT FROM AUTHOR]

This article proposes new symmetric and asymmetric distributions applying methods analogous as the ones in Kim (2005) and Arnold et al. (2009) to the exponentiated normal distribution studied in Durrans (1992), that we call the power-normal (PN) distribution. The proposed bimodal extension, the main focus of the paper, is called the bimodal power-normal model and is denoted byBPN(α) model, where α is the asymmetry parameter. The authors give some properties including moments and maximum likelihood estimation. Two important features of the model proposed is that its normalizing constant has closed and simple form and that the Fisher information matrix is nonsingular, guaranteeing large sample properties of the maximum likelihood estimators. Finally, simulation studies and real applications reveal that the proposed model can perform well in both situations. [ABSTRACT FROM AUTHOR]

This paper focuses on the development of a new extension of the generalized skew-normal distribution introduced in Gómez et al. [Generalized skew-normal models: properties and inference. Statistics. 2006;40(6):495–505]. To produce the generalization a new parameter is introduced, the signal of which has the flexibility of yielding unimodal as well as bimodal distributions. We study its properties, derive a stochastic representation and state some expressions that facilitate moments derivation. Maximum likelihood is implemented via the EM algorithm which is based on the stochastic representation derived. We show that the Fisher information matrix is singular and discuss ways of getting round this problem. An illustration using real data reveals that the model can capture well special data features such as bimodality and asymmetry. [ABSTRACT FROM PUBLISHER]

In this paper, the destructive negative binomial (DNB) cure rate model with a latent activation scheme [V. Cancho, D. Bandyopadhyay, F. Louzada, and B. Yiqi, The DNB cure rate model with a latent activation scheme, Statistical Methodology 13 (2013b), pp. 48–68] is extended to the case where the observations are grouped into clusters. Parameter estimation is performed based on the restricted maximum likelihood approach and on a Bayesian approach based on Dirichlet process priors. An application to a real data set related to a sealant study in a dentistry experiment is considered to illustrate the performance of the proposed model. [ABSTRACT FROM AUTHOR]

In this article, we introduce a new extension of the Birnbaum–Saunders (BS) distribution as a follow-up to the family of skew-flexible-normal distributions. This extension produces a family of BS distributions including densities that can be unimodal as well as bimodal. This flexibility is important in dealing with positive bimodal data, given the difficulties experienced by the use of mixtures of distributions. Some basic properties of the new distribution are studied including moments. Parameter estimation is approached by the method of moments and also by maximum likelihood, including a derivation of the Fisher information matrix. Three real data illustrations indicate satisfactory performance of the proposed model. [ABSTRACT FROM AUTHOR]

In this paper, a new survival cure rate model is introduced considering the Yule–Simon distribution [12] to model the number of concurrent causes. We study some properties of this distribution and the model arising when the distribution of the competing causes is the Weibull model. We call this distribution the Weibull–Yule–Simon distribution. Maximum likelihood estimation is conducted for model parameters. A small scale simulation study is conducted indicating satisfactory parameter recovery by the estimation approach. Results are applied to a real data set (melanoma) illustrating the fact that the model proposed can outperform traditional alternative models in terms of model fitting. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM PUBLISHER]

In this paper, we introduce an extension of the generalized exponential (GE) distribution, making it more robust against possible influential observations. The new model is defined as the quotient between a GE random variable and a beta-distributed random variable with one unknown parameter. The resulting distribution is a distribution with greater kurtosis than the GE distribution. Probability properties of the distribution such as moments and asymmetry and kurtosis are studied. Likewise, statistical properties are investigated using the method of moments and the maximum likelihood approach. Two real data analyses are reported illustrating better performance of the new model over the GE model. [ABSTRACT FROM AUTHOR]

In this paper, we study the statistical inference based on the Bayesian approach for regression models with the assumption that independent additive errors follow normal, Student-t, slash, contaminated normal, Laplace or symmetric hyperbolic distribution, where both location and dispersion parameters of the response variable distribution include nonparametric additive components approximated byB-splines. This class of models provides a rich set of symmetric distributions for the model error. Some of these distributions have heavier or lighter tails than the normal as well as different levels of kurtosis. In order to draw samples of the posterior distribution of the interest parameters, we propose an efficient Markov Chain Monte Carlo (MCMC) algorithm, which combines Gibbs sampler and Metropolis–Hastings algorithms. The performance of the proposed MCMC algorithm is assessed through simulation experiments. We apply the proposed methodology to a real data set. The proposed methodology is implemented in the R package BayesGESM using the functiongesm(). [ABSTRACT FROM PUBLISHER]

In this paper, an algorithm is developed to compute estimates for parameters in destructive weighted Poisson cure rate models. It is shown, analytically, the robustness of the procedure with respect to the maximization of the observed likelihood function. The approach allows a simple implementation of different distributions for non-destroyed concurrent causes. The performance of the method is evaluated through simulation studies and by analysing a real data set related to patients with malignant melanoma. [ABSTRACT FROM AUTHOR]

In this article, we consider the inclusion of random effects in both the survival function for at-risk subjects and the cure probability assuming a bivariate normal distribution for those effects in each cluster. For parameter estimation, we implemented the restricted maximum likelihood (REML) approach. We consider Weibull and Piecewise Exponential distributions to model the survival function for non-cured individuals. Simulation studies are performed, and based on a real database we evaluate the performance of our proposed model. Effect of different follow-up times and the effect of considering independent random effects instead of bivariate random effects are also studied. [ABSTRACT FROM PUBLISHER]

As is the case of many studies, the data collected are limited and an exact value is recorded only if it falls within an interval range. Hence, the responses can be either left, interval or right censored. Linear (and nonlinear) regression models are routinely used to analyze these types of data and are based on normality assumptions for the errors terms. However, those analyzes might not provide robust inference when the normality assumptions are questionable. In this article, we develop a Bayesian framework for censored linear regression models by replacing the Gaussian assumptions for the random errors with scale mixtures of normal (SMN) distributions. The SMN is an attractive class of symmetric heavy-tailed densities that includes the normal, Student-t, Pearson type VII, slash and the contaminated normal distributions, as special cases. Using a Bayesian paradigm, an efficient Markov chain Monte Carlo algorithm is introduced to carry out posterior inference. A new hierarchical prior distribution is suggested for the degrees of freedom parameter in the Student-tdistribution. The likelihood function is utilized to compute not only some Bayesian model selection measures but also to develop Bayesian case-deletion influence diagnostics based on theq-divergence measure. The proposed Bayesian methods are implemented in the R packageBayesCR. The newly developed procedures are illustrated with applications using real and simulated data. [ABSTRACT FROM AUTHOR]

In this paper, we propose the use of Bayesian quantile regression for the analysis of proportion data. We also consider the case when the data present a zero-or-one inflation using a two-part model approach. For the latter scheme, we assume that the response variable is generated by a mixed discrete–continuous distribution with a point mass at zero or one. Quantile regression is then used to explain the conditional distribution of the continuous part between zero and one, while the mixture probability is also modelled as a function of the covariates. We check the performance of these models with two simulation studies. We illustrate the method with data about the proportion of households with access to electricity in Brazil. [ABSTRACT FROM AUTHOR]

In recent years, there has been considerable interest in regression models based on zero-inflated distributions. These models are commonly encountered in many disciplines, such as medicine, public health, and environmental sciences, among others. The zero-inflated Poisson (ZIP) model has been typically considered for these types of problems. However, the ZIP model can fail if the non-zero counts are overdispersed in relation to the Poisson distribution, hence the zero-inflated negative binomial (ZINB) model may be more appropriate. In this paper, we present a Bayesian approach for fitting the ZINB regression model. This model considers that an observed zero may come from a point mass distribution at zero or from the negative binomial model. The likelihood function is utilized to compute not only some Bayesian model selection measures, but also to develop Bayesian case-deletion influence diagnostics based onq-divergence measures. The approach can be easily implemented using standard Bayesian software, such as WinBUGS. The performance of the proposed method is evaluated with a simulation study. Further, a real data set is analyzed, where we show that ZINB regression models seems to fit the data better than the Poisson counterpart. [ABSTRACT FROM PUBLISHER]

Skew scale mixtures of normal distributions are often used for statistical procedures involving asymmetric data and heavy-tailed. The main virtue of the members of this family of distributions is that they are easy to simulate from and they also supply genuine expectation-maximization (EM) algorithms for maximum likelihood estimation. In this paper, we extend the EM algorithm for linear regression models and we develop diagnostics analyses via local influence and generalized leverage, following Zhu and Lee's approach. This is because Cook's well-known approach cannot be used to obtain measures of local influence. The EM-type algorithm has been discussed with an emphasis on the skew Student-t-normal, skew slash, skew-contaminated normal and skew power-exponential distributions. Finally, results obtained for a real data set are reported, illustrating the usefulness of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, authors study properties and inference for the newly introduced skew-normal alpha-power model, generalizing both, the power-normal and skew-normal models. Inference is approached via maximum likelihood. Fisher information matrix is derived and shown to be nonsingular at the whole parametric space. Special emphasis is placed on the special case of the power–skew-normal model. Studies with real data illustrate the fact that the model can be very useful in applications, being able to overfit less general models entertained in the literature. [ABSTRACT FROM PUBLISHER]

Data on fatigue to exceed a critical value (or to grow to a critical level at which failure is likely to occur) is typically adjusted using the Birnbaum–Saunders (BS) distribution [see Birnbaum ZW, Saunders SC. A new family of life distributions. J Appl Probab. 1969a;6:319–327]. Although this type of distribution is asymmetric, in some cases the degree of skewness and/or kurtosis are outside the distributional range allowed by the BS distribution. Thus, a more adequate distribution model for better adjusting such unexpected deviations is called for. With this in mind, the main object of this paper is to propose an extension of the BS distribution based on the asymmetric alpha-power family of distributions [see Pewsey A, Gómez HW, Bolfarine H. Likelihood-based inference for power distributions. Test. 2012;21(4):775–789]. This extension is called the exponentiated BS distribution. We expect that by replacing the normal distribution by such more general family, a more flexible BS family is obtained. Asymmetry in the alpha-power family is controlled by a shape parameter, which also presents a similar performance in the extended BS family. The paper presents the density function for the extended BS and derives closed-form expressions for moments. Estimation is dealt with by using maximum likelihood estimators. Large sample inference can be conducted by using the Fisher information matrix derived in the paper. Estimation performance is studied by using a small scale simulation study. Results of a real application illustrates the good performance of the proposed approach. [ABSTRACT FROM PUBLISHER]

Very often, in psychometric research, as in educational assessment, it is necessary to analyze item response from clustered respondents. The multiple group item response theory (IRT) model proposed by Bock and Zimowski [12] provides a useful framework for analyzing such type of data. In this model, the selected groups of respondents are of specific interest such that group-specific population distributions need to be defined. The usual assumption for parameter estimation in this model, which is that the latent traits are random variables following different symmetric normal distributions, has been questioned in many works found in the IRT literature. Furthermore, when this assumption does not hold, misleading inference can result. In this paper, we consider that the latent traits for each group follow different skew-normal distributions, under the centered parameterization. We named it skew multiple group IRT model. This modeling extends the works of Azevedoet al. [4], Bazánet al. [11] and Bock and Zimowski [12] (concerning the latent trait distribution). Our approach ensures that the model is identifiable. We propose and compare, concerning convergence issues, two Monte Carlo Markov Chain (MCMC) algorithms for parameter estimation. A simulation study was performed in order to evaluate parameter recovery for the proposed model and the selected algorithm concerning convergence issues. Results reveal that the proposed algorithm recovers properly all model parameters. Furthermore, we analyzed a real data set which presents asymmetry concerning the latent traits distribution. The results obtained by using our approach confirmed the presence of negative asymmetry for some latent trait distributions. [ABSTRACT FROM AUTHOR]

In this paper we introduce a modified slash distribution obtained by modifying the usual slash distribution. This new distribution is based on the quotient of two independent random variables, whose distributions are the normal and the power of an exponential distribution of scale parameter equals to two, respectively. In this way, the result is a new distribution whose kurtosis values are greater when compared with that of the slash distribution. We study the density, some properties, moments, kurtosis and make inferences by the method of moments and maximum likelihood. We introduce a multivariate version of this new distribution. Moreover, we provide two illustrations with real data showing that the new distribution fits better the data than the ordinary slash distribution. [ABSTRACT FROM AUTHOR]

In this work, we generalize the controlled calibration model by assuming replication on both variables. Likelihood-based methodology is used to estimate the model parameters and the Fisher information matrix is used to construct confidence intervals for the unknown value of the regressor variable. Further, we study the local influence diagnostic method which is based on the conditional expectation of the complete-data log-likelihood function related to the EM algorithm. Some useful perturbation schemes are discussed. A simulation study is carried out to assess the effect of the measurement error on the estimation of the parameter of interest. This new approach is illustrated with a real data set. [ABSTRACT FROM PUBLISHER]

The main object of this article is to propose an extension of the tobit model for which the error distribution follows the power-normal distribution (Gupta and Gupta, 2008). Inference is dealt with by using the likelihood approach. Simulation studies and application to a real data set are used to demonstrate the usefulness of the extension. [ABSTRACT FROM PUBLISHER]

In this paper we study the cure rate survival model involving a competitive risk structure with missing categorical covariates. A parametric distribution that can be written as a sequence of one-dimensional conditional distributions is specified for the missing covariates. We consider the missing data at random situation so that the missing covariates may depend only on the observed ones. Parameter estimates are obtained by using the EM algorithm via the method of weights. Extensive simulation studies are conducted and reported to compare estimates efficiency with and without missing data. As expected, the estimation approach taking into consideration the missing covariates presents much better efficiency in terms of mean square errors than the complete case situation. Effects of increasing cured fraction and censored observations are also reported. We demonstrate the proposed methodology with two real data sets. One involved the length of time to obtain a BS degree in Statistics, and another about the time to breast cancer recurrence. [ABSTRACT FROM AUTHOR]

In this paper, an alternative skew Student-t family of distributions is studied. It is obtained as an extension of the generalized Student-t (GS-t) family introduced by McDonald and Newey [10]. The extension that is obtained can be seen as a reparametrization of the skewed GS-t distribution considered by Theodossiou [14]. A key element in the construction of such an extension is that it can be stochastically represented as a mixture of an epsilon-skew-power-exponential distribution [1] and a generalized-gamma distribution. From this representation, we can readily derive theoretical properties and easy-to-implement simulation schemes. Furthermore, we study some of its main properties including stochastic representation, moments and asymmetry and kurtosis coefficients. We also derive the Fisher information matrix, which is shown to be nonsingular for some special cases such as when the asymmetry parameter is null, that is, at the vicinity of symmetry, and discuss maximum-likelihood estimation. Simulation studies for some particular cases and real data analysis are also reported, illustrating the usefulness of the extension considered. [ABSTRACT FROM AUTHOR]

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353–365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R.B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362–1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis–Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271–275]. Our algorithm has only one Metropolis–Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146–178; R.J. Patz and B.W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342–366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599–607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models. [ABSTRACT FROM PUBLISHER]

In this article, we introduce an asymmetric extension to the univariate slash-elliptical family of distributions studied in Gómez et al. (2007a). This new family results from a scale mixture between the epsilon-skew-symmetric family of distributions and the uniform distribution. A general expression is presented for the density with special cases such as the normal, Cauchy, Student-t, and Pearson type II distributions. Some special properties and moments are also investigated. Results of two real data sets applications are also reported, illustrating the fact that the family introduced can be useful in practice. [ABSTRACT FROM AUTHOR]

In this article, we study further properties of a skew normal distribution, called the skew-normal-Cauchy (SNC) distribution by Nadarajah and Kotz (2003). A stochastic representation is obtained which allows alternative derivations for moments, moments generating function, and skewness and kurtosis coefficients. Issues related to singularity of the Fisher information matrix are investigated. [ABSTRACT FROM AUTHOR]

In general, the normal distribution is assumed for the surrogate of the true covariates in the classical measurement error model. This paper considers a class of distributions, which includes the normal one, for the variables subject to error. An estimation approach yielding consistent estimators is developed and simulation studies reported. [ABSTRACT FROM AUTHOR]

In this article, we study some results related to a specific class of distributions, called skew-curved-symmetric family of distributions that depends on a parameter controlling the skewness and kurtosis at the same time. Special elements of this family which are studied include symmetric and well-known asymmetric distributions. General results are given for the score function and the observed information matrix. It is shown that the observed information matrix is always singular for some special cases. We illustrate the flexibility of this class of distributions with an application to a real dataset on characteristics of Australian athletes. [ABSTRACT FROM AUTHOR]

We review several asymmetrical links for binary regression models and present a unified approach for two skew-probit links proposed in the literature. Moreover, under skew-probit link, conditions for the existence of the ML estimators and the posterior distribution under improper priors are established. The framework proposed here considers two sets of latent variables which are helpful to implement the Bayesian MCMC approach. A simulation study to criteria for models comparison is conducted and two applications are made. Using different Bayesian criteria we show that, for these data sets, the skew-probit links are better than alternative links proposed in the literature. [ABSTRACT FROM AUTHOR]

In this article, we introduce a semi-parametric Bayesian approach based on Dirichlet process priors for the discrete calibration problem in binomial regression models. An interesting topic is the dosimetry problem related to the dose-response model. A hierarchical formulation is provided so that a Markov chain Monte Carlo approach is developed. The methodology is applied to simulated and real data. [ABSTRACT FROM AUTHOR]

Influence diagnostics methods are extended in this article to the Grubbs model when the unknown quantity x (latent variable) follows a skew-normal distribution. Diagnostic measures are derived from the case-deletion approach and the local influence approach under several perturbation schemes. The observed information matrix to the postulated model and Delta matrices to the corresponding perturbed models are derived. Results obtained for one real data set are reported, illustrating the usefulness of the proposed methodology. [ABSTRACT FROM AUTHOR]

In this article, we study a new class of non negative distributions generated by the symmetric distributions around zero. For the special case of the distribution generated using the normal distribution, properties like moments generating function, stochastic representation, reliability connections, and inference aspects using methods of moments and maximum likelihood are studied. Moreover, a real data set is analyzed, illustrating the fact that good fits can result. [ABSTRACT FROM AUTHOR]

In this article, we consider local influence analysis for the skew-normal linear mixed model (SN-LMM). As the observed data log-likelihood associated with the SN-LMM is intractable, Cook's well-known approach cannot be applied to obtain measures of local influence. Instead, we develop local influence measures following the approach of Zhu and Lee (2001). This approach is based on the use of an EM-type algorithm and is measurement invariant under reparametrizations. Four specific perturbation schemes are discussed. Results obtained for a simulated data set and a real data set are reported, illustrating the usefulness of the proposed methodology. [ABSTRACT FROM AUTHOR]

In this article, we compare three residuals based on the deviance component in generalised log-gamma regression models with censored observations. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. For all cases studied, the empirical distributions of the proposed residuals are in general symmetric around zero, but only a martingale-type residual presented negligible kurtosis for the majority of the cases studied. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for the martingale-type residual in generalised log-gamma regression models with censored data. A lifetime data set is analysed under log-gamma regression models and a model checking based on the martingale-type residual is performed. [ABSTRACT FROM AUTHOR]

In this article, we present EM algorithms for performing maximum likelihood estimation for three multivariate skew-normal regression models of considerable practical interest. We also consider the restricted estimation of the parameters of certain important special cases of two models. The methodology developed is applied in the analysis of longitudinal data on dental plaque and cholesterol levels. [ABSTRACT FROM AUTHOR]

Censoring is a common feature in survival data, usually associated with loss to follow-up. However, when the fraction of censored data is high, it may indicate that part of the experimental units are no longer at risk of presenting the event of interest. In this article we consider the approach of Chen et al. (1999) for such situation, and discuss the case where covariates may be measured with error. Simulations and an application to a real dataset are also presented. [ABSTRACT FROM AUTHOR]

The main object of this article is to discuss maximum likelihood inference for the epsilon-skew-t distribution. Special cases of this distribution include the epsilon-skew-Cauchy and the epsilon-skew-normal distributions. We derive the information matrix for the maximum likelihood estimators. The approach is applied to a data set presenting significant amount of skewness and heavy tails. In the application we consider the epsilon-skew-t distribution with known and unknown degrees of freedom parameter, showing great flexibility in adjusting to skew data with heavy tails. [ABSTRACT FROM AUTHOR]

In this paper, we extend the structural probit measurement error model by considering that the unobserved covariate follows a skew-normal distribution. The new model is termed the structural skew-normal probit model. As in the normal case, the likelihood function is obtained analytically which can be maximized by using existing statistical software. A Bayesian approach using Markov chain Monte Carlo techniques to generate from the posterior distributions is also developed. A simulation study demonstrates the usefulness of the approach in avoiding attenuation which is the case with the naive procedure and it seems to be more efficient than using the structural probit model when the distribution of the covariate (predictor) is skew. [ABSTRACT FROM AUTHOR]

In this paper we propose the use of a multivariate null intercept measurement error model, where the true unobserved value of the covariate follows a mixture of two normal distributions. The proposed model is applied to a dental clinical trial presented in Hadgu and Koch (1999). A Bayesian approach is considered and a Gibbs Sampler is used to perform the computations. [ABSTRACT FROM AUTHOR]

In this article, we introduce a new family of asymmetric distributions, which depends on two parameters namely, α and β, and in the special case where β = 0, the skew-normal (SN) distribution considered by Azzallini [Azzalini, A., 1985, A class of distributions which includes the normal ones. Scandinavian Journal of Statistics, 12, 171–178.] is obtained. Basic properties such as a stochastic representation and the derivation of maximum likelihood and moment estimators are studied. The asymptotic behaviour of both types of estimators is also investigated. Results of a small-scale simulation study is provided illustrating the usefulness of the new model. An application to a real data set is reported showing that it can present better fit than the SN distribution. [ABSTRACT FROM AUTHOR]

The purpose of this article is to investigate hypothesis testing in functional comparative calibration models. Wald type statistics are considered which are asymptotically distributed according to the chi-square distribution. The statistics are based on maximum likelihood, corrected score approach, and method of moment estimators of the model parameters, which are shown to be consistent and asymptotically normally distributed. Results of analytical and simulation studies seem to indicate that the Wald statistics based on the method of moment estimators and the corrected score estimators are, as expected, less efficient than the Wald type statistic based on the maximum likelihood estimators for small n. Wald statistic based on moment estimators are simpler to compute than the other Wald statistics tests and their performance improves significantly as n increases. Comparisons with an alternative F statistics proposed in the literature are also reported. [ABSTRACT FROM AUTHOR]

In this paper we discuss the application of local influence in a measurement error regression model with null intercepts under a Student_t model with dependent populations. The Student_t distribution is a robust alternative to modelling data sets involving errors with longer than Normal tails. We derive the appropriate matrices for assessing the local influence for different perturbation schemes and use real data as an illustration of the usefulness of the application. [ABSTRACT FROM AUTHOR]

In this paper, we present a Bayesian methodology for modelling accelerated lifetime tests under a stress response relationship with a threshold stress. Both Laplace and MCMC methods are considered. The methodology is described in detail for the case when an exponential distribution is assumed to express the behaviour of lifetimes, and a power law model with a threshold stress is assumed as the stress response relationship. We assume vague but proper priors for the parameters of interest. The methodology is illustrated by a accelerated failure test on an electrical insulation film. [ABSTRACT FROM AUTHOR]

Longitudinal data are of great interest in analysis of clinical trials. In many practical situations the covariate can not be measured precisely and a natural alternative model is the errors-in-variables regression models. In this paper we study a null intercept errors-in-variables regression model with a structure of dependency between the response variables within the same group. We apply the model to real data presented in Hadgu and Koch (Hadgu, A., Koch, G. (1999). Application of generalized estimating equations to a dental randomized clinical trial. J. Biopharmaceutical Statistics 9(1):161–178). In that study volunteers with preexisting dental plaque were randomized to two experimental mouth rinses (A and B) or a control mouth rinse with double blinding. The dental plaque index was measured for each subject in the beginning of the study and at two follow-up times, which leads to the presence of an interclass correlation. We propose the use of a Bayesian approach to model a multivariate null intercept errors-in-variables regression model to the longitudinal data. The proposed Bayesian approach accommodates the correlated measurements and incorporates the restriction that the slopes must lie in the (0, 1) interval. A Gibbs sampler is used to perform the computations. [ABSTRACT FROM AUTHOR]

This article discusses a Bayesian analysis of repeated measures pretest/posttest data under null intercepts errors-in-variables regression models. For illustration we consider an example in the field of dentistry involving the comparison of two types of toothbrushes with respect to the efficacy in removing dental plaque. The proposed Bayesian approach accommodates the correlated measurements and incorporates the restriction that the slopes must lie in the [0, 1] interval, a feature not considered in the analysis conducted by Singer & Andrade (1997). The observed values of the (repeated) response and explanatory variables are supposed to follow a Multivariate Student-t distribution. A Gibbs sampler is used to perform the computations. [ABSTRACT FROM AUTHOR]