Параллельное вычисление множественного вменения с использованием пакета мышей R
Я хочу выполнить 150 множественных вменений с помощьюmice
вR
, Однако, чтобы сэкономить некоторое вычислительное время, я хотел бы разделить процесс на параллельные потоки (как предложил Стеф ван Буурен в «Гибком вменении для отсутствующих данных»).
Мой вопрос: как это сделать?
Я могу представить 2 варианта:
opt.1:
imp1<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
imp2<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
imp...<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
imp150<-mice(data, m=1, pred=quicktry, maxit=15, seed=1)
а затем объединить вменения вместе с помощьюcomplete
а такжеas.mids
впоследствии
opt.2:
imp1<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
imp2<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
imp...<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
imp150<-mice(data, m=1, pred=quicktry, maxit=15, seed=VAL_1to150)
добавляяVAL_1to150
в противном случае мне кажется (я могу ошибаться), что если все они будут работать с одним и тем же набором данных и одним и тем же начальным числом, у вас будет 150 раз одинаковый результат.
Есть ли другие варианты?
Спасибо