我目前正在研究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)
1条答案
按热度按时间axzmvihb1#
通过深入研究
str(modelo_multinomial0)
,看起来$random.effects
* 可能 * 是您想要的:我将用一些更大的示例进行实验,以查看此组件的维度和值是否与您的期望相匹配。