Я пытаюсь включить случайный наклон в модель glmmTMB, но получаю следующую ошибку:
Error in KhatriRao(sm, t(mm)) : (p <- ncol(X)) == ncol(Y) is not TRUE
Я прочитал здесь , что это ошибка может происходить из значений NA в нефакторных переменных. Если я удаляю значения NA в моей нефакторной выходной переменной, я получаю ту же ошибку. У меня нет значений NA в моем предсказателе, что является фактором. Я также прочитал здесь , что эта ошибка может быть связана с тем, что формула модели указана неправильно, но я думал, что моя формула модели верна (возможно, я ошибаюсь).
Доступен доступ к моему набору данных здесь и код ниже.
Любая помощь для решения этой ошибки будет принята с благодарностью.
test <- read.csv("test_1.csv", header = TRUE)
glmmTMB(lweight_t ~ (scale_id | 1), data = test, family = gaussian)
test <- filter(long, lweight_t != 'NA')
glmmTMB(lweight_t ~ (scale_id | 1), data = test, family = gaussian)