F分位数自由度转换的直观理解

为了更好理解这个转换自由度之后，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,10,4)
plt.plot(x,y1,'r',label='y=f(x,4,10)')
plt.plot(x,y2,'b',label='y=f(x,10,4)')
plt.fill_between(x,0,y1,where=(x<0.38314),facecolor='red',alpha=0.3)
plt.fill_between(x,0,y2,where=(x>2.61),facecolor='blue',alpha=0.3)
plt.legend()
plt.show()