output |
---|
github_document |
logitnorm
package provides support for the univariate
logit-normal
distribution. In
addition to the usual random, density, percentile, and quantile function, it
helps with estimating distribution parameters from observations statistics.
# From CRAN
install.packages("logitnorm")
# Or the the development version from GitHub:
# install.packages("devtools")
devtools::install_github("bgctw/logitnorm")
See the package vignette for an introduction.
A simple example estimates distribution parameters from observation statistics of mode 0.7 and upper quantile 0.9. Next, the density is computed and plotted across a range of quantiles.
(theta <- twCoefLogitnormMLE(0.7,0.9))
#> Error in twCoefLogitnormMLE(0.7, 0.9): could not find function "twCoefLogitnormMLE"
x <- seq(0,1, length.out=81)
d <- dlogitnorm(x, mu=theta[1,"mu"], sigma=theta[1,"sigma"])
#> Error in dlogitnorm(x, mu = theta[1, "mu"], sigma = theta[1, "sigma"]): could not find function "dlogitnorm"
plot(d~x,type="l")
#> Error in eval(predvars, data, env): object 'd' not found
abline(v=c(0.7,0.9), col="grey")
#> Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...): plot.new has not been called yet