Bootstrap em pares com dados de alfabetização/fertilidade
Usando a função que você acabou de escrever, faça bootstrap em pares para plotar um histograma que descreve a estimativa do coeficiente angular a partir dos dados de analfabetismo/fertilidade. Também informe o intervalo de confiança de 95% para o coeficiente angular. Os dados estão disponíveis para você nos arrays do NumPy illiteracy e fertility.
Relembrando: draw_bs_pairs_linreg() tem a assinatura draw_bs_pairs_linreg(x, y, size=1) e retorna dois valores: bs_slope_reps e bs_intercept_reps.
Este exercício faz parte do curso
Pensamento Estatístico em Python (Parte 2)
Instruções do exercício
- Use sua função
draw_bs_pairs_linreg()para obter1000réplicas bootstrap do coeficiente angular e do intercepto. Os dados do eixo x sãoilliteracye os do eixo y sãofertility. - Calcule e imprima o intervalo de confiança bootstrap de 95% para o coeficiente angular.
- Plote e mostre um histograma das réplicas do coeficiente angular. Não se esqueça de rotular seus eixos. Isso já foi feito para você, então clique em Enviar para ver seu histograma!
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Generate replicates of slope and intercept using pairs bootstrap
bs_slope_reps, bs_intercept_reps = ____
# Compute and print 95% CI for slope
print(np.percentile(____, ____))
# Plot the histogram
_ = plt.hist(bs_slope_reps, bins=50, normed=True)
_ = plt.xlabel('slope')
_ = plt.ylabel('PDF')
plt.show()