[morpheus.git] / patch_Bettina / code.R

library("flexmix")
data("tribolium", package = "flexmix")
set.seed(1234)
## uses FLXMRglm() from package.
tribMix <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species, 
                       k = 2, nrep = 5, data = tribolium,
                       model = FLXMRglm(family = "binomial"))
parameters(tribMix); logLik(tribMix)
source("FLXMRglm.R")
set.seed(1234)
## uses FLXMRglm() from source file which allows to specify link.
tribMixNew <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species, 
                          k = 2, nrep = 5, data = tribolium,
                          model = FLXMRglm(family = "binomial"))
parameters(tribMixNew); logLik(tribMixNew)
set.seed(1234)
tribMixlogit <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species, 
                            k = 2, nrep = 5, data = tribolium,
                            model = FLXMRglm(family = binomial(link = logit)))
parameters(tribMixlogit); logLik(tribMixlogit)
set.seed(1234)
tribMixprobit <- initFlexmix(cbind(Remaining, Total - Remaining) ~ Species, 
                             k = 2, nrep = 5, data = tribolium,
                             model = FLXMRglm(family = binomial(link = probit)))
parameters(tribMixprobit); logLik(tribMixprobit)
Commit	Line	Data
cbd88fe5 BA	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)