2012-04-19 8 views
2

私は相関係数行列Rと部分相関係数行列Rpを得ました。次に、Rでガウスモデルをどのように描くことができますか?Gaussianグラフィカルモデルを描く方法R

Gaussianグラフィカルモデルについていくつかの本を紹介することを推奨しても、実際には何か分かりませんが、最初にやる必要があるのはそれを引き出すことです。どうもありがとう!

#the Correlation coefficient matrix 
R=c(1,0.55,0.55,0.41,0.39,0.55,1,0.61,0.49,0.44,0.55,0.61,1,0.71, 
0.66,0.41,0.49,0.71,1,0.61,0.39,0.44,0.66,0.61,1) 
dim(R)=c(5,5) 

#the Partial correlation coefficient matrix 
library("corpcor") 
Rp=cor2pcor(R) 

次に、ガウスのグラフィカルモデルをどのように描画できますか?

答えて

3

対応するグラフをプロットする場合は、igraphパッケージを使用できます。

library(igraph) 
g <- graph.adjacency(abs(Rp)>.1, mode="undirected", diag=FALSE) 
plot(g, layout=layout.fruchterman.reingold) 
+0

なぜabs(Rp)>?1? – PepsiCo

+0

そうしないと、完全なグラフが常に得られます。これはむしろ情報価値がありません。あるいは、Sachaの答えのように、相関が低い場合にはより暗い色を使用することもできます。 –

2

「Gaussianグラフィカルモデル」という用語には慣れていませんが、私はそうすべきだと感じています(感謝します)。

しかし、(部分的な)相関行列を視覚化するには、これを行うように設計されたqgraphパッケージを使用することができます。たとえば:偏相関を計算

library("qgraph") 
qgraph(round(Rp,5),edge.labels=TRUE) 

enter image description here

は引数で構築されています:

qgraph(round(R,5),edge.labels=TRUE,graph="concentration") 

は同じ結果を与えます。

+0

私はパッケージ "qgraph"を取得できませんでした:library.dynam(lib、package、package.lib)のエラー: DLL'plyrが必要です – PepsiCo

+0

CRAN経由でインストールすると 'plyr'も自動的にインストールされます。 'install.packages(" qgraph "、dep = TRUE)'を実行します。あなたのRバージョンが最新(2.15)であることを確認してください。 –