2016-03-19 30 views
2

私はscipy.stats.rv_continuousのサブクラスを初期化する場合は、例えばscipy.stats.norm連続ランダム変数をscipyで離散化することはできますか?

>>> from scipy.stats import norm 
>>> rv = norm() 

私は範囲の数を提供した後、値の範囲の確率を表す各要素を持つ確率のリストに変換することはできますか? (範囲について - [(-INF、-1)、( - 1,0)、(0,1)、(1、INF)]) - 様 何か

>>> li 
[0.15865525393145707, 0.34134474606854293, 0.34134474606854293, 0.15865525393145707] 

(0.15865525393145707であり範囲内にあるため-1未満と0.34134474606854293 -1 0にし、同様に他の人のためという変数の確率が。

これは、scipyのダウンロードを使用して行うことができますか?そうでない場合はPythonライブラリは、このような変換操作をサポートすることができた?

答えて

2

あなたのコメントに基づいて、CDF

を使用してこれを計算できます
from scipy.stats import norm 
import numpy as np 

>>> norm().cdf(-1) - norm().cdf(-np.inf), \ 
    norm().cdf(0) - norm().cdf(-1), \ 
    norm().cdf(1) - norm().cdf(0), \ 
    norm().cdf(np.inf) - norm().cdf(1) 
(0.15865525393145707, 
0.34134474606854293, 
0.34134474606854293, 
0.15865525393145707) 

これは基本的にCDFの定義に従います。


私は数字の合計が1になることに注意してください。ただし、期待する出力として書き込む数字はありません。私はそれらが正しいものだと言ってあなたの基礎を知らない。私の推測では、非単位の標準偏差を持つNormal変数を暗黙的に使用しています。

+0

私は、次の範囲[(-inf、-1)、(-1,0)、(0,1)、(1、inf)]の確率を求めています。 0.49,0.49,0.01]であり、0.01は変数が-1より小さい確率であり、0.49は-1から0の範囲にあり、同様に他のものである。私はその範囲の総合確率を求めています。上記の方法で総合確率が得られますか? –

+0

@YashuSeth更新を参照してください。 –

+0

ありがとうございます。わかった。そして、私は例を使って説明しました。はい、正しい結果ではありません。 –

関連する問題