FYI this now being logged upstream.
--- Begin Message ---
#16: errbar() does not work correctly on a log axis
--------------------------+-------------------------------------------------
Reporter: dupontct | Owner: dupontct
Type: defect | Status: new
Priority: major | Milestone: 3.5-1 release
Component: Hmisc | Version: 3.4-2
Keywords: Hmisc errbar |
--------------------------+-------------------------------------------------
The errbar() function does not work correctly on a log axis, due
to the code assuming the axis is linear.
Here is a tiny testcase to demonstrate:
{{{
library(Hmisc)
x <- c(0.03,0.1,0.3,1,3,10,30)
y <- c(10,6,2,0.8,0.3,0.06,0.021)
sd <- c(5,2,1,0.3,0.1,0.2,0.08)
par(mfrow=c(2,1))
errbar(x,y,y-sd,y+sd)
errbar(x,y,y-sd,y+sd, log="x")
}}}
It should be possible for the code to check if par xlog=T and then
adjust the x±cap to cope with the log scale.
--
Ticket URL: <http://biostat.mc.vanderbilt.edu/trac/Hmisc/ticket/16>
Hmisc <http://example.org/>
My example project
--- End Message ---
--
Hell, there are no rules here - we're trying to accomplish something.
-- Thomas A. Edison
