Clojure で遷移行列を操作しようとしています。たとえば、年間の債券格付け推移行列を四半期ごとに変換するには、正方行列の 0.25 乗が必要です。
Python では、scipy からの fractional_matrix_power を次のように取得します。
>>> from scipy.linalg import fractional_matrix_power
>>> a = np.array([[1.0, 3.0], [1.0, 4.0]])
>>> b = fractional_matrix_power(a, 0.5)
>>> b
array([[ 0.75592895, 1.13389342],
[ 0.37796447, 1.88982237]])
IncanterとParallel Coltを検索しても、まだ何も見つかりません。javadocs と google 検索を調べても役に立ちませんでしたが、「マトリックス」と「パワー」は一般的すぎて、探しているものにドリルダウンできない可能性があります。
私は本当にpythonまたはR関数をトランスコードする必要がありますか、それとも私が見逃しているクールなcolt docサイトはありますか?