我想使用序数逻辑回归(我的响应变量是序数),该回归逻辑回归适用于2个随机变量,并且用于具有交互作用的定量预测变量(我的公式是:ordinal_variable~ quantitative_variable:habitat + (1|community) + (1|species)
。
我正在使用clmm
分析数据(请参见下面的脚本)并获得了预期的结果,但是我注意到clmm
设计用于在响应和预测变量为阶乘时使用。 / p>
然后我尝试了polmer
(请参见下面的脚本),但是没有得到任何答案。
有人会对如何分析这些数据有任何建议吗?
library(ordinal)
model1 <- clmm(as.factor(ordinal_variable)~
quantitative_variable:habitat + (1|community) + (1|species),
data=baseline)
summary(model1)
library(MPDiR)
library(lme4)
model2 <- polmer(as.factor(ordinal_variable) ~
quantitative_variable:habitat + (community - 1 | Obs) +
(species - 1 | Obs), data=baseline)
summary(model2)