random-betablocker.R
In catdata: Categorical Data

## ----echo=FALSE,eval=FALSE----------------------------------------------------
#  options(width=80)

## ----eval=FALSE---------------------------------------------------------------
#  library(flexmix)
#  data(betablocker)

## ----eval=FALSE---------------------------------------------------------------
#  betablocker$Treatment <- as.factor(betablocker$Treatment)

## ----eval=FALSE---------------------------------------------------------------
#  GlmT <- glm(cbind(Deaths, Total	- Deaths) ~ Treatment, family = "binomial",
#  data = betablocker)
#  summary(GlmT)

## ----eval=FALSE---------------------------------------------------------------
#  GlmTC <- glm(cbind(Deaths, Total - Deaths) ~ Treatment + as.factor(Center),
#               family 	= "binomial", data = betablocker)
#  summary(GlmTC)

## ----eval=FALSE---------------------------------------------------------------
#  library(glmmML)

## ----eval=FALSE---------------------------------------------------------------
#  MixedGH4 <- glmmML(cbind(Deaths, Total - Deaths) ~ Treatment, cluster=Center,
#                     method = c("ghq"), n.points = 4, boot = 0, data=betablocker)
#  
#  summary(MixedGH4)

## ----eval=FALSE---------------------------------------------------------------
#  MixedGH20 <- glmmML(cbind(Deaths, Total - Deaths) ~ Treatment, cluster=Center,
#                      method = c("ghq"), n.points = 20, boot = 0, data=betablocker)
#  
#  summary(MixedGH20)

## ----eval=FALSE---------------------------------------------------------------
#  set.seed(5)

## ----results='hide',eval=FALSE------------------------------------------------
#  detach(package:glmmML)

## ----eval=FALSE---------------------------------------------------------------
#  MixFix3 <-stepFlexmix(cbind(Deaths, Total - Deaths) ~ 1 | Center,	model =
#    FLXMRglmfix(family = "binomial", fixed = ~ Treatment), k = 3, nrep = 5,
#                        data = betablocker)

## ----eval=FALSE---------------------------------------------------------------
#  MixFix3

## ----eval=FALSE---------------------------------------------------------------
#  parameters(MixFix3)

## ----results='hide',eval=FALSE------------------------------------------------
#  library(flexmix)

## ----eval=FALSE---------------------------------------------------------------
#  summary(MixFix3)
#  summary(refit(MixFix3))

## ----eval=FALSE---------------------------------------------------------------
#  set.seed(5)

## ----eval=FALSE---------------------------------------------------------------
#  MixFix4 <-stepFlexmix(cbind(Deaths, Total - Deaths) ~ 1 | Center, model =
#    FLXMRglmfix(family = "binomial", fixed = ~ Treatment), k = 4, nrep = 5,
#                        data = betablocker)

## ----eval=FALSE---------------------------------------------------------------
#  MixFix4

## ----eval=FALSE---------------------------------------------------------------
#  parameters(MixFix4)

## ----eval=FALSE---------------------------------------------------------------
#  summary(MixFix4)
#  summary(refit(MixFix4))

## ----echo=FALSE,results='hide',eval=FALSE-------------------------------------
#  detach(package:flexmix)