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