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.