Thierry, Thanks. This worked. Cheers, Michael
On Fri, Oct 3, 2008 at 2:46 AM, ONKELINX, Thierry <[EMAIL PROTECTED]>wrote: > Michael, > > You get this error because you make the subset at the wrong place. Try > > p <- ggplot(dat[dat$sc_recov %in% c(21,31,41), ], aes(x=bbContag, > y=..density..)) + geom_histogram() > p + facet_grid(. ~ sc_recov) > > or > > subdat <- dat[dat$sc_recov %in% c(21,31,41), ] > p <- ggplot(subdat, aes(x=bbContag, y=..density..)) + geom_histogram() > p + facet_grid(. ~ sc_recov) > 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 > > > ------------------------------ > *Van:* Michael Just [mailto:[EMAIL PROTECTED] > *Verzonden:* vrijdag 3 oktober 2008 6:52 > *Aan:* ONKELINX, Thierry > *CC:* r-help@r-project.org; [EMAIL PROTECTED] > > *Onderwerp:* Re: [R] Multiple hist(ograms) - One plot > > 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. > > 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. >> > > 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.
