Re: [R] plot two columns of a matrix with standard deviation

From: Henrique Dallazuanna <wwwhsd_at_gmail.com>
Date: Thu, 21 Feb 2008 17:41:45 -0300

Try this:

plot(x[,1], ylim=range(x[c(1,3)]), xaxt="n", xlab="", pch=16, col="blue") points(x[,3], pch=16, col="red")
axis(1, at=1:nrow(x), labels=rownames(x), las=2) segments(1:nrow(x), x[,1]-x[,2], 1:nrow(x), x[,1]+x[,2], col="blue") segments(1:nrow(x), x[,3]-x[,4], 1:nrow(x), x[,3]+x[,4], col="red")

On 21/02/2008, Paul Hammer <paul.hammer_at_tfh-wildau.de> wrote:
>
> Henrique Dallazuanna schrieb:
> Perhaps:
>
> plot(x[c(1,3)], pch=16)
> segments(x[,1],x[,3]-x[,2], x[,1], x[,3]+x[,4], col="red")

>
>
> On 21/02/2008, Paul Hammer <paul.hammer_at_tfh-wildau.de> wrote:
>
>
> hi members,
>
> i try to plot two columns of a matrix as points with standard deviation.
> a legend should also be there. the one points should be green and the
> others red...
>
> column 1 and 3 are the values to plot and column 2 and 4 are the
> attendant standard deviation (sd)... as xlab should show the rownames. i
> tried a lot of attempts but without success:(...
>
> R code:
> > PDE1A_int_sd_avg[1:4]
>
> prostate_avg prostate_sd breast_avg breast_sd
> 2590617 7.057866 0.65549017 6.922251 0.1240702
> 2590616 4.407165 0.51983406 4.127747 0.4471652
> 2590615 4.806831 0.58847360 4.309884 0.1884463
> 2590618 3.314213 1.02590469 3.310029 0.6075215
> 2590612 7.878071 0.34838454 6.961996 0.4294136
> 2590613 7.148809 0.63653777 6.013015 0.6484706
> 2590611 5.382736 0.51058425 4.669508 0.5124079
> 2590607 4.448152 1.54271678 3.377348 0.6634425
> 2590608 4.184283 0.44209393 2.767696 0.3458292
> 2590606 6.237610 0.15831587 4.900560 0.3556708
> 2590601 2.871831 0.13629742 2.427491 0.7291913
> 2590600 2.777127 1.50248150 3.181344 0.2703518
> 2590599 3.854438 0.12178512 3.560036 0.4389107
> 2590628 4.065403 0.67835058 3.468860 0.8218848
> 2590619 1.623058 0.07643176 1.946055 0.5335189
> 2590586 1.102444 0.42467133 1.214922 0.5288023
> 2590591 3.842315 0.42584177 2.792450 0.2539151
> 2590645 4.882622 0.75542588 5.073177 0.2141262
>
> help would be great...
>
> thanks
> paul
>
> @ jim thanks for the hint assign :)
>
> ______________________________________________
> 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.
>
>
>
>
> thanks for first but it works not really right. y axis should be the values
> for columns and the x axis just the rowname (2590617, 2590616 etc). it
> belongs always two point together (one for column prostate and one for
> column breast). they should be positioned above. additional the standard
> devation should be shown for every point...
>
> for a better understanding see attached file...
>
> thanks
> paul
>
>

-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O

______________________________________________
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.
Received on Thu 21 Feb 2008 - 20:47:08 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 21 Feb 2008 - 21:30:16 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive