829d82d3363d3d7c726bd1157d9c58a244c3ff47
[morpheus.git] / patch_Bettina / code.R
1 library("flexmix")
2 data("tribolium", package = "flexmix")
3 set.seed(1234)
4 ## uses FLXMRglm() from package.
5 tribMix <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species,
6 k = 2, nrep = 5, data = tribolium,
7 model = FLXMRglm(family = "binomial"))
8 parameters(tribMix); logLik(tribMix)
9 source("FLXMRglm.R")
10 set.seed(1234)
11 ## uses FLXMRglm() from source file which allows to specify link.
12 tribMixNew <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species,
13 k = 2, nrep = 5, data = tribolium,
14 model = FLXMRglm(family = "binomial"))
15 parameters(tribMixNew); logLik(tribMixNew)
16 set.seed(1234)
17 tribMixlogit <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species,
18 k = 2, nrep = 5, data = tribolium,
19 model = FLXMRglm(family = binomial(link = logit)))
20 parameters(tribMixlogit); logLik(tribMixlogit)
21 set.seed(1234)
22 tribMixprobit <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species,
23 k = 2, nrep = 5, data = tribolium,
24 model = FLXMRglm(family = binomial(link = probit)))
25 parameters(tribMixprobit); logLik(tribMixprobit)