1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
module SModel.ASRV (module SModel.ASRV,
module SModel.MixtureModel)
where
import SModel.MixtureModel
import Probability.Distribution.Gamma (gamma)
import Probability.Distribution.Transform (logNormal)
gammaRatesDist alpha = gamma alpha (1/alpha)
gammaRates alpha n base = rateMixture base $ uniformDiscretize (gammaRatesDist alpha) n
logNormalRatesDist sigmaOverMu = logNormal lmu lsigma where x = log(1+sigmaOverMu^2)
lmu = -0.5*x
lsigma = sqrt x
logNormalRates sigmaOverMu n base = rateMixture base $ uniformDiscretize (logNormalRatesDist sigmaOverMu) n
freeRates rateDist base = rateMixture base rateDist
|
