modelsummary()适配rstanarm负二项分布stan_glm模型失败问题求助
modelsummary()适配rstanarm负二项分布stan_glm模型失败问题求助
大家好,我最近把modelsummary更新到2.0.0版本后遇到了一个棘手的问题:用rstanarm包的stan_glm()拟合泊松模型时,modelsummary()能正常生成结果表格,但换成负二项分布模型时就完全失效了。有没有大佬能给我一些排查方向或者解决建议呀?
我用到的包如下:
library(tinytable) library(modelsummary) library(rstanarm) library(tidyverse) library(broom.mixed)
泊松模型可正常运行
先测试了泊松模型,整个流程完全没问题:
# Poisson works poisson <- stan_glm( cyl ~ mpg, data = mtcars, family = poisson(link = "log"), seed = 853 ) modelsummary(poisson)
而且tidy()的输出也符合预期:
tidy(poisson)
负二项模型无法生成表格
但换成负二项分布模型后,modelsummary()就没法生成结果表格了:
# Negative binomial doesn't work neg_binomial <- stan_glm( cyl ~ mpg, data = mtcars, family = neg_binomial_2(link = "log"), seed = 853 ) # modelsummary doesn't built a table modelsummary(neg_binomial)
奇怪的是,tidy()明明能返回带有term列的数据框:
tidy(neg_binomial)
报错信息
运行modelsummary(neg_binomial)时弹出的错误:
modelsummary could not extract the required information from a model object.
备注:内容来源于stack exchange,提问作者Rohan Alexander
相关产品推荐
相关产品推荐

