
80
6 章 確率密度関数
>>> pmf = pdf.MakePmf()
デフォルトでは、結果の
Pmf
は、
mu-3*sigma
から
mu+3*sigma
まで均等にばらま
かれた
101
個の点を含みます。
MakePmf
と
Render
は、
low
,
high
,
n
というキーワード
引数をオプションとして取ることができます。
6.2
カーネル密度推定
カーネル密度推定(
Kernel density estimation
、
KDE
)は、標本を入力として、
データに合致する近似的な円滑
PDF
を求めるアルゴリズムです。詳細は、
http://
ja.wikipedia.org/wiki/
カーネル密度推定
†
で読めます。
scipy
は
KDE
の実装を提供し、
thinkstats2
はその
KDE
を使う
EstimatedPdf
ク
ラスを提供します。
class EstimatedPdf(Pdf):
def __init__(self, sample):
self.kde = scipy.stats.gaussian_kde(sample)
def Density(self, x):
return self.kde.evaluate(x)
__init__
は標本を取って、カーネル密度推定を計算します。結果は
evaluate
メソッ
ドを提供する
gaussian_kde
オブジェクトです。
Density
は値かシーケンスを取り、
gaussian_kde.evaluate
を呼び出し、結果とし
て密度を返します。名前に「
Gaussian
」という語が含まれるのは、 ...