为了更好理解这个转换自由度之后,F与1/F之间分位数的关系,做出了F(4,10)和F(10,4)的图像,当置信度α=0.1时,F(4,10)=2.61,1/F(4,10)=0.38314,所以红色阴影区域和蓝色阴影区域面积相等。
代码实现如下:
import numpy as np
from scipy.stats import f
import matplotlib.pyplot as plt
x = np.linspace(0,4,400)
y1 = f.pdf(x,4,10)
y2 = f.pdf(x