各位，我是一个生物学专业的学生，现在需要用R分析数量性状Qst，但是不专业。

有人给我如下的Code，但是计算出来没有confidence interval，哪位好心的朋友帮忙一下？

如果有相关专业人士，感兴趣的话可以算合作。本人论文即将发表，但是统计分析不够。希望得到大家的帮忙。

邮件rovial@ymail.com

#Variance components for Qst (1)
library(ape) #load functions
hist(At.pheno.Qst.vern$LTL) #distribution of values; this is needed to decide on data transformations

varcomp(lme((LTL)~1, data=At.pheno.Qst.vern[-c(which(is.na(At.pheno.Qst.vern$LTL))),], random=~1|PopID/FamID), scale=TRUE)
#calculates variance components. The first value, the between population variance component, is the Qst.

[此贴子已经被作者于2008-12-19 18:51:38编辑过]