Bootstrap pe perechi

Gândire statistică în Python (Partea 2)

Justin Bois

Lecturer at the California Institute of Technology

Inferență neparametrică

Nu se fac ipoteze despre model sau distribuția de probabilitate a datelor

Alegerile din 2008 în statele decisive din SUA

Rezultatele alegerilor din 2008 în statele decisive din SUA

¹ Date preluate de pe Data.gov (https://www.data.gov/)

Bootstrap pe perechi pentru regresia liniară

Reesantionați datele în perechi
Calculați panta și interceptul din datele reesantionate
Fiecare pantă și intercept reprezintă un replicat bootstrap
Calculați intervalele de încredere din percentilele replicatelor bootstrap

Generarea unui eșantion bootstrap pe perechi

np.arange(7)

array([0, 1, 2, 3, 4, 5, 6])

inds = np.arange(len(total_votes))

bs_inds = np.random.choice(inds, len(inds))

bs_total_votes = total_votes[bs_inds]
bs_dem_share = dem_share[bs_inds]

Calcularea unui replicat bootstrap pe perechi

bs_slope, bs_intercept = np.polyfit(bs_total_votes, 
                                    bs_dem_share, 1)

bs_slope, bs_intercept

(3.9053605692223672e-05, 40.387910131803025)

np.polyfit(total_votes, dem_share, 1)  # fit of original

array([  4.03707170e-05,   4.01139120e+01])

Alegerile din 2008 în statele decisive din SUA

Rezultatele alegerilor din 2008 în statele decisive din SUA

¹ Date preluate de pe Data.gov (https://www.data.gov/)

Să exersăm!

Gândire statistică în Python (Partea 2)