lsmeans R包 提取lsmeans的方法

参考方法：

原文链接

目的

lsmeans生成的结果不是数据框，用str也没有办法提取，怎么提取结果？

方法

代码：

require(lsmeans)
fiber.lm <- lm(strength ~ diameter + machine, data = fiber)
fiber.lsm <- lsmeans (fiber.lm, "machine")
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 

提取方法

使用summary，提取之后是数据框，可以直接保存成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)

结果

result