Pág. 40( parte inicial da página).
Exemplo de programa em Python para explorar a noção de intervalo de confiança. Extraem-se sucessivas amostras, com a mesma dimensão, de uma população com distribuição Normal N(0,1) e para cada uma delas apresenta-se o intervalo de 95% de confiança para o valor médio. Pode verificar-se que este intervalo depende da amostra e que nem sempre contém o valor médio da população. Pode também fazer-se variar a dimensão das amostras e comparar as amplitudes dos intervalos.
Neste caso, fazemos 100 amostras aleatórias de dimensão 20.
(...)
Comentários/Sugestões:
Extraem-se 100 amostras, com a mesma dimensão 20, de uma população com distribuição Normal N(0,1) e para cada uma delas apresenta-se o intervalo de 95% de confiança para o valor médio. (Recorde a fórmula do intervalo de confiança)
Pode verificar-se que este intervalo depende da amostra e que nem sempre contém o valor médio da população.
Pode também fazer-se variar a dimensão das amostras e comparar as amplitudes dos intervalos.
|
Importamos os aleatórios, a matemática e a estatística.
import random
import statistics
import math
|
Fixamos 100 amostras
num_amostras=100
|
Também fixamos que a dimensão de cada amostra é 20
dim_amostra=20
|
Para cada uma das amostras,
for i in range(num_amostras):
|
Chamamos amostra a um conjunto de elementos aleatórios n(0,1)
amostra=[random.normalvariate(0,1) for k in range(dim_amostra)]
|
Calculamos o valor médio e o desvio padrão de cada amostra.
vmedio=statistics.mean(amostra)
desvio_padrao=statistics.stdev(amostra)
|
Chamamos d à habitual margem de erro do intervalo de confiança.
d=1.96*desvio_padrao/math.sqrt(dim_amostra)
|
Apresentamos os 100 intervalos pretendidos.
print('O intervalo de 95% de confiança para o valor médio é [',vmedio-d,',',vmedio+d ,']')
Sugestão-Consulte:
https://pedronoia.pt/python/pyl1.htm
|
https://pedronoia.pt/python/pyl2.htm
|
https://pedronoia.pt/python/pyl3.htm
|