A sondagem para as eleições europeias

Estimados amigos.
Bem sei que há muito tempo que não escrevo mas tal deveu-se ao trabalho mas também a algo de que não quero falar por agora.
Por agora vou falar da sondagem da RTP sobre as europeias.
Pegando nos dados dessa sondagem e fazendo 10000 re-amostragens, os resultados foram:

1 = Número de deputados mais provável
   PS = 9 (86,3%)  = Aumenta um
   PSD = 6 (91,9%) = Mantém
   BE = 2 ( 98,7%) = Aumenta um
   CDU = 2 ( 95,7%) = Perde um
   CDS = 2 (95,7%)  = Aumenta um
   PAN = 0 (93,5%) = Mantém zero
   Aliança = 0 (93,5%) = Mantém zero
   Marinho e Pinho = 0 = Perde dois (a sondagem não explicitou este caso)

2 - As probabilidades detalhadas
   PS = 8 (11,0%) // 9 (86,3%) // 10 (2,8%)
   PSD = 5 (2,8%) // 6 (91,9%) // 7 (5,2%)
   BE = 1 (0,2%) // 2 ( 98,7%) // 3 (1,1%)
   CDU = 1 (4,3%) // 2 (95,7%)
   CDS = 1 (4,3%) // 2 (95,7%) 
   PAN = 0 (93,5%) // 1 (6,5%)
  Aliança = 0 (93,5%) // 1 (6,5%)

#O código que utilizei
dados<-c(621, 433, 169, 151, 151, 56, 56, 245)
perc<-dados/sum(dados)
D<-matrix(1,7,1,FALSE)
for (j in 1:10000)
 {resul<-sample(1:8,1882,prob=perc,replace=TRUE)
 resul.t<-table(resul)[1:7]
 deputados<-rep(1,7)
 for (i in 1:21)
  {resul.x <- resul.t/deputados
  x<-as.numeric(which.max(resul.x))
  deputados[x]<-deputados[x]+1
  }
 deputados<-deputados-1
 D<-matrix(c(D,deputados),7,j+1,FALSE)
 }
D<-D[,2:ncol(D)]
for (i in 1:7)
 print(table(D[i,])/ncol(D))

Read more »