Displaying 1 result from an estimated 1 matches for "modelgammas".
2011 Feb 28
0
Gamma mixture models with flexmix
...#39;, freq=FALSE, ylim=ylim,
main=paste(KK, 'component(s) in model\nmodel prior(s) =',
paste(round(pis, 2), collapse=', ')))
cols <- rainbow(KK)
for (ii in 1:KK) {
lines(x, y[[ii]], col=cols[ii])
}
}
##
## Model dat as a mixture of Gammas then plot.
##
modelGammas <- function(dat, which='BIC') {
set.seed(939458)
fmo <- stepFlexmix(dat ~ 1, k=1:3,
model=FLXMRglm(family='Gamma'))
mdl <- getModel(fmo, which=which)
print(smry <- summary(mdl))
print(prm <- parameters(mdl))
plotGammaMixture(dat, prm[&...