  Generalized linear mixed models (GLMM) encompass large class of statistical
models, with a vast range of applications areas. GLMM extends the linear mixed
models allowing for different types of response variable. Three most common
data types are continuous, counts and binary and standard distributions for
these types of response variables are Gaussian, Poisson and Binomial,
respectively. Despite that flexibility, there are situations where the response
variable is continuous, but bounded, such as rates, percentages, indexes and
proportions. In such situations the usual GLMM's are not adequate because
bounds are ignored and the beta distribution can be used. Likelihood and
Bayesian inference for beta mixed models are not straightforward demanding a
computational overhead. Recently, a new algorithm for Bayesian inference called
INLA (Integrated Nested Laplace Approximation) was proposed.INLA allows
computation of many Bayesian GLMMs in a reasonable amount time allowing
extensive comparison among models. We explore Bayesian inference for beta mixed
models by INLA. We discuss the choice of prior distributions, sensitivity
analysis and model selection measures through a real data set. The results
obtained from INLA are compared with those obtained by an MCMC algorithm and
likelihood analysis. We analyze data from an study on a life quality index of
industry workers collected according to a hierarchical sampling scheme. Results
show that the INLA approach is suitable and faster to fit the proposed beta
mixed models producing results similar to alternative algorithms and with
easier handling of modeling alternatives. Sensitivity analysis, measures of
goodness of fit and model choice are discussed.
