In [12]:
%matplotlib inline
import pandas as pd
import matplotlib.pyplot as plt
import numpy as np
import os
from scipy import stats
lambp=10
dist= stats.poisson(lambp)
x= stats.poisson.rvs(mu=lambp, loc=0, size=1000000)
media= np.mean(x)
var= np.var(x)
In [23]:
media2=[]
for i in range (0,10001):
valores = np.random.choice(x,2)
media2.append(np.mean(valores))
plt.hist(media2)
plt.title("n valendo 2")
media = np.mean(media2)
variancia = np.var(media2)
print("Media:{}".format(media))
print("Variancia:{}".format(variancia))
In [15]:
mediaVal=[]
for i in range (0,10001):
valores = np.random.choice(x,4)
mediaVal.append(np.mean(valores))
plt.hist(mediaVal)
media = np.mean(mediaVal)
plt.title("n valendo 4")
variancia = np.var(mediaVal)
print("Media:{}".format(media))
print("Variancia:{}".format(variancia))
In [21]:
media20=[]
for i in range (0,10001):
valores = np.random.choice(x,20)
media20.append(np.mean(valores))
plt.hist(media20)
plt.title("n valendo 20")
media = np.mean(media20)
variancia = np.var(media20)
print("Media:{}".format(media))
print("Variancia:{}".format(variancia))
In [28]:
media50=[]
for i in range (0,10001):
valores = np.random.choice(x,50)
media50.append(np.mean(valores))
plt.hist(media50)
plt.title("n valendo 50")
media = np.mean(media50)
variancia = np.var(media50)
print("Media:{}".format(media))
print("Variancia:{}".format(variancia))
In [45]:
media50=[]
primeiro= []
nono = []
for i in range (0,10001):
valores = np.random.choice(x,50)
primeiro.append(valores[0])
nono.append(valores[9])
media50.append(np.mean(valores))
media = np.mean(media50)
variancia = np.var(media50)
print("Media:{}".format(media))
print("Variancia:{}".format(variancia))
plt.hist(media50)
plt.title("n valendo 50")
plt.show()
dp1= (np.var(primeiro))**1/2
dp9=(np.var(nono))**1/2
cov = np.cov(primeiro,nono)
plt.hist(primeiro)
plt.title("primeiro")
plt.show()
plt.hist(nono)
plt.title("nono")
plt.show()
print(cov/(dp1*dp9))
https://www.youtube.com/watch?v=YGhL0ltpuo0
In [ ]: