## Trondheim, Norway. On a Hybrid Data Cloning Method and Its Application in Generalized Linear Mixed Models

### BibTeX

@MISC{Baghishani_trondheim,norway.,

author = {Hossein Baghishani and Mohsen Mohammadzadeh and Hossein Baghishani A and H˚avard Rue B and Mohsen Mohammadzadeh A},

title = {Trondheim, Norway. On a Hybrid Data Cloning Method and Its Application in Generalized Linear Mixed Models},

year = {}

}

### OpenURL

### Abstract

Data cloning method is a new computational tool for computing maximum likelihood estimates in complex statistical models such as mixed models. The data cloning method is synthesized with integrated nested Laplace approximation to compute maximum likelihood estimates efficiently via a fast implementation in generalized linear mixed models. Asymptotic normality of the hybrid data cloning based distribution is established aided by modification of Stein’s Identity. The results are illustrated through a series of well known examples. It is shown that the proposed method as well as normal approximation perform very well and justify the theory.

### Citations

1807 |
Generalised Linear Modelling
- McCullagh, Nelder
- 1989
(Show Context)
Citation Context ...ces such as biology, ecology, epidemiology and medicine. The Generalized Linear Mixed Models are flexible models for modeling these types of data. As an extension of generalized linear models (GLMs) (=-=McCullagh and Nelder, 1989-=-), a GLMM assumes that the response variable follows a distribution from the exponential family and is conditionally independent given latent variables, while the latent variables are modeled by rando... |

358 |
Approximate inference in generalized linear mixed models
- Breslow, Clayton
- 1993
(Show Context)
Citation Context ...iable follows a distribution from the exponential family and is conditionally independent given latent variables, while the latent variables are modeled by random effects that are typically Gaussian (=-=Breslow and Clayton, 1993-=-). ∗ To whom correspondence should be addressed. Email addresses: hbaghishani@modares.ac.ir, hossein.baghishani@math.ntnu.no (Hossein Baghishani), hrue@math.ntnu.no (H˚avard Rue), mohsen m@modares.ac.... |

90 | Generalized linear models with random effects: A Gibbs sampling approach - Zeger, Karim - 1991 |

88 | Approximte bayesian inference for latent gaussian models using integrated nested laplace approximations (with discussion
- Rue, Martino, et al.
- 2009
(Show Context)
Citation Context ... main aim is to compute the posterior marginals π(ψl|y), l=1,...,q and π(θv|y), v = 1,...,d. It is well known, however, that MCMC methods tend to exhibit poor performance when applied to such models (=-=Rue et al., 2009-=-). INLA is a new tool for Bayesian inference on latent Gaussian models introduced by Rue et al. (2009). The method combines Laplace approximations and numerical integration in a very efficient manner.... |

85 |
Maximum likelihood algorithms for generalized linear mixed models
- McCulloch
- 1997
(Show Context)
Citation Context ...30, 2010Statistical inferences in such models has been the subject of a great deal of research over the past decade. Both frequentist and Bayesian methods have been developed for inference in GLMMs (=-=McCulloch, 1997-=-). Due to the advances in computation, the most commonly used approach for inference in these models is based on the Bayesian paradigm, especially Markov chain Monte Carlo (MCMC) algorithms. But Bayes... |

40 | Some covariance models for longitudinal count data with overdispersion - Thall, Vail - 1990 |

38 |
Beta-binomial Anova for proportions
- Crowder
- 1978
(Show Context)
Citation Context ...ion This example concerns data on the proportion of seeds that germinated on each of m =21 plates arranged according to a 2 × 2 factorial design with respect to seed variety and type of root extract (=-=Crowder, 1978-=-). The sampling model is Yi|β,pi ∼ Bin(ni,pi) where, for plate i, yi is the number of germinating seeds and ni is the total number of seeds for i =1,...,m. To account for the extraneous between plate ... |

20 | Approximate Bayesian inference for hierarchical Gaussian Markov random field models - Rue, Martino - 2007 |

14 | A survey of Monte Carlo algorithms for maximizing the likelihood of a two-stage hierarchical model. Statistical Modelling
- Booth, Hobert, et al.
- 2001
(Show Context)
Citation Context ...tational time required and convergence assessment, e.g. penalized quasi likelihood (Breslow and Clayton, 1993), composite likelihood (Bellio and Varin, 2005) and Monte Carlo expectation maximization (=-=Booth et al., 2001-=-). 23A disadvantage of the our work is that, according to INLA methodology, the prior distributions for fixed effects β must be Gaussian. However, we can use Gaussian priors with high variances to co... |

8 |
A Pairwise Likelihood Approach to Generalized Linear Models with Crossed Random Effects. Statistical Modelling 5(3
- Bellio, Varin
- 2005
(Show Context)
Citation Context ...s in the sense of consistent estimation, loss of efficiency, computational time required and convergence assessment, e.g. penalized quasi likelihood (Breslow and Clayton, 1993), composite likelihood (=-=Bellio and Varin, 2005-=-) and Monte Carlo expectation maximization (Booth et al., 2001). 23A disadvantage of the our work is that, according to INLA methodology, the prior distributions for fixed effects β must be Gaussian.... |

8 | Data cloning: easy maximum likelihood estimation for complex ecological models using Bayesian Markov chain Monte Carlo methods - Lele, Dennis, et al. - 2007 |

6 |
Bayesian inference for Generalized Linear Mixed Models
- Fong, Rue, et al.
- 2009
(Show Context)
Citation Context ... main aim is to compute the posterior marginals π(ψl|y), l=1,...,q and π(θv|y), v = 1,...,d. It is well known, however, that MCMC methods tend to exhibit poor performance when applied to such models (=-=Rue et al., 2009-=-). INLA is a new tool for Bayesian inference on latent Gaussian models introduced by Rue et al. (2009). The method combines Laplace approximations and numerical integration in a very efficient manner.... |

2 |
A data cloning algorithm for computing maximum likelihood in spatial generalized linear mixed models
- Baghishani, Mohammadzadeh
- 2009
(Show Context)
Citation Context ...terior distribution, but it is constructed from two functions which are not in fact a prior distribution and a likelihood. However, considering them as prior and likelihood can be mimicked virtually (=-=Baghishani and Mohammadzadeh, 2009-=-). The trick in the DC is generating samples from a DC-based distribution constructed by duplicating the original data set enough times, k say, such that the sample mean as well as the scaled sample v... |

2 | Hierarchical models in ecology: confidence intervals, hypothesis testing, and model selection using data cloning - Ponciano, Taper, et al. - 2009 |

2 |
Asymptotic posterior normality for multiparameter problems, Journal of statistical planning and inference
- Weng, Tsai
- 2008
(Show Context)
Citation Context ... ℜr−j ∫ ∞ bj {hj(b1,...,bj−1,c) − hj−1(b1,...,bj−1)}e 1 − 2 e2 for −∞ <b1,...,br < ∞ and j =1,...,r.ThenletUh =(g1,...,gr) T . Following lemma states a modified version of Stein’s Identity. Lemma 1. (=-=Weng and Tsai, 2008-=-). Let s be a nonnegative integer and let dΓ =fdΦr, wheref is differentiable on ℜr such that ∫ |f|dΦr + ℜ r ∫ ℜ r (1 + ‖a‖ s )‖∇f(a)‖Φr(da) < ∞. dc, Then for all h ∈ Hs. ∫ Γh − Γ1 · Φrh = (Uh(a)) T ∇f... |

1 | Asymptotic normality of posterior distributions for generalized linear mixed models - Baghishani, Mohammadzadeh - 2010 |

1 | Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations - Tierney, Kadane - 2009 |