Я использую HMI на двухуровневых данных (учащиеся на курсах) с отсутствующими данными на уровне учащихся. Код выдает следующую ошибку (Ошибка: использование стека C 7969776 слишком близко к пределу), когда я включаю двоичное значение (пол) с отсутствующими данными.
Ниже приведен воспроизводимый пример.
library(MCMCglmm)
library(hmi)
library(mice)
df <- data.frame(post = rtnorm(100,60,20,0,100),
pre = rtnorm(100,40,20,0,100),
gender = rbinom(n=100, size=1, prob=0.20),
course = rep( c("A","B","C","D"), 100*c(0.1,0.2,0.65,0.05) ))
df<span class="math-container">$post[1:round((0.3)*length(df$</span>post),digits = 0)] <- NA
df<span class="math-container">$gender[round((0.2)*length(df$</span>post),digits = 0):round((0.5)*length(df$post),
digits = 0)] <- NA
hmi_test <- hmi(df, model_formula = post ~ 1 + pre + gender + (1|course ),
M = 2, maxit = 5, list_of_types = NULL, nitt = 3000, burnin = 1000)
list_of_types_maker(df)