R 使用对比软件包的lm对比度
我正在使用R 使用对比软件包的lm对比度,r,contrast,R,Contrast,我正在使用contrast包为R中的lm构建对比。通过下面的代码,我对Trt1和Trt3进行了对比 Y <- c(10, 12, 14, 16, 9, 8) Trt <- gl(n=3, k=2, length=3*2) Data1 <- data.frame(Y, Trt) Data1.lm <- lm(Y~Trt, data = Data1) library(contrast) Contrs1 <- contrast(Data1.lm, a=list(Trt
contrast
包为R
中的lm
构建对比。通过下面的代码,我对Trt1和Trt3进行了对比
Y <- c(10, 12, 14, 16, 9, 8)
Trt <- gl(n=3, k=2, length=3*2)
Data1 <- data.frame(Y, Trt)
Data1.lm <- lm(Y~Trt, data = Data1)
library(contrast)
Contrs1 <- contrast(Data1.lm, a=list(Trt="1"), b=list(Trt="3"), type = "average")
print(Contrs1, X=TRUE)
Y您应该将包含的治疗水平指定为向量(Trt=c(“1”、“2”)
),而不是列表。我通过查看?contrast.lm
中的示例了解到了这一点(尽管不可否认,了解您在寻找什么会有所帮助):
contracs2 | t |)
## 1 4.5 1.06066 1.124506 7.875494 4.24 3 0.024
Contrs2 <- contrast(Data1.lm, a=list(Trt="1", Trt="2"), b=list(Trt="3"), type = "average")
print(Contrs2, X=TRUE)
lm model parameter contrast
Contrast S.E. Lower Upper t df Pr(>|t|)
1 6.5 1.224745 2.602315 10.39768 5.31 3 0.0131
Contrast coefficients:
(Intercept) Trt2 Trt3
1 0 1 -1
Contrs2 <- contrast(Data1.lm, a=list(Trt=c("1","2")), b=list(Trt="3"),
type = "average")
print(Contrs2, X=TRUE)
## lm model parameter contrast
## Contrast S.E. Lower Upper t df Pr(>|t|)
## 1 4.5 1.06066 1.124506 7.875494 4.24 3 0.024