Estou tentando ajustar um modelo de efeitos mistos considerando o pacote gamlss
.
O modelo que eu estou ajustando é definido a seguir:
library(gamlss)
ModelGamlss <- gamlss(log(Var1)~re(fixed=~log(Var2)+log(Var3)+
(Var4)+(Var5),random = list(Var6=pdDiag(~Var7))),
data=Data2)
No entanto, o seguinte erro é apresentado:
Error in gamlss(log(Var1) ~ re(fixed = ~log(Var1)~re(fixed=~log(Var2)+log(Var3)+(Var4)+ :
The data contains NA's, use data = na.omit(mydata)
Após me deparar com este erro, realizei a seguinte modificação na sintaxe do modelo considerado:
ModelGamlss <- gamlss(log(Var1)~re(fixed=~log(Var2)+log(Var3)+
(Var4)+(Var5),random = list(Var6=pdDiag(~Var7))),
data=na.omit(Data2))
Todavia, um novo erro é apresentado:
Error in lm.wfit(x, z, w, method = "qr", ...) : 0 (non-NA) cases
Como poderia solucionar os problemas enfrentados anteriormente? Visto que o primeiro erro informa que existem NA's na base de dados, porém, ao observar a base de dados nenhum NA é observado, e mais, qual seria o motivo dessa segunda mensagem de erro estar aparecendo ao tentar ajustar o modelo acima?