非统计学家试图用Python复制一些代码。你知道吗
R具有函数qchisq
。你知道吗
qchisq(c(0.5, 0.1, 0.99), df=1, lower.tail=F)
# [1] 0.4549364231 2.7055434541 0.0001570879
它可以在Python中复制,如下所示:
from scipy.special import chdtri
import pandas as pd
chdtri(1, pd.Series([0.5, 0.1, 0.99]))
# Out[57]:
# 0 0.454936
# 1 2.705543
# 2 0.000157
但是,Rsqchisq
也可以有标志log.p
,它允许我执行以下操作:
qchisq(-c(0.5, 0.1, 0.99) * log(10), df=1, lower.tail=F, log.p=T)
# [1] 1.00448472 0.06796154 2.66885996
我在Python中找不到任何方法来实现这一点。我怎样才能达到同样的效果?你知道吗
也许在不设置R参数的情况下重新生成
log.p
的效果是有用的:对于
log(10)
部分:看来
qchisq(exp(x)) == qchisq(x, log.p = TRUE)
。你知道吗相关问题 更多 >
编程相关推荐