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.
