我想使用 r 包 mlogit 运行多项混合模型分析。 不幸的是我的数据结构有问题，因为 r 似乎找不到数据（下面的错误）。 Chid 是 Choiceid，CASE 是主题 ID。我还想添加一个组变量，因为我有一个主题因素组，这意味着参与者被分成 3 个不同的组，选择场景略有不同。有没有人知道如何将其添加到模型中。一般来说，我仍然是选择建模 r 的初学者。我想在我的模型中考虑个人级别和组级别。
> mixdata <- dfidx(df_model2, choice = "Choice", varying = 4:19, sep = "_",idx = list(c("chid", "CASE")), idnames = c("chid", "alt"))
> head(mixdata, 10)
~~~ indexes ~~~~
indexes: 1, 1, 2
>
> #radomisierte effekte bestimmen
> rdm <- c(ArtHausarzt = "n", ArtBeratung = "n", ArtEhrenamt = "n", ArtSelbsthilfe = "n", SettingGruppe = "n")
> #modellschätzung
> ml_mix <-mlogit(formula = Choice ~Art + Setting, data = mixdata, rpar = rdm, panel = TRUE, method = "nr")
Error in array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), :
'data' must be of a vector type, was 'NULL'type here
我的数据具有以下结构：
| 案例 | 团体 | 选择 | 艺术_1 | 艺术_2 | 格式_1 | 格式 2 | ... | 孩子 |
|---|---|---|---|---|---|---|---|---|
| 1 | 2 | 1 | 心理治疗 | 自私 | Einzeln | 集团 | ... | 1 |
| 1 | 2 | 2 | Beratung | 豪萨尔茨 | 集团 | 集团 | ... | 2 |
我尝试使用小插图作为模板：https://cran.r-project.org/web/packages/mlogit/vignettes/c5.mxl.html varying 中的列是正确的，因此变量的所有名称都是正确的。