Joel Fürstenberg-Hägg
2009-Dec-10 12:53 UTC
[R] Get rid of automatic title from Anova plots
Hi all, I'm having a problem when putting the four plots (Residuals vs Fitted, Normal Q-Q, Scale - Location, and Residuals vs Leverage) from an Anova using plot(aov()) into a pdf. A string "aov(df[,i] ~ cat))" is added as the main title, so when I use mtext() the two titles are overlapping. Does anyone know how to get rid of that title? I've tried plot(aov(df[,i] ~ cat), main="") but without succes. pdf("Anova plots 0809.pdf", height=10, width=10) par(mfrow=c(2,2), oma=c(2,2,4,2)) # Anova test for LT50 with different anthocyanin scores df=fieldTrial[idx0809, c(31:32)] cat=fieldTrial[idx0809, c(14)] for (i in colnames(df)) { print(i) print(summary(aov(df[,i] ~ cat))) plot(aov(df[,i] ~ cat)) mtext(text=paste("Anova test for", i, "with different anthocyanin scores"), cex=1.5, side=3, outer=TRUE) } dev.off() Best regards, Joel _________________________________________________________________ Hitta kärleken i vinter! http://dejting.se.msn.com/channel/index.aspx?trackingid=1002952 [[alternative HTML version deleted]]