如何运行R中具有混合效应的多项式logistic回归?

wpx232ag  于 2023-01-06  发布在  其他
关注(0)|答案(1)|浏览(147)

我正在尝试运行一个混合效应的多项式逻辑回归,假设我有以下变量:
参与者(10名参与者,每人有10个观察结果)词(10个不同的词,每位参与者看一次)IV(一些两水平分组变量)DV(可以是1、2或3)
我该如何用ppt和word作为随机变量进行多项式逻辑回归?
以下是一些示例数据:

ppt <- rep(1:10, each = 10)
word <- rep(1:10, times = 10)
IV <- rep(1:2, times = 50)
DV <- sample(x = c(1,2,3), size = 100, replace = TRUE)

d <- as.data.frame(cbind(ppt, word, IV, DV))
mccptt67

mccptt671#

我认为brms是最适合这类模型的软件包,但是这个软件包适合贝叶斯模型,所以你必须指定先验。我已经拟合了一个示例模型,但没有指定先验,所以它使用了你不应该使用的默认值。

library(brms)

d$DV <- as.factor(d$DV)
d$IV <- as.factor(d$IV)

model_fit <- brm(DV ~ IV + (1|ppt) + (1|word),
                 data = d,
                 family = categorical())

相关问题