參考方法:
目的
lsmeans生成的結(jié)果不是數(shù)據(jù)框没隘,用str也沒(méi)有辦法提取贰拿,怎么提取結(jié)果烁落?
方法
代碼:
require(lsmeans)
fiber.lm <- lm(strength ~ diameter + machine, data = fiber)
fiber.lsm <- lsmeans (fiber.lm, "machine")
fiber.lsm
結(jié)果:
machine lsmean SE df lower.CL upper.CL
A 40.38241 0.7236252 11 38.78972 41.97510
B 41.41922 0.7444169 11 39.78077 43.05767
C 38.79836 0.7878785 11 37.06426 40.53247
Confidence level used: 0.95
提取方法
使用summary,提取之后是數(shù)據(jù)框膜赃,可以直接保存成csv形式挺邀。
> (ls <- summary(fiber.lsm))
machine lsmean SE df lower.CL upper.CL
A 40.38241 0.7236252 11 38.78972 41.97510
B 41.41922 0.7444169 11 39.78077 43.05767
C 38.79836 0.7878785 11 37.06426 40.53247
Confidence level used: 0.95
write.csv(ls,"d:/ls.csv",row.names = F)
結(jié)果
result