r - mcp2matrix(模型,linfct = linfct)中的错误:“整数”类的变量“方法”是/不包含在“模型”中
问题描述
我不明白为什么它在这里不起作用是我的代码:
library(lme4,multcomp)
tijd <- c(7.667, 9.5, 6.467, 10.033, 9.767, 29.133, 15.933, 8.1, 8.6,
12.3, 5.867, 13.667, 12.933, 11.933, 27.2, 11.567, 10.9, 12.567,
9.033, 12.6, 13.133, 17.067, 14.333, 15.2, 9.1, 9.167, 8.4, 11.4,
18.467, 15.633, 10.4, 10.767, 12.7, 7.667, 9.867, 36.167, 20.4,
33.8, 15.767, 20.1, 10.533, 17.7, 62.667, 21.133, 40.133, 9.2,
19.467, 7.667, 8.867, 20.333, 14.367, 10.167, 28.4, 11.367, 10.867,
22.1, 24.8, 13.3, 17.067, 10.967, 12.333, 28.867, 21.1)
week <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3)
Methode <- c(1, 1, 2, 2, 3, 3, 4, 4, 1, 1, 2, 2, 3, 3, 4,
4, 1, 1, 2, 2, 3, 3, 4, 4, 1, 1, 2, 2, 3, 3, 4,
1, 1, 2, 2, 3, 3, 4, 1, 1, 2, 2, 3, 3, 4, 1, 1,
2, 2, 3, 3, 1, 1, 2, 2, 3, 3, 1, 1, 2, 2, 3, 3)
freq <- c(c(5.9, 5.9, 5.9, 5.9, 5.9, 5.9, 5.9, 5.9, 8, 8, 8, 8, 8, 8, 8,
8, 10, 10, 10, 10, 10, 10, 10, 10, 5.9, 5.9, 5.9, 5.9, 5.9, 5.9,
5.9, 8, 8, 8, 8, 8, 8, 8, 10, 10, 10, 10, 10, 10, 10, 5.9, 5.9,
5.9, 5.9, 5.9, 5.9, 8, 8, 8, 8, 8, 8, 10, 10, 10, 10, 10, 10))
ID <-c(1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 7,
8, 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 8,
1, 2, 3, 4, 5, 6, 8, 1, 2, 3, 4, 5, 6, 8, 1, 2,
3, 4, 5, 6, 1, 2, 3, 4, 5, 6, 1, 2, 3, 4, 5, 6)
week <- as.factor(week)
Methode <- as.factor(Methode)
ID <- as.factor(ID)
M1 <- lmer(log(tijd) ~ Methode + (1|freq) + (1|week) + (1|ID), data = data)
summary(M1)
anova(M1, type = "III", test = "F")
plot(M1)
plot(cooks.distance(M1))
summary(glht(M1, linfct = mcp(Methode = "Tukey")))
这是我的数据集的第一行:https ://imgur.com/RYr9nNp 我目前正在写我的学士论文,如果你能帮助我,那就太好了。
提前致谢 :)
解决方案
在您的 lmer 代码中,您依赖data
但您没有转换data$Methode
为因子,那么您的代码应该是:
data$week <- as.factor(data$week)
data$Methode <- as.factor(data$Methode)
data$ID <- as.factor(data$ID)
M1 <- lmer(log(tijd) ~ Methode + (1|freq) + (1|week) + (1|ID), data = data)
您在图像中显示的数据不包含与 lmer 代码相同的列名,因此您必须注意(例如:frequentie 和 freq)
推荐阅读
- php - 如何从 php 的 .gz 存档中获取原始文件名
- python - 使用 Pandas 的计算 apply & lambda
- c++ - C++ Boost Asio 池线程与 lambda 函数并通过引用变量传递
- python - 如何在 Python 制表中合并/连接表格?
- c++ - 移动 shared_ptr
进入 shared_ptr > - sql - 过滤掉不是这个日期格式的记录oracle
- java - 用 java8 中引入的 Map 方法替换 containsKey 检查的有效方法
- reactjs - 在三元运算符语句中值始终为空
- gradle - 在应用程序重启时调用自己的 API(Grails 3.3.9)(Tomcat 重启)
- sql - CASE WHEN foo.bar IS NULL THEN 0 ELSE 1 END 范式是否有更简洁的版本?