[R] Histograms by two factors

[Peter Dalgaard]

"Briggs, Meredith M" <Meredith.Briggs at team.telstra.com> writes:

> Hi
> 
> 
> I am trying to print out means, STDs and histograms under two sets
> of factors. I can manage it for one set - see below but not for two
> sets. That is, I want ot print out the mean STD and Histogram for
> each ITEM code within each DELIVERABLE code. 

Tapply accepts a list of factors, tapply(y,list(Df1,Df2),mean), and so
forth (and the help page does say so...).

> In addition I can only get to view the histogram for the last item.

You do get to view the others, just not for very long ;-)

> How do you get R to stop overriding the histogram for eg level 1 for
> factor 1 by level 2 of factor 1?

par(ask=TRUE) is one possible answer. 
Others include 
- turning on plot history (in the Windows GUI) or use
  recordPlot/replayPlot
- switching graphics device to a file and viewing it later
- opening a separate graphics device for each histogram
- setting up for multiple plots on the same page.


> 
> thanks
> 
> 
> 
> 
> X11()
> 
> Indat <- read.table ("C:/testdata.txt",header=T)
> 
> B <- c(Indat[,1],Indat[,2],Indat[,3],Indat[,4],Indat[,5])
> y <- Indat[,5]/Indat[,4]
> DELIV <- Indat[,1]
> ITEM <- Indat[,3]
> Df1 <- factor(DELIV)
> Df2 <- factor(ITEM)
> 
> d1 <- tapply(y,Df1, mean)
> d2 <- tapply(y, Df1,sd)
> d3 <- tapply(y, Df1,max)
> d4 <- tapply(y, Df1,hist)
> 
> print(d1)
> print(d4)
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
> 

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907