Dear list,
I would like to extract the model details (coefficients, p-values etc.)
for the regression lines within a coplot. How are these results
accessible? A code example is as follows:
dat <- data.frame(a=rnorm(111), b=rnorm(111), c=rnorm(111))
coplot(dat$a ~ dat$b | dat$c, data=dat,
cex=1, number=3, columns=3, col="black", pch=16, overlap = 0.1,
panel=function(x,y,...) {
points(x,y,pch=16)
panel.smooth(x,y,span=.8,iter=5,...)
abline(lm(y ~ x), col="blue")} )
Thank you.
Marco