Mcp2matrix(型号，链接=链接) glht出错

Question

我想对我的数据做一个Tukey测试，找到一些code for it。我花了一些时间将它应用于我的数据，当我到达glht时，我得到了上面的错误。我想可能是我的数据输入错误，所以我用演示数据对其进行了测试，但我得到了相同的错误：

mcp2matrix(model，linfct = linfct)中出错：‘character’类的变量‘X’未作为因子包含在‘model’中。

以下是演示数据：

           Y X
1  0.2855104 A
2  0.3385240 A
3  0.0883132 A
4  0.2059308 A
5  0.3632401 A
6  0.5217391 B
7  0.7633588 B
8  0.3254679 B
9  0.4253057 B
10 0.3780718 B
11 0.9891197 C
12 1.1927181 C
13 0.7882883 C
14 0.5491762 C
15 0.5445882 C
16 1.2670565 D
17 1.6253208 D
18 1.2661090 D
19 1.1541876 D
20 1.2684989 D
21 1.0695187 D

这是完整的代码：

install.packages("multcomp")
expt1 <- read.table("demo.txt",header=T)
amod <- aov(Y~X,data=expt1)
summary(amod)
library("multcomp")
tmod <- glht(amod,linfct=mcp(X="Tukey"))

这是它在R中的显示方式(在安装软件包之后)：

> expt1 <- read.table("demo.txt",header=T)
> amod <- aov(Y~X,data=expt1)
> summary(amod)
            Df Sum Sq Mean Sq F value   Pr(>F)    
X            3  3.244  1.0815   27.59 9.26e-07 ***
Residuals   17  0.666  0.0392                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> library("multcomp")
Loading required package: mvtnorm
Loading required package: survival
Loading required package: TH.data
Loading required package: MASS

Attaching package: ‘TH.data’

The following object is masked from ‘package:MASS’:

    geyser

> tmod <- glht(amod,linfct=mcp(X="Tukey"))
Error in mcp2matrix(model, linfct = linfct) : 
  Variable(s) ‘X’ of class ‘character’ is/are not contained as a factor in ‘model’.

我完全是个新手，完全被难住了。我做错了什么？

Answer 1

试试这个：

data$X <- ordered(data$X, levels = c("A", "B", "C", "D"))