[R] coplot: extract regression output

[Marco Helbich]

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