David,
I would do this with the HH package.
Rich
install.packages("HH") ## if you don't already have it
library(HH)
## the first call is to learn the number .8162 and see what the scaling is.
normal.and.t.dist(mu.H0=5.01, obs.mean=log(300), std.dev=.77,
Use.alpha.left=TRUE, Use.alpha.right=FALSE,
Use.obs.mean=TRUE)
ticks <- c(outer(c(1, 2, 5), c(10,100,1000)))
axis(side=1, at=log(ticks), labels=ticks, line=7)
## the second call makes the display much prettier, including putting in the
log axis.
normal.and.t.dist(mu.H0=5.01, obs.mean=log(300), std.dev=.77,
Use.alpha.left=TRUE, Use.alpha.right=FALSE,
alpha.left=.8162)
ticks <- c(outer(c(1, 2, 5), c(10,100,1000)))
axis(side=1, at=log(ticks), labels=ticks, line=7)
axis(side=1, at=log(300), labels="", line=7, col="blue", tck=-.06)
axis(side=1, at=log(300), labels=300, line=4.5, tck=FALSE, col.axis="blue",
cex.axis=1.3)
axis(side=1, at=1.5, labels=expression(ufc/m^3), xpd=NA, line=5.4,
tck=FALSE, tick=FALSE)
## the third call suppresses the xbar axis
normal.and.t.dist(mu.H0=0, obs.mean=(log(300)-5.01)/.77,
Use.alpha.left=TRUE, Use.alpha.right=FALSE,
alpha.left=.8162)
ticks <- c(outer(c(1, 2, 5), c(10,100,1000)))
axis(side=1, at=(log(ticks)-5.01)/.77, labels=ticks, line=7)
axis(side=1, at=(log(300)-5.01)/.77, labels="", line=7, col="blue",
tck=-.06)
axis(side=1, at=(log(300)-5.01)/.77, labels=300, line=4.5, tck=FALSE,
col.axis="blue", cex.axis=1.3)
axis(side=1, at=(1.5-5.01)/.77, labels=expression(ufc/m^3), xpd=NA,
line=5.4, tck=FALSE, tick=FALSE)
[[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.
