
200
■
15
章
次元を扱う
サンプル数
図 15-5 シミュレーションした検体B1242 の希薄化曲線
15.11
ジョイント事後確率
これらのシミュレーションを使って、
num_new
と
k
のジョイント分布を推定して、
k
の
任意の値を条件とする
num_new
の分布を求めることができる。
def MakeJointPredictive(curves):
joint = thinkbayes.Joint()
for curve in curves:
for k, num_new in curve:
joint.Incr((k, num_new))
joint.Normalize()
return joint
MakeJointPredictive
は、値がタプル(
tuple
)の
Pmf
であるジョイントオブジェクトを
作る。
curves
は、
RunSimulation
で作られた希薄化曲線のリストである。各曲線には、
k
と
num_new
の対のリストが含まれる。