The answer to my post is yes (which I just figured out).
Solution:
#super small version of R code for pd plot using apply
a <- rbind(c(0:1111)*(max(m$x1)-min(m$x1))/1111 +
min(m$x1),c(0:1111)*0-99999)
b <- matrix(rep(c(0:1111)*(max(m$x1)-min(m$x1))/1111 + min(m$x1), nrow(m)),
nrow(m), 1112, byrow=T)
a[2,] <- apply(b,2,FUN=function(x)
{mean(predict(lm1,cbind(m[,-match("x1",names(m))],x1=x))) })
plot(a[1,],a[2,],xlab="x1",ylab="Response",type="l",main="Partial Dependence
Plot")
Mike Dugas
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
