微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

为什么`performance :: model_parameters`对负二项式`rstanarm`模型抛出错误

如何解决为什么`performance :: model_parameters`对负二项式`rstanarm`模型抛出错误

我正在尝试为stan_glmer.nbrstanarm输出创建表,但是来自软件包model_parameters的{​​{1}}抛出一个奇怪的错误,我不确定如何解决。也许这是一个错误

可复制的示例:

performance

在这里我将为您保留输出,因为它并不重要,但是该功能有效。 现在是负二项式模型:

library(rstanarm)
library(parameters)

dat<-data.frame(x<-rnorm(500),z=rep(c("A","B","C","D","E"),100),y=.2+x*.7)

mod1<-stan_glmer(y~x+(x|z),data=dat)

model_parameters(mod1,effects="all")

现在出现错误消息:

dat.nb<-data.frame(x=rnorm(500),y=rnbinom(500,size=1,prob = .5))

mod2<-stan_glmer.nb(y~x+(x|z),data=dat.nb)

model_parameters(mod2,effects="all")

函数似乎是从错误的地方拉出的,或者我丢失了什么?据我在帮助文档中所看到的,没有任何东西表明需要指定一个负二项式模型。此外,该功能可用于Error in `$<-.data.frame`(`*tmP*`,"parameter",value = c("(Intercept)",: replacement has 3 rows,data has 1 模型:

lme4

此模拟数据存在模型收敛性等问题,但是library(lme4) mod1<-lmer(y~x+(x|z),effects="all") mod2<-glmer.nb(y~x+(x|z),effects="all") 适用于model_parameters模型,但不适用于glmer.nb模型。知道这里发生了什么吗?

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。