Hi Basil,
I think you actually want par(new=TRUE) rather than plot.new()
y1<-c(min(Y), max(Y))
x1<-c(min(X), max(X))
pp <- plot(effect('X', m1), xlim = x1, ylim = y1)
par(new=TRUE)
plot(Y~X, xlim = x1, ylim = y1, pch = 16, col = "black")
but this then has another problem - effects alters the par
Dear R users,
I am analyzing data using a gls model from the nlme package
The model is:
m2<-gls(C~LDH + SI3, method = "ML", weights = vpDist) ## I will later fit
this model using REML
C = total soil carbon, LDH = log of a distance that I am interested in, and
SI3 = is a index of invasion severi