# NOT RUN {
# start by estimating lognormal parameters from moments
.mean <- 1
.var <- c(1.3,2)^2
parms <- getParmsLognormForMoments(.mean, .var)
#
# computed moments must equal previous ones
(ans <- getLognormMoments(parms[,"mu"], parms[,"sigma"]))
cbind(.var, ans[,"var"])
#
getLognormMedian(mu = log(1), sigma = log(2))
getLognormMode(mu = log(1), sigma = c(log(1.2),log(2)))
# }
Run the code above in your browser using DataLab