在lsmeans
包中,您可以通过默认的lsmobj (ref.grid)
函数绘制plot()
对象。
# Uses supplied dataset 'fiber'
fiber.lm <- lm(strength ~ diameter + machine, data = fiber)
# adjusted means and comparisons, treating machine C as control
fiber.lsm <- lsmeans (fiber.lm, "machine")
plot(fiber.lsm)
使用标准xlab
和ylab
参数自定义x和y轴标签。
plot(fiber.lsm, xlab="Estimated Strength", ylab="Machine Type")
但是,使用horiz=F
将分类变量放在x轴上似乎会抑制xlab
参数(尽管不是ylab
参数)。
plot(fiber.lsm, horiz=F, ylab="Estimated Strength", xlab="Machine Type")
在允许x轴标签定制的同时,如何使用horiz=F
进行绘图的任何想法?这可能是lsmeans
包或底层lattice
包的问题。我很高兴能够完全抑制x轴标签并稍后将其添加到title
,但这似乎不适用于lattice
图。
答案 0 :(得分:1)
嗯,您更喜欢使用自己的功能,但使用ggplot2:
tmp<-data.frame(summary(fiber.lsm))
ggplot(tmp, aes(x=machine, y=lsmean)) +
geom_errorbar(aes(ymin=lower.CL, ymax=upper.CL), width=.1) +
geom_point()+
labs(x="Machine Type", y="Estimated Strength")
我认为默认是您想要的方向,但如果您想翻转图表,请使用+ coord_flip()
答案 1 :(得分:1)
我发现错误,它可以在版本2.13中正常工作 - 也许在11月中旬左右的CRAN上。
> plot(fiber.lsm, horiz=F, ylab="Estimated Strength", xlab="Machine Type")
感谢您报告此事。