【学习笔记】第四章概率论与数理统计

4.1 随机变量的概率计算和数字特征

4.1.1 随机变量的概率计算

例4.1 设 $X\sim N(3,5^2)$ (1）求P{2<X<6}；（2）确定c，使P{-3c<X<2c}=0.6

from scipy.stats import norm
from scipy.optimize import fsolve
print("p=",norm.cdf(6,3,5)-norm.cdf(2,3,5))#做差，后减前
f=lambda c: norm.cdf(2*c,3,5)-norm.cdf(-3*c,3,5)-0.6
print("c=",fsolve(f,0))

例4.2 设 $X\sim N(0,1)$ ,若 $z_\alpha$ 满足条件 $P\left \{X > z_\alpha \right \}= \alpha ,0<\alpha <1$ ，则称 $z_\alpha$ 为标准正态分布的上α分位数。试计算几个常用的 $z_\alpha$ 的值,并画出 $z_{0.1}$ 的示意图.
计算得到几个常用的 $z_\alpha$ 的值见下表， $z_{0.1}$ 的示意图见下图：

from scipy.stats import norm
from pylab import plot,fill_between,show,text,savefig,rc 
from numpy import array, linspace, zeros
alpha=array([0.001, 0.005, 0.01, 0.025, 0.05, 0.10])#表中给的alpha的值
za=norm.ppf(1-alpha,0,1)  #求上alpha分位数  #定义4.2
print("上alpha分位数分别为", za)
x=linspace(-4, 4, 100); y=norm.pdf(x, 0, 1)
rc('font',size=16); rc('text',usetex=True)
plot(x,y)  #画标准正态分布密度曲线
x2=linspace(za[-1],4,100)#在0.1对应的上alpha分位数和4之间取100个点
y2=norm.pdf(x2);
y1=[0]*len(x2)
fill_between(x2, y1, y2, color='r')  #y1,y2对应的点之间填充
plot([-4,4],[0,0])  #画水平线    # [-4,0]|[4,0]
text(1.9, 0.07, "$\\leftarrow\\alpha$=0.1")  #标注
savefig("figure4_2.png", dpi=500); show()