x <- c(rep(0,20), exp(rnorm(1000,0.05)))
fx <- f.log(x)
## Not run:
# oldpar <- par(mfrow = c(2, 3))
# plot(x,main="exp(normal)+zeros")
# qqnorm(x)
# T3plot(x)
# plot(fx,main="optimized offset")
# qqnorm(fx)
# T3plot(fx)
# par(oldpar)
# ## End(Not run)
Run the code above in your browser using DataLab