Подмножество объекта mids перед использованием функции pool.scalar
Мои данные продольные, и у меня есть 2 подгруппы (в зависимости от статуса исхода). Я выполнил многократное вменение, используя пакет MICE, чтобы решить мои проблемы с отсутствующими данными.
Впоследствии я использовал pool.scalar функция для оценки объединенной медианы и (Q1-Q3) (см. кодирование ниже).
Тем не менее, я хотел бы сделать это для каждого timepoint а также subgroup отдельно. Когда я использую subset_datlist функция для подмножества моего объекта импутации сред, тогда я теряю свою информацию вменения, поэтому я не могу использовать pool.scalar функция потом (которая требует объект среднего).
Используемое кодирование для объединенных оценок Медиана (q1-q3)
new <- complete(imputedobs,'long') m<-(imputedobs$m)
CRPnew<-new$CRP
Q1 <- rep(NA, m)
U <- rep(NA, m)
Q3 <- rep(NA, m)
Q2 <- rep(NA, m)
for (i in 1:m) {
Q1[i] <- quantile(complete(imputedobs, i)$CRP, c(0.25))
Q2[i] <- median(complete(imputedobs, i)$CRP)
Q3[i] <- quantile(complete(imputedobs, i)$CRP, c(0.75))
U[i] <- var(complete(imputedobs, i)$CRP) / nrow(nhanes)
}
pool.scalar(Q1, U, method = "rubin") # Rubin 1987
pool.scalar(Q2, U, method = "rubin")
pool.scalar(Q3, U, method = "rubin")