Pravdepodobnostný prieskum: septembrový FOCUS II

Prieskum verejnej mienky nie je statický. Ukazuje nám vzorku z celej populácie, ktorú sme nemuseli vybrať tak, že dokonale reprezentuje všetkých. Má teda chybu. Štatistická odchýlka vo výstupoch od agentúr je arbitrárna (bude to 95 %, 99 %?) a neukazuje na pravdepodobnosť jednotlivých možností v rozptyle.

Nahodíme si výsledky prieskumu v percentách. Nesmieme zabudnúť na možnosti "nezúčastním sa" a "neviem", lebo zverejňované percentá sú iba z rozhodnutých.

N = 1017 neviem = 0.032 neucast = 1 - 0.744 - neviem viem = (1-neucast-neviem) NR = N * viem POLL = { 'neucast': neucast, 'neviem': neviem, 'HLAS': viem * 0.137, 'SMER': viem * 0.180, 'PS': viem * 0.166, 'REPUBLIKA': viem * 0.077, 'SaS': viem * 0.058, 'SME RODINA': viem * 0.041, 'OLANO': viem * 0.082, 'KDH': viem * 0.065, 'Aliancia': viem * 0.035, 'Demokrat': viem * 0.040, 'SNS': viem * 0.064, 'LSNS': viem * 0.017, 'MF': viem * 0.006, 'Modri': viem * 0.007, 'Ini': viem * 0.025 } K = len(POLL) Z = 2.576 # confidence interval for 99 %

Percentá prevedieme na počty respondentov.

Takto vyzerá rozptyl možností (modrou) pri 99% spoľahlivosti. Nevidím tam ale čo z toho je najpravdepodobnejšie ani ako sa to premietne do mandátov.

Spravíme si preto 10 000 "hodení kockou" aby sme videli čo pri takomto prieskume bude najčastejší výsledok.

Z toho vytiahneme iba rozhodnutých voličov, aby sme potom mohli počítať mandáty

Takto vyzerá pravdepodobnostná distribúcia volebných výsledkov v %

Potom to prehodíme na mandáty (s 5 % hranicou, resp. 7 % pre OLANO) a vyzerá to takto:

# OTAZKY # VYHRA PS? vyhra_ps = (MANDATY.max(axis=1) == MANDATY.PS).astype(int).sum() / ROLLS print(f'PS vyhra so sancou {vyhra_ps:.1%}') # ZOSTAVI PS VLADU BEZ HLASU? # Uvazujeme PS, SaS, Demokratov, KDH vlada_bez_hlasu = (MANDATY.PS + MANDATY.SaS + MANDATY.Demokrat + MANDATY.KDH > 75).astype(int).sum() / ROLLS print(f'PS zostavi vladu bez Hlasu so SaS, Demokratmi a KDH so sancou {vlada_bez_hlasu:.1%}') # ZOSTAVI PS VLADU S HLASOM? # Uvazujeme PS, SaS, Demokratov, KDH, Hlas vlada_s_hlasom = (MANDATY.PS + MANDATY.SaS + MANDATY.Demokrat + MANDATY.KDH + MANDATY.HLAS > 75).astype(int).sum() / ROLLS print(f'PS zostavi vladu s Hlasom, SaS, Demokratmi a KDH so sancou {vlada_s_hlasom:.1%}') # ZOSTAVI SMER VLADU S HLASOM A SNS? smer_sns_hlas = (MANDATY.SMER + MANDATY.HLAS + MANDATY.SNS > 75).astype(int).sum() / ROLLS print(f'Smer zostavi vladu s Hlasom a SNS so sancou {smer_sns_hlas:.1%}') # ZOSTAVI SMER VLADU S REPUBLIKOU A SNS? smer_sns_rep = (MANDATY.SMER + MANDATY.REPUBLIKA + MANDATY.SNS > 75).astype(int).sum() / ROLLS print(f'Smer zostavi vladu s Republikou a SNS so sancou {smer_sns_rep:.1%}') # ZOSTAVI SMER VLADU S HLASOM, SNS a KDH? smer_sns_hlas_kdh = (MANDATY.SMER + MANDATY.HLAS + MANDATY.SNS + MANDATY.KDH > 75).astype(int).sum() / ROLLS print(f'Smer zostavi vladu s Hlasom, SNS a KDH so sancou {smer_sns_hlas_kdh:.1%}') # PREJDU DEMOKRATI? demo = (MANDATY.Demokrat > 0).astype(int).sum() / ROLLS print(f'Demokrati prejdu do parlamentu so sancou {demo:.1%}') # PREJDE SNS? sns = (MANDATY.SNS > 0).astype(int).sum() / ROLLS print(f'SNS prejde do parlamentu so sancou {sns:.1%}') # PREJDE KDH? kdh = (MANDATY.KDH > 0).astype(int).sum() / ROLLS print(f'KDH prejde do parlamentu so sancou {kdh:.1%}') # PREJDE SME RODINA? smr = (MANDATY['SME RODINA'] > 0).astype(int).sum() / ROLLS print(f'Sme Rodina prejde do parlamentu so sancou {smr:.1%}') # PREJDE SaS? sas = (MANDATY.SaS > 0).astype(int).sum() / ROLLS print(f'SaS prejde do parlamentu so sancou {sas:.1%}') # PREJDE OLANO? olano = (MANDATY.OLANO > 0).astype(int).sum() / ROLLS print(f'OLANO prejde do parlamentu so sancou {olano:.1%}')

koalicie = pd.DataFrame({ 'PS+bezHlasu': (MANDATY.PS + MANDATY.SaS + MANDATY.Demokrat + MANDATY.KDH).values, 'PS+shlasom': (MANDATY.PS + MANDATY.SaS + MANDATY.Demokrat + MANDATY.KDH + MANDATY.HLAS).values, 'SMER+HLAS+SNS': (MANDATY.SMER + MANDATY.HLAS + MANDATY.SNS).values, 'PS+bezhlasu+olano': (MANDATY.PS + MANDATY.SaS + MANDATY.Demokrat + MANDATY.KDH + MANDATY.OLANO).values, 'SMER+HLAS+SNS+Rodina': (MANDATY.SMER + MANDATY.HLAS + MANDATY.SNS + MANDATY['SME RODINA']).values }) fig = ff.create_distplot(koalicie.values.T, group_labels = koalicie.columns, bin_size=1) fig.update_yaxes(range=[0,0.1]) fig.update_xaxes(range=[40,90]) fig.add_vline(x=75.) fig.show()