代码:
# -*- coding: utf-8 -*-'''Created on 5月15日@author: user@attention: dirichret distribution'''from scipy.stats import dirichletimport matplotlib.pyplot as pltimport numpy as npTOL = 1.48e-8def test_dirichret_distribution():print (dirichlet.pdf([0.6,0.3,0.1],[3,2,1]))print (dirichlet.logpdf([0.6,0.3,0.1],[1,2,3]))if __name__ == '__main__': #scipy.stats库test_dirichret_distribution()#beta分布#pip install dirichlet #/ericsuh/dirichlet
结果:
6.479999999999995-1.7147984280919264