Ошибка в X% *% fixef (объект): несовместимые аргументы

Question

Я получаю эту ошибку, когда пытаюсь предсказать nlme. И я не знаю, как ее решить, мой код выглядит следующим образом:

mod1.lme <- lme(pRNFL~ status +  Time:status, random= ~1|DBID/Eye_ID, data = long, na.action = na.omit)
summary(mod1.lme)


plot(fitted(mod1.lme), resid(mod1.lme), col = long$status) 
dim(long)
length(fitted(long))
length(fitted(mod1.lme))

long2 <- long[!is.na(long$pRNFL),]  
dim(long2)
plot(long2$pRNFL, fitted(mod1.lme))
plot(long2$pRNFL, fitted(mod1.lme), pch = 16, cex = 0.8, xlab = "Measured pRNFL", ylab = "Model fitted pRNFL")
abline(a = 0, b = 1, col = "red", lty = 2)

A_AQP4 <- subset(long, status=='A-AQP4')
F_AQP4 <- subset(long, status=='F-AQP4')
N_AQP4 <- subset(long, status=='N-AQP4')
A_MOG <- subset(long, status=='A-MOG')
F_MOG <- subset(long, status=='F-MOG')
N_MOG <- subset(long, status=='N-MOG')


x <- c(0:800)/100
A_AQP4data <- data.frame(Time = x, status = "A-AQP4")
F_AQP4data <- data.frame(Time = x, status = "F-AQP4")
N_AQP4data <- data.frame(Time = x, status = "N-AQP4")
A_MOGdata <- data.frame(Time = x, status = "A-MOG")
F_MOGdata <- data.frame(Time = x, status = "F-MOG")
N_MOGdata <- data.frame(Time = x, status = "N-MOG")

yA_AQP4 <- predict(mod1.lme, newdata = A_AQP4data, level = 0)
yF_AQP4 <- predict(mod1.lme, newdata = F_AQP4data, level = 0)
yN_AQP4 <- predict(mod1.lme, newdata = N_AQP4data, level = 0)
yA_MOG <- predict(mod1.lme, newdata = A_MOGdata, level = 0)
yF_MOG <- predict(mod1.lme, newdata = F_MOGdata, level = 0)
yN_MOG <- predict(mod1.lme, newdata = N_MOGdata, level = 0)

Любая идея, как я могу решить это, пожалуйста? Ошибка следующая:

> yA_AQP4 <- predict(mod1.lme, newdata = A_AQP4data, level = 0)
Error in X %*% fixef(object) : non-conformable arguments

Это сводка mod1.lme (называемый здесь mod2) на случай, если это будет полезно [] [2

Уровни моих переменных разные

Что я делаю не так?

Спасибо !!