Dear all,

 

I was using the lme4 package in R to implement "two-level logistic mixed
model" for my data. In particular, my model is:

 

yij ~ Binomial (nij,pij)

logit(pij)=xij*beta+vi+uij

 

where we have I areas (i=1,.,I) and ni units in each area (j=1,.,ni), and vi
and uij are area-specific and unit level Normal random effects,
respectively.

 

I used glmer(.) and successfully was able to get the results including model
parameter estimates and corresponding log-likelihood. Here is the code:

"glmer(cbind(y, n-y) ~  x+(1 | area)+(1|unit),family = binomial, data =
dat)"

 

 

Now, I am interested in evaluating "different parameters" to get the
corresponding "log-likelihood". In particular, I am interested to see if I
use other "known and fixed parameters" whether I can use lme4 package to get
the corresponding "log-likelihood".

 

Thanks in advance for your time.

 

Sincerely,

Mahmoud Torabi

 


	[[alternative HTML version deleted]]


