Python SciPy special.fdtr用法及代码示例

本文简要介绍 python 语言中 scipy.special.fdtr 的用法。

用法:

scipy.special.fdtr(dfn, dfd, x, out=None) = <ufunc 'fdtr'>#

F累积分布函数。

返回 F-distribution 的累积分布函数的值，也称为 Snedecor 的 F-distribution 或 Fisher-Snedecor 分布。

带有参数\(d_n\)和\(d_d\)的F-distribution是随机变量的分布，

\[X = \frac{U_n/d_n}{U_d/d_d},\]

其中 \(U_n\) 和 \(U_d\) 是随机变量分布 \(\chi^2\) ，分别具有 \(d_n\) 和 \(d_d\) 自由度。

参数 ：：

dfn： array_like

第一个参数(正浮点数)。

dfd： array_like

第二个参数(正浮点数)。

x： array_like

参数(非负浮点数)。

out： ndarray，可选

函数值的可选输出数组

返回 ：：

y： 标量或 ndarray

F-distribution 的 CDF，参数 dfn 和 dfd 在 x 处。

注意：

采用正则化不完全beta函数，根据公式，

\[F(d_n, d_d; x) = I_{xd_n/(d_d + xd_n)}(d_n/2, d_d/2).\]

Cephes [1] 例程的包装器 fdtr 。 F 分布也可用作 scipy.stats.f 。与 scipy.stats.f 的 cdf 方法相比，直接调用 fdtr 可以提高性能(请参见下面的最后一个示例)。

参考：

[1]

Cephes 数学函数库，http://www.netlib.org/cephes/

例子：

计算 x=1 处的 dfn=1 和 dfd=2 的函数。

>>> import numpy as np
>>> from scipy.special import fdtr
>>> fdtr(1, 2, 1)
0.5773502691896258

通过为 x 提供 NumPy 数组来计算多个点的函数。

>>> x = np.array([0.5, 2., 3.])
>>> fdtr(1, 2, x)
array([0.4472136 , 0.70710678, 0.77459667])

绘制多个参数集的函数。

>>> import matplotlib.pyplot as plt
>>> dfn_parameters = [1, 5, 10, 50]
>>> dfd_parameters = [1, 1, 2, 3]
>>> linestyles = ['solid', 'dashed', 'dotted', 'dashdot']
>>> parameters_list = list(zip(dfn_parameters, dfd_parameters,
...                            linestyles))
>>> x = np.linspace(0, 30, 1000)
>>> fig, ax = plt.subplots()
>>> for parameter_set in parameters_list:
...     dfn, dfd, style = parameter_set
...     fdtr_vals = fdtr(dfn, dfd, x)
...     ax.plot(x, fdtr_vals, label=rf"$d_n={dfn},\, d_d={dfd}$",
...             ls=style)
>>> ax.legend()
>>> ax.set_xlabel("$x$")
>>> ax.set_title("F distribution cumulative distribution function")
>>> plt.show()

F 分布也可用作 scipy.stats.f 。直接使用 fdtr 比调用 scipy.stats.f 的 cdf 方法要快得多，特别是对于小型数组或单个值。为了获得相同的结果，必须使用以下参数化：stats.f(dfn, dfd).cdf(x)=fdtr(dfn, dfd, x)。

>>> from scipy.stats import f
>>> dfn, dfd = 1, 2
>>> x = 1
>>> fdtr_res = fdtr(dfn, dfd, x)  # this will often be faster than below
>>> f_dist_res = f(dfn, dfd).cdf(x)
>>> fdtr_res == f_dist_res  # test that results are equal
True