Displaying 1 result from an estimated 1 matches for "popc100kpc".
2011 Oct 21
4
plotting average effects.
...pc100[which(dat3$popc100 < min(dat$popc100))] <- min(dat$popc100)
> pred1 <- predict(mod, type="response")
> pred2 <- predict(mod, newdata=dat2, type="response")
> pred3 <- predict(mod, newdata=dat3, type="response")
> pop.group <- cut(dat$popc100kpc, breaks=quantile(dat$popc100kpc,
> seq(0,1,by=.3)), include.lowest=T)
> means <- by(cbind(pred1, pred2, pred3), list(pop.group), apply, 2, mean)
> means <- do.call(rbind, means)
> par(mar=c(7,4,4,2))
> plot(c(1,10), range(c(means)), type="n", xlab="",
+...