rainning
请教高手指点:
下面是我的一些结果,
mon1<-by(Dd,Dd$mon,function(dat) lm(monrtn~mkrtn+monmgsy+montover+rate+as.factor(cxsid)+as.factor(jx),data=dat))
sapply(mon1,sigma)
结果出来的是:
01 02 03 04
call Expression Expression Expression Expression
terms Expression Expression Expression Expression
residuals Numeric,4997 Numeric,5037 Numeric,5063 Numeric,5112
coefficients Numeric,32 Numeric,32 Numeric,32 Numeric,32
aliased Logical,8 Logical,8 Logical,8 Logical,8
sigma 0.1097819 0.09411328 0.1084145 0.1496733
df Integer,3 Integer,3 Integer,3 Integer,3
r.squared 0.4884857 0.3504603 0.5252814 0.3936332
adj.r.squared 0.487768 0.3495562 0.524624 0.3928016
fstatistic Numeric,3 Numeric,3 Numeric,3 Numeric,3
cov.unscaled Numeric,64 Numeric,64 Numeric,64 Numeric,64
na.action Integer,42 Integer,34 Integer,46 Integer,53
没有标准差,从而不知系数是否显著,不知怎么办才好了!
请高手指点