library(agricolae)
x<-c(1,3,5,2,3,1)
y<-c(2,5,4,3,2,3)
z<-c(4,10,NA,6,7,NA)
# original data
cbind(x,y,z)
modelo<-lm(z~x+y)
# worked data
zz<-wxyz(modelo,x,y,z)
zz
# Use function gxyz
xyz<-gxyz(zz)
par(mfrow=c(2,2),cex=0.6)
# startgraph
# original data
plot(x,y,main="original data")
#endgraph
# startgraph
# worked data
plot(xyz$x,xyz$y,xlab="x",ylab="y",main="worked data
with model")
#endgraph
Run the code above in your browser using DataLab