问问太笼统了,最好可以给出可以重复的代码,你想要哪个,这样大家才会有目的的帮助。
以?lm的例子为例:
ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
你要的标准误是不是这个:
as.data.frame(summary(lm.D9)$coefficients[,2])


雷达卡



京公网安备 11010802022788号







