Displaying 1 result from an estimated 1 matches for "comptab".
Did you mean:
compta
2011 Feb 28
0
Gamma mixture models with flexmix
...model=FLXMRglm(family='Gamma'))
mdl <- getModel(fmo, which=which)
print(smry <- summary(mdl))
print(prm <- parameters(mdl))
plotGammaMixture(dat, prm['shape', ],
prm['shape', ]*prm['coef.(Intercept)', ],
smry at comptab[, 'prior'])
}
##
## Works well for a single Gamma distribution.
##
set.seed(78483)
dat1 <- rgamma(6000, shape=2, rate=0.5)
modelGammas(dat1)
set.seed(78483)
dat2 <- rgamma(6000, shape=5, rate=0.3)
modelGammas(dat2)
##
## Please help me get it to work for mixtures of
## two Gamma...