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

## ----eval=FALSE---------------------------------------------------------------
#  partypref <- matrix(data=c(114, 10, 53,224,134,9,42,226,114,8,23,174,339,30,13,
#  414,42,5,44,161,88,10, 60,171,90,8,31,168,413,23,14,375), nrow=8, byrow=TRUE)
#  
#  
#  partydat<-data.frame(
#  party=c(rep("CDU",sum(partypref[,1])),rep("SPD",sum(partypref[,4])),
#  rep("The Liberals",sum(partypref[,2])),rep("The Greens",sum(partypref[,3]))),
#  sex=c(rep(0,sum(partypref[1:4,1])),rep(1,sum(partypref[5:8,1])),
#  rep(0,sum(partypref[1:4,4])),rep(1,sum(partypref[5:8,4])),
#  rep(0,sum(partypref[1:4,2])),rep(1,sum(partypref[5:8,2])),
#  rep(0,sum(partypref[1:4,3])),rep(1,sum(partypref[5:8,3]))),
#  age=c(rep(c(1:4,1:4), partypref[,1]),rep(c(1:4,1:4), partypref[,4]),
#  rep(c(1:4,1:4), partypref[,2]),rep(c(1:4,1:4), partypref[,3])))
#  

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

## ----eval=FALSE---------------------------------------------------------------
#  datmat<-as.matrix(table(partydat$sex,partydat$party))
#  tparty<-data.frame("CDU"=datmat[,1],"SPD"=datmat[,2],"Green"=datmat[,3],
#  "Liberals"=datmat[,4],"sex"=0:1)
#  tparty
#  
#  logitParty <- multinom(cbind(CDU,SPD,Green,Liberals)~sex, data=tparty)
#  
#  summary(logitParty)
#  exp(coef(logitParty))

## ----eval=FALSE---------------------------------------------------------------
#  coefSPD <- matrix(data = c(-coefficients(logitParty)[3,1],
#  coefficients(logitParty)[1,1] - coefficients(logitParty)[3,1],
#  coefficients(logitParty)[2,1] - coefficients(logitParty)[3,1],
#  -coefficients(logitParty)[3,2],
#  coefficients(logitParty)[1,2] - coefficients(logitParty)[3,2],
#  coefficients(logitParty)[2,2] - coefficients(logitParty)[3,2]),
#  nrow=3, ncol=2)
#  
#  coefSPD
#  
#  exp(coefSPD)