読者です 読者をやめる 読者になる 読者になる

tak0kadaの何でもノート

発声練習、生存確認用。

医学関連は 医学ノート

KLダイバージェンス、JSダイバージェンス

分布 統計
# Kallback-Leiblerダイバージェンス
def kld(p, q):
    p = np.array(p)
    q = np.array(q)
    return np.sum(p * np.log(p/q))

# Jensen-Shannonダイバージェンス
def jsd(p, q):
    p = np.array(p)
    q = np.array(q)
    m = 0.5 * (p + q)
    return 0.5 * (kld(p, m) + kld(q, m))