Displaying 1 result from an estimated 1 matches for "vsire".
Did you mean:
sire
2002 Aug 28
0
Extracting variance component estimates from lme
...omponents
so that I can then calculate a statistic called Qst from them. So I
have this:
reg1 <- lme(y ~ temp*food, random =~1|POP/SIRE/DAM,
na.action=na.exclude))
summary(reg1)
which shows me the random effects as standard deviations. I then want
the SIRE and POP variance components (as Vsire and Vpop respectively) to
put into the calculation of Qst as
Qst <- Vpop/(Vpop + 8*Vsire)
And obviously writing down the components and then plugging them in is
not elegant. Especially as I then want to calculate the standard error
for Qst, which leads me to my supplementary question...
I c...