Rev Language Reference


dnBimodalLognormal - Bimodal Lognormal Distribution

A bimodal lognormal distribution, that is, with probability p a value is distributed according to the first lognormal distribution and with probability 1-p from the second lognormal distribution.

Usage

dnBimodalLognormal(Real mean1, Real mean2, RealPos sd1, RealPos sd2, Probability p)

Arguments

mean1 : Real (pass by const reference)
The mean (in log-space) of the first lognormal distribution.
mean2 : Real (pass by const reference)
The mean (in log-space) of the second lognormal distribution.
sd1 : RealPos (pass by const reference)
The standard deviation of the first lognormal distribution.
sd2 : RealPos (pass by const reference)
The standard deviation of the secind lognormal distribution.
p : Probability (pass by const reference)
The probability to belong to the first distribution.

Domain Type

Example

p ~ dnBeta(1.0,1.0)
x ~ dnBimodalLognormal(mean1=-1,mean2=1,sd1=0.1,sd2=0.1,p=p)
x.clamp( exp(1) )
moves[1] = mvSlide(p, delta=0.1, weight=1.0)
monitors[1] = screenmonitor(printgen=1000, separator = "        ", x)
mymodel = model(p)
mymcmc = mcmc(mymodel, monitors, moves)
mymcmc.burnin(generations=20000,tuningInterval=100)
mymcmc.run(generations=200000)

See Also