Question

我正在使用mvn.mvnun()模块中的函数scipy.stats计算我的论文中给定多元正态分布的CDF。审阅者问我，由于CDF没有封闭格式，我如何估算CDF。我猜可能会使用抽样方法。为了了解其工作原理，我搜索了scipy源代码存储库。但是，在scipy.stats folder中，我没有看到mvn.py，而是看到了mvn.pyf。

在文件中，mvn.mvnun的定义如下，

!    -*- f90 -*-
! Note: the context of this file is case sensitive.

python module mvn ! in 
    interface  ! in :mvn
        subroutine mvnun(d,n,lower,upper,means,covar,maxpts,abseps,releps,value,inform) ! in :mvn:mvndst.f
            integer intent(hide) :: d=shape(means,0)
            integer intent(hide) :: n=shape(means,1)
            double precision dimension(d) :: lower
            double precision dimension(d) :: upper
            double precision dimension(d,n) :: means
            double precision dimension(d,d) :: covar
            integer intent(optional) :: maxpts=d*1000
            double precision intent(optional) :: abseps=1e-6
            double precision intent(optional) :: releps=1e-6
            double precision intent(out) :: value
            integer intent(out) :: inform
        end subroutine mvnun

但是，没有该功能的详细定义。我想知道在哪里可以找到mvnun的源代码，它告诉我如何计算CDF？

Answer 1

该实现是用Fortran编写的。该pyf文件在mvndst.f

中定义了Fortran函数的接口。

如何找到scipy.stats.mvn的源代码？

1 个答案: