*Dear R-help:

My name is Rodrigo and I have a question with nlme package
in R to fit a mixed beta regression model. The details of the model are:

Suppose that:*

*j in {1, ..., J}* *(level 1)*
*i in {1, ..., n_j}* *(level 2)*

*y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij})

y_{ij} = mu_{ij} + w_{ij}
*
*with*

*logit(mu_{ij}) = Beta_{0i} + Beta_{1i} * x1_{ij} + b2 * x2_{ij}
log(phi_{ij}) = Gamma_{0i} + Gamma_{1i} * z1_{ij} + c2 * z2_{ij}
*
*Beta_{0i} = b_0 + u_{0i}
Beta_{1i} = b_1 + u_{1i}
Gamma_{0i} = c_0 + v_{0i}
Gamma_{1i} = c_1 + v_{1i}
*
*The vector* *(u_{0i}, u_{1i})'* *has normal distribution with mean*
*(0, 0)'* *and covariance matrix*

*sigma_{00} sigma_{01}
sigma_{10} sigma_{11}
*
*The vector* *(v_{0i}, v_{1i})'* *has normal distribution with mean*
*(0, 0)'* *and covariance matrix*

*delta_{00} delta_{01}
delta_{10} delta_{11}
*
*The* *w_{ij}s are independents. Each* *w_{ij}* *has mean* *0* *and variance
* *W*.

*The unknown parameters are:*
*
b_0, b_1, b_2, c_0, c_1, c_2, sigma_{00}, sigma_{10},
sigma_{11}, delta_{00}, delta_{10}, **and* *delta_{11}*.

*The question is:
How can I use nlme package in R to fit this model?
If you want to know additional information, send me a mail, please.

In advance, thanks a lot for your help.

Kind regards,
Rodrigo.*

	[[alternative HTML version deleted]]

