Deepayan,
Thanks for the advice.

-Micahel

On Fri, Oct 3, 2008 at 1:24 AM, Deepayan Sarkar
<[EMAIL PROTECTED]>wrote:

> On 10/2/08, Michael Just <[EMAIL PROTECTED]> wrote:
> > Dieter and Thierry:
> >
> >  Per you suggestions I have tried:
> >
> >  ggplot2 from Thierry:
> >
> > > p <- ggplot(dat, aes(x=bbContag, y=..density..)) + geom_histogram()
> >
> > > p + facet_grid(. ~ sc_recov %in% c(21,31,41))
> >
> >  But get the followinng error:
> >
> > Error in check_formula(formula, varnames) :
> >   Formula contains variables not in list of known variables
> >
> >
> > If I don't make selections it works, but makes that very busy plot.
> >
> >  histogram from Dieter:
> >  > histogram(~bbContag | sc_recov %in% c(21,31,41), data=dat)
> >  This produces a plot with only two histograms side by side, but I have
> >  chosen three groups, where is the third? I tried only choosing two
> groups
> >  and it plotted two histograms different than when I tried choosing
> three.
>
> Actually, Dieter's suggestion was to use the subset argument, which
> would go something like
>
> histogram(~bbContag | sc_recov, data=dat,
>          subset = (sc_recov %in% c(21,31,41)))
>
> Given that your 'sc_recov' seems to be a numeric variable, you will
> probably get better annotation with
>
> histogram(~bbContag | factor(sc_recov), data=dat,
>          subset = (sc_recov %in% c(21,31,41)))
>
> -Deepayan
>
> >  Please advise,
> >  Thank you for your continued help,
> >  Michael Just
> >
> >
> >
> >  On Thu, Oct 2, 2008 at 5:21 AM, ONKELINX, Thierry
> >
> > <[EMAIL PROTECTED]>wrote:
> >
> >  > Michael,
> >  >
> >  > Use %in% to select multiple cases:
> >  >
> >  > dat[dat$sc_recov %in% c(21, 31, 41), ]
> >  >
> >  >
> >  >
> ------------------------------------------------------------------------
> >  > ----
> >  > 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
> >  > [EMAIL PROTECTED]
> >  > 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: [EMAIL PROTECTED] [mailto:
> [EMAIL PROTECTED]
> >  > Namens Michael Just
> >  > Verzonden: donderdag 2 oktober 2008 10:59
> >  > Aan: r-help@r-project.org
> >  > CC: [EMAIL PROTECTED]
> >  > Onderwerp: Re: [R] Multiple hist(ograms) - One plot
> >  >
> >  > 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 <[EMAIL PROTECTED]>
> 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 <
> >  > > [EMAIL PROTECTED]> 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
> >  > >> [EMAIL PROTECTED]
> >  > >> 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: [EMAIL PROTECTED]
> >  > [mailto:[EMAIL PROTECTED]
> >  > >> 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]]
> >  >
> >  > ______________________________________________
> >  > 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]]
> >
> >  ______________________________________________
> >  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.
> >
>

        [[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.

Reply via email to