> From: Abhijit Dasgupta <adasgu...@araastat.com> > Date: October 31, 2010 1:30:02 AM EDT > To: Matt Shotwell <shotw...@musc.edu> > Subject: Re: [R] ForestPlot or similar > > I just did something very similar using ggplot's pointrange geom. In the > following, I'm plotting hazard ratios, for which the nominal value is 1 and > not 0. x has 5 columns: drug, hr, hr.lcb, hr.ucb, and group, and I'm faceting > by group. If you want the plots horizontal, add coord_flip() to the command > --- as it stands the plots are vertically oriented. > > require(ggplot2) > p <- ggplot(x, aes(x=drug, y = hr, ymin=hr.lcb, ymax=hr.ucb))+ > geom_pointrange()+ facet_grid(.~group) > p <- p + xlab('Drug') + ylab('Hazard ratio')+ > geom_hline(y=1, col='red', lty=2) > > > Abhijit > On Oct 30, 2010, at 5:31 PM, Matt Shotwell wrote: > >> Here is a small function for forest plots in R, with an example: >> >> http://biostatmatt.com/wiki/r-credplot >> >> -Matt >> >> On Sat, 2010-10-30 at 11:40 -0400, Mestat wrote: >>> Here is one example: >>> I have three vectors (mean,lower interval, upper interval) >>> mean<-c(2,4,6,8) >>> l<-c(1,2,3,4) >>> u<-c(4,8,12,16) >>> How would I plot that if I want to use the FORESTPLOT function. I dont need >>> to use the TABLETEXT option. >>> I am working in something like this: >>> tabletext<-c(NA,NA,NA,NA,NA) >>> mean<-c(NA,2,4,6,8) >>> l<-c(NA,1,2,3,4) >>> u<-c(NA,4,8,12,16) >>> forestplot(tabletext,mean,l,u,zero=0) >>> But I am having a problem with the length of the dimension... >>> Thanks in advance, >>> Marcio >>> >> >> -- >> Matthew S. Shotwell >> Graduate Student >> Division of Biostatistics and Epidemiology >> Medical University of South Carolina >> >> ______________________________________________ >> 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.