UA MATH564 概率论VI 数理统计基础5 F分布

UA MATH564 概率论VI 数理统计基础5 F分布

假设 X χ m , δ 2 , Y χ n 2 X \sim \chi^2_{m,\delta},Y \sim \chi^2_{n} 且二者独立,则
Z = X / m Y / n F m , n , δ Z = \frac{X/m}{Y/n} \sim F_{m,n,\delta}
其中 δ \delta 被称为非中心参数, δ = 0 \delta=0 则称 Z Z 服从中心化的F分布,简称F分布。注意到 Z Z 也是两个随机变量的商,可以用上一讲介绍的引理1计算 Z Z 的概率密度,此处省略过程:
f ( x m , n , δ ) = e δ 2 2 i = 0 ( δ 2 2 ) i i ! n n 2 m m 2 + i c i x m 2 1 + i ( n + m x ) m + n 2 + i c i = Γ ( m + n 2 + i ) Γ ( n 2 ) Γ ( m 2 + i ) f(x|m,n,\delta) = e^{-\frac{\delta^2}{2}}\sum_{i=0}^{\infty} \frac{(\frac{\delta^2}{2})^i}{i!}n^{\frac{n}{2}}m^{\frac{m}{2}+i}c_i \frac{x^{\frac{m}{2}-1+i}}{(n+mx)^{\frac{m+n}{2}+i}} \\ c_i = \frac{\Gamma(\frac{m+n}{2}+i)}{\Gamma(\frac{n}{2})\Gamma(\frac{m}{2}+i)}

δ = 0 \delta=0 时,概率密度为
f ( x m , n ) = Γ ( m + n 2 ) Γ ( n 2 ) Γ ( m 2 ) n n 2 m m 2 x m 2 1 ( n + m x ) m + n 2 f(x|m,n) = \frac{\Gamma(\frac{m+n}{2})}{\Gamma(\frac{n}{2})\Gamma(\frac{m}{2})}n^{\frac{n}{2}}m^{\frac{m}{2}}\frac{x^{\frac{m}{2}-1}}{(n+mx)^{\frac{m+n}{2}}}
下面介绍三条性质:

性质1: t n , δ = d F 1 , n , δ t_{n,\delta} =_d F_{1,n,\delta}

性质2: X i N ( a , σ 2 ) , Y j N ( b , σ 2 ) , i = 1 , , m , j = 1 , , n X_i \sim N(a,\sigma^2),Y_j \sim N(b,\sigma^2),i=1,\cdots,m,j=1,\cdots,n 所有的 X , Y X,Y 均互相独立,则
Z = 1 m 1 i = 1 n ( X i X ˉ ) 2 1 n 1 j = 1 n ( Y j Y ˉ ) 2 F m 1 , n 1 Z = \frac{\frac{1}{m-1}\sum_{i=1}^n (X_i - \bar{X})^2}{\frac{1}{n-1}\sum_{j=1}^n (Y_j - \bar{Y})^2} \sim F_{m-1,n-1}
性质3:记 X n F m , n , δ X_n \sim F_{m,n,\delta} ,则 X n d 1 m χ m , δ 2 X_n \to_d \frac{1}{m}\chi^2_{m,\delta}

性质1根据概率密度的表达式就可以证明,性质2也可以根据定义写出来。性质1的意义是说明在线性模型中,针对单个系数做t检验和做partial F检验是等价的;性质2的意义是给双正态总体的方差检验提供依据。下面证明性质3:

证明
根据定义,可以将 X n X_n 写成
X n = K / m Y / n F m , n , δ X_n = \frac{K/m}{Y/n} \sim F_{m,n,\delta}

K χ m , δ 2 , Y χ n 2 K \sim \chi^2_{m,\delta},Y \sim \chi^2_{n} 且二者独立。其中 Y / n Y/n 又可以写成
Y n = 1 n j = 1 n Z j 2 , Z j i i d N ( 0 , 1 ) \frac{Y}{n} = \frac{1}{n}\sum_{j=1}^n Z_j^2,Z_j \sim_{iid} N(0,1)

根据大数定律,
Y n = 1 n j = 1 n Z j 2 P E [ Z 1 2 ] = 1 \frac{Y}{n} = \frac{1}{n}\sum_{j=1}^n Z_j^2 \to_P E[Z_1^2] = 1

因此 X n d K / m 1 m χ m , δ 2 X_n \to_d K/m \sim \frac{1}{m}\chi^2_{m,\delta}
证毕

猜你喜欢

转载自blog.csdn.net/weixin_44207974/article/details/106825831