search for: vsire

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...