2b5abd2e6175017

http://nhiro.org.s3.amazonaws.com/9/0/90331315436b070e267185f0c511e1a5.jpg https://gyazo.com/90331315436b070e267185f0c511e1a5

(OCR text)

2次元ガウス分布の負の対数度

def nll2d (args, group)

(mu1, mu2, sig11, sig12, sig22)

det np.linalg.det (Sigma)

if det 0.0: return np.inf

invsigma

np.log (det) len (g roup) 2

for y in group:

diff

args

np. linalg.inv (Sigma)

ret =

y - mu

ret np.dot (np.dot (diff. transpose (), invsigma), diff) 2

if np.isnan (ret): return np.inf

return ret

wilo la

7

式をそのまんまコードにしてる:

2

毎回データを紙めるなんて遅そう?でもとりあえずやってみる