1

orangeslsmeans で提供されるデータを使用しています。

library(lsmeans)
oranges.rg1<-lm(sales1 ~ price1 + price2 + day + store, data = oranges)
days.lsm <- lsmeans(oranges.rg1, "day")
days_contr.lsm <- contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6))

信頼区間は をプロットすることで視覚化できますがplot(contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6)))、days_contr.lsm

> days_contr.lsm
 contrast       estimate       SE df t.ratio p.value
 1 - avg(5,6) -7.8538769 2.194243 23  -3.579  0.0058
 2 - avg(5,6) -6.9234858 2.127341 23  -3.255  0.0125
 3 - avg(5,6)  0.2462789 2.155529 23   0.114  0.9979
 4 - avg(5,6) -4.6760034 2.110761 23  -2.215  0.1184

信頼区間を data.frame に抽出するにはどうすればよいですか?

> days_contr.lsm
 contrast       estimate       SE df t.ratio p.value lower.CL upper.CL
 1 - avg(5,6) -7.8538769 2.194243 23  -3.579  0.0058  ?       ?
 2 - avg(5,6) -6.9234858 2.127341 23  -3.255  0.0125  ?       ?
 3 - avg(5,6)  0.2462789 2.155529 23   0.114  0.9979  ?       ?
 4 - avg(5,6) -4.6760034 2.110761 23  -2.215  0.1184  ?       ?
4

2 に答える 2

2

confint(contrast(days.lsm, "trt.vs.ctrl", ref = c(5,6)))うまくいった

于 2016-07-07T21:17:18.100 に答える