2016-12-08 2 views
-1

パッケージmpmiを使用して2つの連続変数間の相互情報を計算しようとします。私はGutHubに置くソースコードの混乱しています:https://github.com/cran/mpmi/blob/master/src/cminjk.f95Rパッケージ "mpmi" 2つの連続変数間の相互情報を計算する

ans = ans + log(s12(i)/(s1(i) * s2(i))) 

end do 

ans = ans/lv + log(dble(lv)) 

S12はp(x,y)のように思えるし、S1、S2はp(x),p(y)のように思えます。 なぜp(x,y)log(s12(i)/(s1(i) * s2(i)))前に、式のためにMIがp(x,y)*log(p(x,y)/(p(x)p(y)))

ある計算に掛けていないと、なぜ和を終えた後ans = ans/lv + log(dble(lv))ありますか?この参照して

+0

私は、これらの質問がパッケージ作成者/管理者でよりうまく対処できると思います。 –

+0

私は保守担当者に連絡しましたが、返事はまだ届いていません。 –

答えて

0

参照「遺伝子 制御ネットワークの再構築のためのペアごとの相互情報の高速計算」を参照してください

0

ケア:「遺伝子制御ネットワークの再構築のためのペアごとの相互情報の高速計算」

この戦略では、ということを前提としてい2つの確率変数のそれぞれは正規分布に従う。両方の変数が正規分布に従っていることを確実に知っている場合を除いて、私は分布の仮定がないパッケージの使用を提案します。私はRパッケージ "minerva"をお勧めします。 MIC(X、Y)$ MICを指定するだけです。

関連する問題