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

## ----results='hide',eval=FALSE------------------------------------------------
#  library(catdata)
#  
#  data(birth)
#  attach(birth)

## ----eval=FALSE---------------------------------------------------------------
#  table1 <- table(Sex, Membranes, Cesarean, Induced)
#  
#  ftable(table1)

## ----eval=FALSE---------------------------------------------------------------
#  m4 <- loglin(table1, margin=list(c(1,2,3,4)), fit=TRUE)
#  cat("deviance(m4)=", m4$lrt, "df(m4)=", m4$df, "\n")
#  
#  m3 <- loglin(table1, margin=list(c(1,2,3), c(1,2,4), c(1,3,4), c(2,3,4)), fit=TRUE)
#  cat("deviance(m3)=", m3$lrt, "df(m3)=", m3$df, "\n")
#  
#  m2 <- loglin(table1, margin=list(c(1,2), c(1,3), c(1,4), c(2,3), c(2,4), c(3,4)),
#               fit=TRUE)
#  cat("deviance(m2)=", m2$lrt, "df(m2)=", m2$df, "\n")
#  
#  m1 <- loglin(table1, margin=list(c(1), c(2), c(3), c(4)), fit=TRUE)
#  cat("deviance(m1)=", m1$lrt, "df(m1)=", m1$df, "\n")

## ----eval=FALSE---------------------------------------------------------------
#  (df34 <- m3$df - m4$df)
#  (dev34 <- m3$lrt - m4$lrt)
#  1-pchisq(dev34, df34)
#  
#  
#  (df23 <- m2$df - m3$df)
#  (dev23 <- m2$lrt - m3$lrt)
#  1-pchisq(dev23, df23)
#  
#  (df12 <- m1$df - m2$df)
#  (dev12 <- m1$lrt - m2$lrt)
#  1-pchisq(dev12, df12)

## ----eval=FALSE---------------------------------------------------------------
#  m2.GM <- loglin(table1, margin=list(c(1,3), c(1,4), c(2,3), c(2,4), c(3,4)),
#                  fit=TRUE)
#  cat("deviance(m2.GM)=", m2.GM$lrt, "df(m2.GM)=", m2.GM$df, "\n")
#  
#  m2.MC <- loglin(table1, margin=list(c(1,2), c(1,3), c(1,4), c(2,4), c(3,4)),
#                  fit=TRUE)
#  cat("deviance(m2.MC)=", m2.MC$lrt, "df(m2.MC)=", m2.MC$df, "\n")
#  
#  m2.CI <- loglin(table1, margin=list(c(1,2), c(1,3), c(1,4), c(2,3), c(2,4)),
#                  fit=TRUE)
#  cat("deviance(m2.CI)=", m2.CI$lrt, "df(m2.CI)=", m2.CI$df, "\n")
#  
#  m2.GI <- loglin(table1, margin=list(c(1,2), c(1,3), c(2,3), c(2,4), c(3,4)),
#                  fit=TRUE)
#  cat("deviance(m2.GI)=", m2.GI$lrt, "df(m2.GI)=", m2.GI$df, "\n")
#  
#  m2.GC <- loglin(table1, margin=list(c(1,2), c(1,4), c(2,3), c(2,4), c(3,4)),
#                  fit=TRUE)
#  cat("deviance(m2.GC)=", m2.GC$lrt, "df(m2.GC)=", m2.GC$df, "\n")
#  
#  m2.MI <- loglin(table1, margin=list(c(1,2), c(1,3), c(1,4), c(2,3), c(3,4)),
#                  fit=TRUE)
#  cat("deviance(m2.MI)=", m2.MI$lrt, "df(m2.MI)=", m2.MI$df, "\n")

## ----eval=FALSE---------------------------------------------------------------
#  1 - pchisq(m2.GM$lrt - m2$lrt, 1)
#  
#  1 - pchisq(m2.MC$lrt - m2$lrt, 1)
#  
#  1 - pchisq(m2.CI$lrt - m2$lrt, 1)
#  
#  1 - pchisq(m2.GI$lrt - m2$lrt, 1)
#  
#  1 - pchisq(m2.GC$lrt - m2$lrt, 1)
#  
#  1 - pchisq(m2.MI$lrt - m2$lrt, 1)

## ----eval=FALSE---------------------------------------------------------------
#  m2.GM.GI.GC<- loglin(table1, margin=list(c(1), c(2,3), c(2,4), c(3,4)), fit=TRUE)
#  cat("deviance(m2.GM.GI.GC)=", m2.GM.GI.GC$lrt, "df(m2.GM.GI.GC)=", m2.GM.GI.GC$df,
#      "\n")
#  
#  1 - pchisq(m2.GM.GI.GC$lrt - m2$lrt, m2.GM.GI.GC$df - m2$df)

## ----eval=FALSE---------------------------------------------------------------
#  m2.G<- loglin(table1, margin=list(c(2,3), c(2,4), c(3,4)), fit=TRUE)
#  cat("deviance(m2.G)=", m2.G$lrt, "df(m2.G)=", m2.G$df, "\n")
#  
#  1 - pchisq(m2.G$lrt - m2$lrt, m2.G$df - m2$df)