问题描述
我正在尝试为stan_glmer.nb
(rstanarm
)输出创建表,但是来自软件包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 (将#修改为@)