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]]

