Anaconda-Spyder中使用Scipy拟合F分布时遇RuntimeWarning问题

最近我在Anaconda-Spyder里碰到个挺奇怪的问题，因为要做F分布拟合，就从Scipy里导入了F分布相关工具，写了下面这段代码：

import numpy as np
from scipy.stats import f
a = f.rvs(29, 18, size = 1000) # 官方文档中使用29和18作为参数
f.fit(a)

运行之后收到了这样的RuntimeWarning：

C:\Users\rpaca\Anaconda3\lib\site-packages\scipy\stats_distn_infrastructure.py:2303: RuntimeWarning: invalid value encountered in double_scalars
Lhat = muhat - Shat*mu

这警告看得我有点困惑，明明是照着官方文档的参数生成的样本，怎么拟合的时候就出现无效值了呢？

内容的提问来源于stack exchange，提问作者Bobesh