![](/img/trans.png)
[英]How to get confidence intervals after extracting robust standard errors in R?
[英]how to pool MI confidence intervals of robust mixed model in r?
我可以運行rlmer
model 和由mice
產生的 object ,但是當我嘗試合並結果時,會收到消息Error: No tidy method for objects of class rlmerMod
。 有替代方案嗎?
下面是我的數據和模型的可重現示例:
set.seed(1)
library(data.table)
library(robustlmm)
library(mice)
library(miceadds)
dt <- data.table(id = rep(1:10, each=3),
group = rep(1:2, each=15),
time = rep(1:3, 10),
sex = rep(sample(c("F","M"),10,replace=T), each=3),
x = rnorm(30),
y = rnorm(30))
setDT(dt)[id %in% sample(1:10,4) & time == 2, `:=` (x = NA, y = NA)][
id %in% sample(1:10,4) & time == 3, `:=` (x = NA, y = NA)]
# Multiple imputation -------------------------------------------------------------------
pm <- make.predictorMatrix(dt)
pm[,c('x','y')] <- 0
pm[c('x','y'), 'id'] <- -2
imp <- mice(dt, pred = pm, meth = "2l.pmm", seed = 1, m = 2, print = FALSE, maxit = 20)
# Modelling -----------------------------------------------------------------------------
m <- with(imp, rlmer(y ~ 1 + time * group + sex + (1 | id), REML=F))
pool.fit <- pool(m)
> pool.fit <- pool(m)
Error: No tidy method for objects of class rlmerMod
In addition: Warning message:
In get.dfcom(object, dfcom) : Infinite sample size assumed. # I don't get this warning using my real data
謝謝!
編輯:
正如@BenBolker 所評論的那樣, library(broom.mixed)
讓pool.fit
可以正常運行。 然而, summary(pool.fit,conf.int = TRUE)
返回估計值,但NaN
表示自由度、p 值和置信區間。
library(broom.mixed)
pool.fit <- pool(m)
summary(pool.fit,conf.int = TRUE)
term estimate std.error statistic df p.value 2.5 % 97.5 %
1 (Intercept) -1.31638288 1.2221584 -1.07709683 NaN NaN NaN NaN
2 time 0.02819273 0.4734632 0.05954578 NaN NaN NaN NaN
3 group 1.49581955 0.8776475 1.70435124 NaN NaN NaN NaN
4 sexM -0.61383469 0.7137998 -0.85995356 NaN NaN NaN NaN
5 time:group -0.25690287 0.3005254 -0.85484573 NaN NaN NaN NaN
我不知道是否需要另一個參數(例如,用於定義 df 方法)。
現在,我嘗試tbl_regression(m)
但它也不起作用:
> tbl_regression(m)
pool_and_tidy_mice(): Tidying mice model with
`mice::pool(x) %>% mice::tidy(exponentiate = FALSE, conf.int = TRUE, conf.level = 0.95)`
Error in match.call() : ... used in a situation where it does not exist # how to correct this?
In addition: Warning message:
In get.dfcom(object, dfcom) : Infinite sample size assumed. # again, this warning don't occur with my original data
任何提示?
只需加載broom.mixed
package,它有rlmerMod
對象的 tidiers。 ( broom.mixed
get_methods()
function:
remotes::install_github("bbolker/broom.mixed")
library(broom.mixed)
print(get_methods(), n = Inf)
# A tibble: 22 × 4
class tidy glance augment
<chr> <lgl> <lgl> <lgl>
1 allFit TRUE TRUE FALSE
2 brmsfit TRUE TRUE TRUE
3 gamlss TRUE TRUE FALSE
4 gamm4 TRUE TRUE TRUE
5 glmmadmb TRUE TRUE TRUE
6 glmmTMB TRUE TRUE TRUE
7 gls TRUE TRUE TRUE
8 lme TRUE TRUE TRUE
9 lmList4 TRUE FALSE FALSE
10 mcmc TRUE FALSE FALSE
11 mcmc.list TRUE FALSE FALSE
12 MCMCglmm TRUE FALSE FALSE
13 merMod TRUE TRUE TRUE
14 MixMod TRUE FALSE FALSE
15 ranef.mer FALSE FALSE TRUE
16 rjags TRUE FALSE FALSE
17 rlmerMod TRUE FALSE FALSE
18 stanfit TRUE FALSE FALSE
19 stanreg TRUE TRUE FALSE
20 TMB TRUE FALSE FALSE
21 varComb TRUE FALSE FALSE
22 varFunc TRUE FALSE FALSE
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.