Как я могу найти исходный код для scipy.stats.mvn?

Я использовал функцию mvn.mvnun() в scipy.stats Модуль для расчета CDF заданного многомерного нормального распределения в моей статье. Рецензент спросил меня, как я оцениваю CDF, так как у CDF нет закрытой формы. Я думаю, что метод выборки может быть использован. Чтобы узнать, как это работает, я искал scipy исходный код репо. Однако в scipy.stats папка, я не видел 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?

Источник

1 ответ

Решение

Реализация написана на Фортране. Тот pyf файл определяет интерфейс для функции Fortran в mvndst.f

Источник
Другие вопросы по тегам