[R] Histograms by two factors
Peter Dalgaard
p.dalgaard at biostat.ku.dk
Fri Feb 6 11:34:53 CET 2004
"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
More information about the R-help
mailing list