Pythonで正方行列の擬似行列式を得る方法

Question

私は素朴なベイズ分類子のために計算している特異なテストに失敗した行列を持っています。私は式のln(det(sigma))部分を処理しています。

共分散行列が特異である場合、私は擬似行列式を見つけなければなりません。どうやってやるの？

Answer 1

擬似逆行列を計算するためにnumpyのpinvを使用して、それを使って行列式を計算できますか？ http://www.sosmath.com/matrix/inverse/inverse.html

Answer 2

1. 最初に行列の固有値を

 

    eig_values = np.linalg.eig(covarianceMatrix) 
 

を計算

、非ゼロの固有値の積を計算する（これは擬似に等しいです行列の決定値）、

 

    pseudo_determinent = np.product(eig_values[eig_values > 1e-12])