Accent on Teaching Materials

Estimation of covariance components in the multivariate random-effect model with nested covariance structure is discussed. There are two covariance matrices to be estimated, namely, the between-group and the within-group covariance matrices. These two covariance matrices are most often estimated by forming a multivariate analysis of variance and equating mean square matrices to their expectations. Such a procedure involves taking the difference between the between-group mean square and the within-group mean square matrices, and often produces an estimated between-group covariance matrix that is not nonnegative definite. We present estimators of the two covariance matrices that are always proper covariance matrices. The estimators are the restricted maximum likelihood estimators if the random effects are normally distributed. The estimation procedure is extended to more complicated models, including the twofold nested and the mixed-effect models. A numerical example is presented to illustrate the use of the estimation procedure.     

The power series skew normal class of distributions

A class of power series skew normal distributions is introduced by generalizing the geometric skew normal distribution of Kundu. Various mathematical properties are derived and estimation addressed by the method of maximum likelihood. The data application of Kundu [Sankhyā B, 76, 2014, 167–189] is revisited and the proposed class is shown to provide a better fit.     

Reply

Selection processes that are inherent in litigation complicate policy-oriented research of medical malpractice tort reforms. In regard to their deterrent impact, the range of potential inferences based on analyses of claim frequency is limited because plaintiffs only file a subset of potential claims. In regard to their impact on litigation costs, researchers often analyze data on claim disposition, but it is difficult to determine whether effects attributed to tort reforms reflect changes in litigant behavior or their influence on the selection of claims. In this article, we evaluate these problems and report results of our study of the effects of medical malpractice reforms on claim disposition.     

The Extended Skew-Exponential Power Distribution and Its Derivation

We consider an extended family of asymmetric univariate distributions generated using a symmetric density, f, and the cumulative distribution function, G, of a symmetric distribution, which depends on two real-valued parameters λ and β and is such that when β = 0 it includes the entire class of distributions with densities of the form g(z | λ) = 2 G(λ z) f(z). A key element in the construction of random variables distributed according to the family is that they can be represented stochastically as the product of two random variables. From this representation we can readily derive theoretical properties, easy-to-implement simulation schemes, as well as extensions to the multivariate case and an explicit procedure for obtaining the moments. We give special attention to the extended skew-exponential power distribution. We derive its information matrix in order to obtain the asymptotic covariance matrix of the maximum likelihood estimators. Finally, an application to a real data set is reported, which shows that the extended skew-exponential power model can provide a better fit than the skew-exponential power distribution.     

Fourth Quarter Survey of the Economic Outlook

Written mainly for its pedagogical interest, this note deals with the computational formulas for the recursive updating of weighted least squares parameter estimates and the residual sum of squares in the general linear model under the assumption that the errors have a multivariate normal distribution. This approach simplifies considerably the derivations of Haslett (1985).     

Maximum likelihood estimation of parameter structures for the Wishart distribution using constraints

Maximum likelihood estimation under constraints for estimation in the Wishart class of distributions, is considered. It provides a unified approach to estimation in a variety of problems concerning covariance matrices. Virtually all covariance structures can be translated to constraints on the covariances. This includes covariance matrices with given structure such as linearly patterned covariance matrices, covariance matrices with zeros, independent covariance matrices and structurally dependent covariance matrices. The methodology followed in this paper provides a useful and simple approach to directly obtain the exact maximum likelihood estimates. These maximum likelihood estimates are obtained via an estimation procedure for the exponential class using constraints.     

Classification Rules for Multivariate Repeated Measures Data with Equicorrelated Correlation Structure on Both Time and Spatial Repeated Measurements

We study the problem of classification for multivariate repeated measures data with structured correlations on both time and spatial repeated measurements. This is a very important problem in many biomedical as well as in engineering field. Classification rules as well as the algorithm to compute the maximum likelihood estimates of the required parameters are given.     

Acceleration of the EM Algorithm by using Quasi-Newton Methods

The EM algorithm is a popular method for maximum likelihood estimation. Its simplicity in many applications and desirable convergence properties make it very attractive. Its sometimes slow convergence, however, has prompted researchers to propose methods to accelerate it. We review these methods, classifying them into three groups: pure , hybrid and EM-type accelerators. We propose a new pure and a new hybrid accelerator both based on quasi-Newton methods and numerically compare these and two other quasi-Newton accelerators. For this we use examples in each of three areas: Poisson mixtures, the estimation of covariance from incomplete data and multivariate normal mixtures. In these comparisons, the new hybrid accelerator was fastest on most of the examples and often dramatically so. In some cases it accelerated the EM algorithm by factors of over 100. The new pure accelerator is very simple to implement and competed well with the other accelerators. It accelerated the EM algorithm in some cases by factors of over 50. To obtain standard errors, we propose to approximate the inverse of the observed information matrix by using auxiliary output from the new hybrid accelerator. A numerical evaluation of these approximations indicates that they may be useful at least for exploratory purposes.     

Large Sample Properties of ML Estimator of the Parameters of Multivariate O–U Random Fields

Computer modeling is having a profound effect on scientific research, replacing direct physical experimentation by computer simulation of complex models. In this research, the computer output, X(t), is assumed to be a multivariate, three-dimensional (time) Ornstein-Uhlenbeck (O–U) process with parametric covariance function. It is shown that the ML estimates of the parameters are strongly consistent and asymptotically normal when the observations are taken from a complete lattice, not necessarily equally spaced.     

Covariance Matrix Formula for Exponential Family Nonlinear Models

This article gives a matrix formula for second-order covariances of maximum likelihood estimators in exponential family nonlinear models, thus generalizing the result of Cordeiro (2004 Cordeiro , G. M. ( 2004 ). Second-order covariance matrix of maximum likelihood estimates in generalized linear models . Statist. Probab. Lett. 66 : 153 – 160 .[Crossref], [Web of Science ®] , [Google Scholar]) valid for generalized linear models with known dispersion parameter. Some simulations show that the second-order covariances for exponential family nonlinear models can be quite pronounced in small to moderate sample sizes.     

On the Estimation of Parameters of Kumaraswamy-G Distributions

Cordeiro and de Castro proposed a new family of generalized distributions based on the Kumaraswamy distribution (denoted as Kw-G). Nadarajah et al. showed that the density function of the new family of distributions can be expressed as a linear combination of the density of exponentiated family of distributions. They derived some properties of Kw-G distributions and discussed estimation of parameters using the maximum likelihood (ML) method. Cheng and Amin and Ranneby introduced a new method of estimating parameters based on Kullback–Leibler divergence (the maximum spacing (MSP) method). In this article, the estimates of parameters of Kw-G distributions are obtained using the MSP method. For some special Kw-G distributions, the new estimators are compared with ML estimators. It is shown by simulations and a real data application that MSP estimators have better properties than ML estimators.     

The Minimum Density Power Divergence Estimation for the Lognormal Density

In this article, we implement the minimum density power divergence estimation for estimating the parameters of the lognormal density. We compare the minimum density power divergence estimator (MDPDE) and the maximum likelihood estimator (MLE) in terms of robustness and asymptotic distribution. The simulations and an example indicate that the MDPDE is less biased than MLE and is as good as MLE in terms of the mean square error under various distributional situations.     

Limit expressions for the risk of james-stein estimators

Limit expressions (as the dimension p ← ∞ ) are derived for the relative risk of the James-Stein estimator and its positive-part version. The limit is simple to evaluate, and gives the amount of improvement in risk that is possible. The technique used is to bound the risk, both above and below. with bounds that converge to the same limit. For the James-Stein estimator these bounds are simple to calculate, and are quite accurate even for moderate dimensions.     

Non–zero inflated modified power series distributions

In this paper, we study the class of inflated modified power series distributions (IMPSD) where inflation occurs at any of the support points. This class include among other the generalized Poisson, the generalized negative binomial, the generalized logarithmic series and the lost games distributions. We give expressions for the moments, factorial moments and central moments of the IMPSD. The maximum likelihood estimation of the parameters of the IMPSD and the variance – covariance matrix of the estimators is obtained. We derive these estimators and their information matrices for mentioned above particular members of IMPSD class. The second part of this paper deals with the distribution of sum of independent and identically distributed random variables taking values s, s+1. s + 2, …, s ≥ 0, with modified power series distributions inflated at the point s.     

The Two-Piece Normal-Laplace Distribution

This article considers the two-piece normal-Laplace (TPNL) distribution, a split skew distribution consisting of a normal part, and a Laplace part. The distribution is indexed by three parameters, representing location, scale, and shape. As illustrated with several examples, the TPNL family of distributions provides a useful alternative to other families of asymmetric distributions on the real line. However, because the likelihood function is not well behaved, standard theory of maximum-likelihood (ML) estimation does not apply to the TPNL family. In particular, the likelihood function can have multiple local maxima. We provide a procedure for computing ML estimators, and prove consistency and asymptotic normality of ML estimators, using non standard methods.     

On Jointly Estimating Parameters and Missing Data by Maximizing the Complete-Data Likelihood

One approach to handling incomplete data occasionally encountered in the literature is to treat the missing data as parameters and to maximize the complete-data likelihood over the missing data and parameters. This article points out that although this approach can be useful in particular problems, it is not a generally reliable approach to the analysis of incomplete data. In particular, it does not share the optimal properties of maximum likelihood estimation, except under the trivial asymptotics in which the proportion of missing data goes to zero as the sample size increases.