아래 문제를 풀어서 제출하시오.
import math
import scipy
import scipy.stats
p=0.1
n=1000
np= n * p
print("Average = ", np)
sd = math.sqrt(n * p * (1-p))
print("Stddev = ", sd)
[avg, v] = scipy.stats.binom.stats(n=1000, p=0.1)
print(avg, math.sqrt(v))
#Average = 100.0
#Stddev = 9.486832980505138
#100.0 9.486832980505138
import scipy.stats as stats
print("P[X>=30] = ", 1- stats.norm.cdf(30, loc=20, scale=5))
print("P[X<15] = ", stats.norm.cdf(15, loc=20, scale=5))
print("P[10<X<15] = ", stats.norm.cdf(25, loc=20, scale=5) - stats.norm.cdf(10, loc=20, scale=5))
#P[X>=30] = 0.02275013194817921
#P[X<15] = 0.15865525393145707
#P[10<X<15] = 0.8185946141203637
print("P[X<=5] = ", stats.poisson.cdf(5, mu=10))
print("P[7<=X<14] = ", stats.poisson.cdf(13, mu=10) - stats.poisson.cdf(6, mu=10))
print("P[X>=9] = ", 1-stats.poisson.cdf(8, mu=10) )
#P[X<=5] = 0.06708596287903189
#P[7<=X<14] = 0.7343230017368279
#P[X>=9] = 0.6671803212492808
stats.expon.cdf(1/6, scale=1/10)
#0.8111243971624381