Re: [R] how to use 'points' function to plot two curves with errbar

From: Yogesh Tiwari <yogesh.mpi_at_googlemail.com>
Date: Tue, 15 Jun 2010 13:50:18 +0530

Dear John,
Great thanks,
I would appreciate if any expert can throw more light on it.

Thanks in advance,

Best Regards,
Yogesh

On 6/14/10, John Kane <jrkrideau_at_yahoo.ca> wrote:
> Hi Yogesh,
>
> I think you accidentelly replied only to me and not to the mailing list.
> I'd suggest that you copy your post to me to the list to get input from some
> of the
> real experts.
>
> Currently I don't see any provision for differently coloured error bars in
> the the help page but I easilly could be missing something.
>
> Also you have a fairly busy plot there, you might want to consider having
> the cuves in separate panels. Perhaps something like this: you can probably
> do this much better in lattice or ggplot but can give you an idea.
> #=========================================================
> x=1:7
> y=6:12
>
> uperror=y+2
> downerror=y-2
>
> opar <- par(mfrow=c(2,2))
> errbar(x,y, uperror, downerror,type="l")
> errbar(x,y, uperror, downerror, type="l")
> errbar(x,y, uperror, downerror, type="l")
> errbar(x,y, uperror, downerror, type="l")
> par(opar)
> #==========================================================
> Also we really should have some sample data., It have to be the entire data
> set or even the real data as
> long as it has the same format and layout.
>
> Have a look at ?dput for way to include some easily useable data.
>
> Sorrry I'm not more help.
>
> --- On Sun, 6/13/10, Yogesh Tiwari <yogesh.mpi_at_googlemail.com> wrote:
>
> From: Yogesh Tiwari <yogesh.mpi_at_googlemail.com>
> Subject: Re: [R] how to use 'points' function to plot two curves with errbar
> To: "John Kane" <jrkrideau_at_yahoo.ca>
> Received: Sunday, June 13, 2010, 11:21 PM
>
> Dear John,
> Thanks,
> I am using errbar in the package Hmisc,
>
> par(new=TRUE), worked between two errbar ploting commands but the new
> problem is now 'color of arrow (bars)', its always black for all the four
> curves whereas I have specified separate color for each curve. I want to
> have similar color of arrow as curve. My command which plots four curves
> and respective errbar subsequently the generated plot is attached herewith.
> Kindly help.
>
>
>> x.val <- as.integer(names(co2mean))
>> errbar(x.val, co2mean, co2mean + co2sd, co2mean - co2sd, xaxt='n',
>> col="black", xlab=NA,ylab=NA)# obs error bar
>> lines(x.val, co2mean, col="black", lwd=2)
>
>
>> par(new=TRUE)
>
>> errbar(x.val, co2tm3.month.mean, co2tm3.month.mean + co2sd.tm3,
>> co2tm3.month.mean - co2sd.tm3, xaxt='n', yaxt='n', col="red",
>> xlab=NA,ylab=NA)# tm3 error bar
>
>> lines(x.val, co2tm3.month.mean, col="red", lwd=2)
>
>> par(new=TRUE)
>
>> errbar(x.val, co2lsce.month.mean, co2lsce.month.mean + co2sd.lsce,
>> co2lsce.month.mean - co2sd.lsce, xaxt='n',yaxt='n', col="green",
>> xlab=NA,ylab=NA)# lsce error bar
>
>> lines(x.val, co2lsce.month.mean, col="green", lwd=2)
>
>> par(new=TRUE)
>
>> errbar(x.val, co2patra1.month.mean, co2patra1.month.mean + co2sd.actm,
>> co2patra1.month.mean - co2sd.actm, xaxt='n',yaxt='n', col="blue",
>> xlab=NA,ylab=NA)# actm error bar
>
>> lines(x.val, co2patra1.month.mean, col="blue", lwd=2)
>
>
>
> Best regards,
> Yogesh
>
>
>

> On Sat, Jun 5, 2010 at 11:36 PM, John Kane <jrkrideau_at_yahoo.ca> wrote:
>
> The first thing we need to know is what "errbar" are you using?
>
> There are at least two, on in the Hmisc package and one in sfsmisc.
>
>
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html and provide commented, minimal,
> self-contained,reproducible code.
>
>
> I'm not sure what you are trying to do here but perhaps looking at "new" in
> ?par might hrlp you.
>
> Example:
> plot(1:10, col="red")
> par(new=TRUE)
> plot(10:1, col="blue")
>
>
>
>
>
> --- On Sat, 6/5/10, Yogesh Tiwari <yogesh.mpi_at_googlemail.com> wrote:
>
>> From: Yogesh Tiwari <yogesh.mpi_at_googlemail.com>
>
>> Subject: [R] how to use 'points' function to plot two curves with errbar
>> To: r-help_at_stat.math.ethz.ch
>> Received: Saturday, June 5, 2010, 10:06 AM
>
>
>
>
>> Dear R Users,
>> I am using R on windows.
>>
>> how to use 'points' function to plot two curves with
>> errbar
>>
>> I am doing like:
>>
>> x.val <- as.integer(names(co2mean))
>
>> errbar(x.val, co2mean, co2mean + co2sd, co2mean - co2sd,
>> xaxt='n', col=1,
>> xlab=NA,ylab=NA)# obs error bar
>> lines(x.val, co2mean, col=1, lwd=2)
>>
>> errbar(x.val, co2tm3.month.mean, co2tm3.month.mean +
>
>> co2sd.tm3,
>> co2tm3.month.mean - co2sd.tm3, xaxt='n', col=2,
>> xlab=NA,ylab=NA)# model
>> error bar
>> lines(x.val, co2tm3.month.mean, col=2, lwd=2)
>>
>> Above code replaces first error bar and its line with
>
>> second. Kindly help,
>> how to use 'points' function here so both above can get
>> over ploted
>>
>> Thanks,
>> Regards,
>> Yogesh
>
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>
>> and provide commented, minimal, self-contained,
>> reproducible code.
>>
>
>
>
>
>
> --
> Yogesh K. Tiwari (Dr.rer.nat),
> Scientist,
> Indian Institute of Tropical Meteorology,
>
> Homi Bhabha Road,
> Pashan,
> Pune-411008
> INDIA
>
> Phone: 0091-99 2273 9513 (Cell)
> : 0091-20-25904350 (O)
> Fax : 0091-20-258 93 825
>
>
>
>

-- 
Yogesh K. Tiwari (Dr.rer.nat),
Scientist,
Indian Institute of Tropical Meteorology,
Homi Bhabha Road,
Pashan,
Pune-411008
INDIA

Phone: 0091-99 2273 9513 (Cell)
         : 0091-20-25904350 (O)
Fax    : 0091-20-258 93 825

______________________________________________
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 Tue 15 Jun 2010 - 13:08:46 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 Tue 15 Jun 2010 - 14:00:33 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