Dear list:
Before posting this message, I read the posting guide, examined the manuals,
searched the R help files, and examined in detail the books that I have
available on R.
I am using R version 2.8.0 and WinXP.
I want to pause the R Graphics window to permit the use of the File,
History, and Resize options between each graphic in a series of graphics. If
I use the functions
par(ask=T)
or
devAskNewPage(ask=T)
between the graphics there is a pause but one cannot use the File, History,
and Resize options. For an example use
demo(Graphics)
What I want to do is approximated by the following code:
trat <- factor(c(rep("A",5),rep("B",5),rep("C",5)))
trat
resp <-
c(rnorm(5,mean=5,sd=1.75),rnorm(5,mean=6.5,sd=1.75),rnorm(5,mean=15,sd=1.75)
)
aovOBJ <- aov(resp~trat)
summary(aovOBJ)
AOV <- function(aovOBJ){
print(tmp <- summary.lm(aovOBJ))
boxplot(tmp$resid,col="grey")
reply <- readline(prompt="Next:")
hist(tmp$resid,freq=F,col="grey",breaks="fd",xlab="",ylab="",main="")
xfit <- seq(min(tmp$resid),max(tmp$resid),length=40)
yfit <- dnorm(xfit,mean=mean(tmp$resid),sd=sd(tmp$resid))
lines(xfit,yfit,col="red")
reply <- readline(prompt="Next:")
require(car,quietly=T)
qq.plot(tmp$resid,ylab="empirical",main="",col="black")
detach("package:car")
reply <- readline(prompt="Next:")
plot(density(tmp$resid),main="")
reply <- readline(prompt="Next:")
plot.design(aovOBJ$model)
reply <- readline(prompt="Next:")
tmp <- TukeyHSD(aovOBJ)
plot(tmp)
}
AOV(aovOBJ)
There ought to exist a lot easier user friendly way to do this. Does anyone
know one?
Thank you very much,
pbm
Peter B. Mandeville
cel: 444 860 3204
tel: 52 444 826 2346-49 ext 532
fax: 52 444 826 2350
P.D. Favor de confirmar la llegada de este correo. Gracias.
[[alternative HTML version deleted]]
______________________________________________
[email protected] 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.
