terça-feira, 21 de maio de 2019

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))

4 comentários:

Anónimo disse...

Meu caro, força. A vida continua

Anónimo disse...

Bem regreesado

Anónimo disse...

volte sempre

Henrique Silva disse...

Regresse rápido Sr. Professor. As suas "crónicas" são muito instrutivas. Um abraço

Enviar um comentário

Twitter Delicious Facebook Digg Stumbleupon Favorites More

 
Design by Free WordPress Themes | Bloggerized by Lasantha - Premium Blogger Themes | Best Hostgator Coupon Code