在R中对mblogit模型使用函数ranef时出现问题

siv3szwd  于 2023-04-03  发布在  Git
关注(0)|答案(1)|浏览(112)

我目前正在研究R中的logit模型,我正在使用mclogit库中的mblogit模型。我想使用lme4库中的函数ranef来研究随机效应。
但是,我得到了以下错误
UseMethod(“ranef”)中的错误:没有适用于“c('mmblogit ',' mblogit','mmclogit',' mclogit','lm')”类对象的“ranef”的方法
我想知道是否有人可以向我指出一个函数,它可以在mblogit对象上执行类似于ranef的操作。
下面的代码是这个问题的一个玩具示例。程序以相同的错误结束。

library(mclogit)
library(tidyverse)
options(warn=-1)

mi_df <- data.frame(
  "factor" = c("a", "b", "c", "c", "d", "a"),
  "category" = c("1", "1","1", "2", "2", "2"),
  "province" = c("aa", "bb", "aa", "bb", "aa", "bb"),
  "number1" = c(1.2, 3.4, 4.5, 5.6, 4.5, 3.6),
  "number2" = c(6, 4, 9, 4, 2, 6)
  )

mi_df$category <- relevel(as.factor(mi_df$category), ref="1")

modelo_multinomial0 <- mblogit(category ~  factor + number1+number2 , random = ~ 1  | (province), 
                               data=mi_df, maxit = 30, epsilon = 1e-04)
summary(modelo_multinomial0) 
confint(modelo_multinomial0) 

library(lme4)
ranef(modelo_multinomial0)
axzmvihb

axzmvihb1#

通过深入研究str(modelo_multinomial0),看起来$random.effects * 可能 * 是您想要的:

modelo_multinomial0$random.effects[1,1][[1]]
            [,1]
[1,] 5.418201e-41
[2,] 6.020224e-42

我将用一些更大的示例进行实验,以查看此组件的维度和值是否与您的期望相匹配。

相关问题