參考方法：

原文鏈接

目的

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