Dieter Menne:
Thanks for the suggestion and link, it looks good. I think my trouble now is
lack of basic R knowledge.

Cheers,
Michael

For example,

I tried:
histogram(~bbContag | sc_recov, data=dat)

This would work fine if I could select cases from sc_recov. How can I select
more than one case? I know I can do:

sc_recov.21 <- dat[dat$sc_recov=21,]

but how could I select all cases where sc_recov = 21, 31 or 41?


On Thu, Oct 2, 2008 at 3:33 AM, Michael Just <mgjust@gmail.com> wrote:

> Hello,
> If I use:
>
> p <- ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram()
>
> and then:
>
> p + facet_grid(. ~ sc_recov)
>
> Its a little crazy because I have 48 different values in 'sc_recov'.
> Instead I want to select cases from 'sc_recov' and only use three at a time:
>
> I tried:
> > p + facet_grid(sc_recov=="21" ~.)
> Error in check_formula(formula, varnames) :
>   Formula contains variables not in list of known variables
>
> How can I select cases?
>
> Any ideas or suggestions?
> Thanks,
> M Just
>
> On Thu, Oct 2, 2008 at 2:42 AM, ONKELINX, Thierry <
> Thierry.ONKELINX@inbo.be> wrote:
>
>> Dear Michael,
>>
>> Try ggplot2. Use something like
>>
>> install.packages("ggplot2")
>> library(ggplot2)
>> recov <- 0:2
>> n <- 1000
>> all <- data.frame(bbED = rnorm(3 * n, mean = recov), recov =
>> factor(rep(recov, n)))
>> ggplot(data = all, aes(x = bbED)) + geom_histogram() + facet_grid(. ~
>> recov)
>> ggplot(data = all, aes(x = bbED)) + geom_histogram() + facet_grid(recov
>> ~ .)
>> ggplot(data = all, aes(x = bbED, colour = recov)) + geom_density()
>>
>> You'll find more information on ggplot2 at http://had.co.nz/ggplot2/
>>
>> HTH,
>>
>> Thierry
>>
>> ------------------------------------------------------------------------
>> ----
>> ir. Thierry Onkelinx
>> Instituut voor natuur- en bosonderzoek / Research Institute for Nature
>> and Forest
>> Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
>> methodology and quality assurance
>> Gaverstraat 4
>> 9500 Geraardsbergen
>> Belgium
>> tel. + 32 54/436 185
>> Thierry.Onkelinx@inbo.be
>> www.inbo.be
>>
>> To call in the statistician after the experiment is done may be no more
>> than asking him to perform a post-mortem examination: he may be able to
>> say what the experiment died of.
>> ~ Sir Ronald Aylmer Fisher
>>
>> The plural of anecdote is not data.
>> ~ Roger Brinner
>>
>> The combination of some data and an aching desire for an answer does not
>> ensure that a reasonable answer can be extracted from a given body of
>> data.
>> ~ John Tukey
>>
>> -----Oorspronkelijk bericht-----
>> Van: r-help-bounces@r-project.org [mailto:r-help-bounces@r-project.org]
>> Namens Michael Just
>> Verzonden: donderdag 2 oktober 2008 7:29
>> Aan: r-help@r-project.org
>> Onderwerp: [R] Multiple hist(ograms) - One plot
>>
>> Hello,
>> I am trying to plot multiple histograms with the same scales, etc into
>> one
>> plot. The commands below produce a 3 page PDF with each histogram
>> occupying
>> the upper right quadrant. And use slightly different scales on the X and
>> Y
>> axes.
>>
>> > s21 <- dat[dat$sc_recov=="21",]
>> > s21.ED <- subset(s21, select=(bbED))
>> > s31 <- all[all$sc_recov=="31",]
>> > s31.ED <- subset(s31, select=(bbED))
>> > s41 <- all[all$sc_recov=="41",]
>> > s41.ED <- subset(s41, select=(bbED))
>> > pdf("234_1_ED.pdf")
>> > par(mfrow = c(2, 2), oma = c(0, 0, 2, 0))
>> > hist(s21.ED) + hist(s31.ED) + hist(s41.ED)
>> > dev.off()
>>
>> Also, for example if I add something like col="x", or main="x" to the
>> hist(x, ...) nothing appears on the plot.  ?
>>
>> Please advise,
>> Thank you kindly,
>> M Just
>>
>>         [[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.
>>
>> Dit bericht en eventuele bijlagen geven enkel de visie van de schrijver
>> weer
>> en binden het INBO onder geen enkel beding, zolang dit bericht niet
>> bevestigd is
>> door een geldig ondertekend document. The views expressed in  this message
>> and any annex are purely those of the writer and may not be regarded as
>> stating
>> an official position of INBO, as long as the message is not confirmed by a
>> duly
>> signed document.
>>
>
>

	[[alternative HTML version deleted]]

