Estou tentando realizar o gráfico de resíduos do modelo de efeitos mistos por meio da função ggplot2. No entanto, após realizar uma busca encontrei algumas funções disponíveis mas o que me parece é que para a função nlme elas não estão funcionando. Os gráficos que eu pretendo realizar são os do exemplo abaixo: ![](https://i.imgur.com/3Wg3HEd.png[/img]) As rotinas computacionais que eu tentei inicialmente se encontram abaixo, vejam os erros que estão aparecendo ao executar a função em ggplot2. ```` library(splines) library(ggplot2) library(nlme) library(gridExtra) setwd("C:\\Users\\Desktop") datanew1 = read.table("dadosnew.csv", header = T, sep=";", dec = ",") datanew1$DummyVariable = as.factor(datanew1$DummyVariable) datanew1$Variable2 = as.factor(datanew1$Variable2) datanew1$Variable3 = as.factor(datanew1$Variable3) ############################################################################# ############################## Model ######################################## ############################################################################# model <- lme(Response~(bs(Variable1, df=3)) + DummyVariable, random=~1|Variable2/Variable3, datanew1, method="REML") completemodel <- update(model, weights = varIdent(form=~1|DummyVariable)) p1 <- qplot(.fitted, .resid, data = completemodel) + geom_hline(yintercept = 0) + geom_smooth(se = FALSE) Erro: `data` must be a data frame, or other object coercible by `fortify()`, not an S3 object with class lme Run `rlang::last_error()` to see where the error occurred. p2 <- qplot(sample =.stdresid, data = completemodel, stat = "qq") + geom_abline() grid.arrange(p1,p2) Erro: `data` must be a data frame, or other object coercible by `fortify()`, not an S3 object with class lme Run `rlang::last_error()` to see where the error occurred. Além disso: Warning message: `stat` is deprecated ```` Uma outra maneira na qual eu tentei realizar o gráfico foi com a função abaixo, mas não obtive sucesso. ```` ggplot(completemodel, aes(.fitted, .resid)) + geom_point() Erro: `data` must be a data frame, or other object coercible by `fortify()`, not an S3 object with class lme Run `rlang::last_error()` to see where the error occurred. ````