[R program]ベイジアンネット

良くわからんがベイジアンネットをやってみる。

library(G1DBN)

data(arth800line)
S1 <- DBNScoreStep1(arth800line[,1:10], method="ls")
G1 <- BuildEdges(S1$S1ls, threshold=0.05, dec=F)
Step1InferredNet <- BuildNetwork(G1, 1:10)


attach(Step1InferredNet)

a<- gplot(abs(AdjMatrix), vertex.cex=1.5, diag=TRUE, displaylabel=TRUE, usecurv=TRUE,boxed.label=FALSE, main="Inferred network - Step 1")

detach(Step1InferredNet)

#coor=でノードの位置を固定する

S1 <- DBNScoreStep1(arth800line[,1:10], method="ls")
G1 <- BuildEdges(S1$S1ls, threshold=0.01, dec=F)
Step1InferredNet <- BuildNetwork(G1, 1:10)

attach(Step1InferredNet)

gplot(abs(AdjMatrix), vertex.cex=1.5, diag=TRUE, displaylabel=TRUE, coor=a, usecurv=TRUE,boxed.label=FALSE, main="Inferred network - Step 1")

detach(Step1InferredNet)

結果はこんな感じ↓

自分に矢印がいくのはどういうときなんだろう？

とりあえず使ってみたものの、使い勝手や解釈は他のパッケージの方が良いかも。