偏相関・無向グラフメモ

偏相関係数にもとづいて無向グラフを手っ取り早く描く方法。

使うパッケージ

corpcor

qgraph

library(corpcor)

library(qgraph)

dat=read.table("clipboard") #読み込み，データフレーム形式

cor(dat) #相関係数行列の算出

cor2pcor(cor(dat)) #（相関係数行列に基づく）偏相関係数行列の算出

qgraph(cor2pcor(cor(dat)),edge.labels=T)

これでどん。

qgraph(cor2pcor(cor(dat)),edge.labels=T,minimum=.2,edge.color="black",filetype="png")

オプションなどはこんな感じで。

うむ。

別に偏相関係数行列でなくてもいい。qgraphは色々サポートしているらしい。

相関係数なら，

qgraph(cor(dat),edge.labels=T,minimum=.2,edge.color="black",filetype="png")

うん。