乗算代入データセットでパッケージのcoxph()関数を使用し、結果をプールしようとすると警告が発生しました。survival警告メッセージには次のように記載されています。"In mice.df(m, lambda, dfcom, method) : Large sample assumed.
mice()再現可能な例を以下に示します (公開されているデータを使用しcoxph()、これらのデータを使用することの適切性についてあまり心配する必要はありません)。
library(mice)
library(survival)
#load publically available data
data(pbc)
#select variables for the reproducable example
pbc.select <- pbc[pbc$status %in% c(0,1) , c("id", "time", "status", "trt")]
imp <- mice(pbc.select) #impute trt
fit <- with(imp, coxph(Surv(time, status) ~ trt)) #fit coxph in each imp
pool(fit) #pool the models; get's the error
この警告は、pool()関数が dfcom を要求しようとしているために発生したようです。
dfcom <- df.residual(object)
wheredf.residual()は、このコンテキストで参照されるオブジェクトには適用されません。coxph
class(fit) # "mira" "matrix"
class(fit$analyses[[1]]) "coxph"
私の質問は
- 目的に合わせて正しい構文を使用していますか?
- もしそうなら
pool()、適切な情報を提供する方法はありますか? - この仮定は結果にどのように影響しますか?