Great! Thank you, Brian. To answer your question about intercept and slopes, I got them from a covariance analysis that I had already conducted. It seems like I can not use the regressions command for the model that I used to get the intercepts and slopes. I guess 2 factors are the maximum.
+ ddply(test, c("country","treatment"), + function(x) { + coef(lm(total~ treatment+ year+ country + treatment:year, x)) + }) Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") : contrasts can be applied only to factors with 2 or more levels Any way to get around this? How about getting back to the abline command? Thank you. -- View this message in context: http://r.789695.n4.nabble.com/grey-colored-lines-and-overwriting-labels-i-qqplot2-tp3657119p3666473.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.