正規分布の積分計算について

Question

実現したいこと

正規分布における積分計算を行いたい

前提

scipyを用いて正規分布N(1, 1)の区間[0,∞]の積分を行いたいです。

発生している問題・エラーメッセージ

ValueError: invalid callable given

該当のソースコード

import math
import numpy as np
from scipy import stats
from scipy import integrate

x = np.arange(start = -1, stop = 3, step = 0.01)
y = stats.norm.pdf(x = x, loc = 1, scale = 1)
s1 = integrate.quad(y, 0, delta_mu)
S  = 0.5 + s1

Comments

[melian]

おそらく以下を、

x = np.arange(start = -1, stop = 3, step = 0.01)
y = stats.norm.pdf(x = x, loc = 1, scale = 1)
s1 = integrate.quad(y, 0, delta_mu)

こちらに変更することになりますが、これは期待する結果でしょうか。

y = lambda x: stats.norm.pdf(x = x, loc = 1, scale = 1)
s1 = integrate.quad(y, 0, np.inf)

[geant-kami]

できました、ありがとうございます！

Answer 1

cdf(累積密度関数)はpdf(確率密度関数)の積分関数なので，
pdfの0から∞までの積分はcdf(∞) - cdf(0)となり次のように計算できるのでは?

import numpy as np
from scipy import stats

pdf_sum_0_inf = stats.norm.cdf(np.inf, loc=1, scale=1) \
    - stats.norm.cdf(0, loc=1, scale=1)

Comments

[geant-kami]

解決しました。ありがとうございます！