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

## ----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---------------------------------------------------------------
#  x1 <- partypref/rowSums(partypref)
#  
#  x1 <- x1[c(1,4,5,8),]
#  
#  x1 <- cbind(x1[,2],x1[,3],x1[,1],x1[,4])
#  
#  x1 <- x1*6

## ----eval=FALSE---------------------------------------------------------------
#  loc1 <- matrix(data=c(0,0,15,0,0,15,15,15),ncol=2,byrow=T)
#  library(grDevices)

## ----eval=FALSE---------------------------------------------------------------
#  stars(x1, scale=FALSE,key.loc = c(25,6),len=2,cex=1.2,lwd=1.2,
#  xlim=c(-10,30),ylim=c(-5,10),key.labels=c("Liberals","Green","CDU","SPD"),
#  location=loc1, labels = c("male young","male old","female young","female old"))

## ----eval=FALSE---------------------------------------------------------------
#  partydat$age <- as.factor(partydat$age)
#  
#  library(nnet)
#  
#  partymult <- multinom(party ~ sex + age , data=partydat)
#  
#  summary(partymult)

## ----eval=FALSE---------------------------------------------------------------
#  x2 <- matrix(data=c(rep(1,5),as.matrix(t(exp(coefficients(partymult))))),nrow=5)
#  
#  x2 <- cbind(x2[,4],x2[,3],x2[,1],x2[,2])
#  
#  loc2 <- matrix(data=c(0,9,9,9,0,0,9,0,18,0),ncol=2,byrow=T)

## ----eval=FALSE---------------------------------------------------------------
#  stars(x2, scale=FALSE,key.loc = c(18,9),len=2,cex=1.2,lwd=1.2,
#  xlim=c(-2,23),ylim=c(-3,10),key.labels=c("Liberals","Green","CDU","SPD"),
#  location=loc2, labels = c("reference","gender","age 2","age 3","age 4"))