A likelihood reformulation method in non-normal random effects models.
In this paper, we propose a practical computational method to obtain the maximum likelihood estimates (MLE) for mixed models with non-normal random effects. By simply multiplying and dividing a standard normal density, we reformulate the likelihood conditional on the non-normal random effects to that conditional on the normal random effects. Gaussian quadrature technique, conveniently implemented in SAS Proc NLMIXED, can then be used to carry out the estimation process. Our method substantially reduces computational time, while yielding similar estimates to the probability integral transformation method (J. Comput. Graphical Stat. 2006; 15:39-57). Furthermore, our method can be applied to more general situations, e.g. finite mixture random effects or correlated random effects from Clayton copula. Simulations and applications are presented to illustrate our method.
Liu
Lei
L
Division of Biostatistics and Epidemiology, Department of Public Health Sciences, The University of Virginia, Charlottesville, VA 22908-0717, USA. liulei@virginia.edu
Yu
Zhangsheng
Z
D003198:Computer Simulation; D003930:Diabetic Retinopathy / Q000628:therapy; D006801:Humans; D017075:Laser Coagulation; D016013:Likelihood Functions; D016014:Linear Models; D015233:Models, Statistical; D016016:Proportional Hazards Models; D011897:Random Allocation