I'm having problems with the identify function.
I am trying to identify the points on the residuals graph of the adjusted model, however the identify function
is giving an error.
library(mgcv)
require(gamm4)
fit4.gamm <- gamm(log(massaseca)~factor(Trat)+s(Tempo,k=10,bs="ps",m=2,
by=factor(Trat)),
random=list(id=pdSymm(~Tempo)),data=dados)
comp_lme = fit4.gamm$lme
x11()
plot(comp_lme)
Error:
identify(comp_lme)
Error in xy.coords(x, y, setLab = FALSE) :
'x' is a list, but does not have components 'x' and 'y'
Take a reproducible example from the documentation:
What you are trying:
Instead, make your own base plot:
This will work.