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


