Hi:

I am attempting to extract the RE standard deviations from my glmer object. I've used the following code, but it gives cval "NULL". 

pf1 <- glmer(PfPrev ~ 1 + Age + I(Age^2) + SinWeek + CosWeek + ITN  + (1|vblock) + (1|cblock), family = binomial, data = all.prev)

ChildRE <- VarCorr(pf1)[1]
cval <- attr(ChildRE, "stddev")

I also tried:

ChildRE <- VarCorr(pf1)
cval <- attr(ChildRE, "sc")

Thanks for your help.

Kind regards,

Katie Benton






	[[alternative HTML version deleted]]

