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]]