Hi folks, Can someone please help me in identifying the type of plot shown here? Sample R code or specific package name would be of help as well. Thanks, -Girish http://www.nabble.com/file/p24477714/example_plot.jpg -- View this message in context: http://www.nabble.com/Help-needed-in-identifying-type-of-plot-tp24477714p24477714.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]]
On 14 Jul 2009, at 12:30, Girish A.R. wrote:> Can someone please help me in identifying the type of plot shown here?Its a Forest Plot Chris. ----------------------- Gynaecological Cancer Research Laboratories, UCL EGA Institute for Women's Health, University College London, Paul O'Gorman Building, 72 Huntley Street, London WC1E 6DD United Kingdom Telephone; 020 3108 2007 Fax; 020 3108 2010 [[alternative HTML version deleted]]
It is a plot with pointranges. Here is an example with the ggplot2
package
library(ggplot2)
#use the diamond data set from ggplot2
diamonds$A <- diamonds$depth < 60
dmod <- lm(price ~ cut * A, data=diamonds)
cuts <- unique(diamonds[, c("cut", "A")])
cuts <- cbind(cuts, predict(dmod, cuts,
se=T)[c("fit","se.fit")])
#the actual plot
ggplot(cuts, aes(x = cut, y = fit, ymin = fit - se.fit, ymax=fit +
se.fit, colour = factor(A))) + geom_pointrange(position position_dodge(width =
0.2)) + coord_flip()
More details on the ggplot2 website: 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 at 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 at r-project.org [mailto:r-help-bounces at r-project.org]
Namens Girish A.R.
Verzonden: dinsdag 14 juli 2009 13:31
Aan: r-help at r-project.org
Onderwerp: [R] Help needed in identifying type of plot
Hi folks,
Can someone please help me in identifying the type of plot shown here?
Sample R code or specific package name would be of help as well.
Thanks,
-Girish
http://www.nabble.com/file/p24477714/example_plot.jpg
--
View this message in context:
http://www.nabble.com/Help-needed-in-identifying-type-of-plot-tp24477714
p24477714.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]
______________________________________________
R-help at 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.
Thanks, Chris and Thierry! I think I can explore both your solutions.
best,
-Girish
It is a plot with pointranges. Here is an example with the ggplot2
package
library(ggplot2)
#use the diamond data set from ggplot2
diamonds$A <- diamonds$depth < 60
dmod <- lm(price ~ cut * A, data=diamonds)
cuts <- unique(diamonds[, c("cut", "A")])
cuts <- cbind(cuts, predict(dmod, cuts,
se=T)[c("fit","se.fit")])
#the actual plot
ggplot(cuts, aes(x = cut, y = fit, ymin = fit - se.fit, ymax=fit +
se.fit, colour = factor(A))) + geom_pointrange(position position_dodge(width =
0.2)) + coord_flip()
More details on the ggplot2 website: 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 at 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
--
View this message in context:
http://www.nabble.com/Help-needed-in-identifying-type-of-plot-tp24477714p24478359.html
Sent from the R help mailing list archive at Nabble.com.
It looks like a forest plot normally used in meta-analysis work. See http://www.bmj.com/cgi/content/full/322/7300/1479 . See http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=114 for an example from rmeta package with accompanying code. I think that package metafor also gives forest plots. --- On Tue, 7/14/09, Girish A.R. <garamach at gmail.com> wrote:> From: Girish A.R. <garamach at gmail.com> > Subject: [R] Help needed in identifying type of plot > To: r-help at r-project.org > Received: Tuesday, July 14, 2009, 7:30 AM > > Hi folks, > > Can someone please help me in identifying the type of plot > shown here? > Sample R code or specific package name would be of help as > well. > > Thanks, > -Girish > > http://www.nabble.com/file/p24477714/example_plot.jpg > > > -- > View this message in context: http://www.nabble.com/Help-needed-in-identifying-type-of-plot-tp24477714p24477714.html > Sent from the R help mailing list archive at Nabble.com. > > ??? [[alternative HTML version deleted]] > > ______________________________________________ > R-help at 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. >__________________________________________________________________ Looking for the perfect gift? Give the gift of Flickr! http://www.flickr.com/gift/