Как нарисовать гауссовскую графическую модель в R

Question

У меня есть матрица коэффициентов корреляции R и матрица коэффициентов парциальной корреляции Rp, то как я могу нарисовать гауссову графическую модель в R?

Было бы лучше, если бы рекомендовал введение некоторых книг о графической модели Гаусса, действительно, я не знаю, что это такое, но первое, что мне нужно сделать, - это вытащить его. Большое спасибо!

#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) 

Тогда как я мог бы нарисовать графическую модель Гаусса?

Answer 1

Если вы хотите построить соответствующий график, вы можете использовать пакет igraph.

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

Answer 2

Я не знаком с термином «гауссовой графической модели», хотя я чувствую, что должно быть (я буду читать на нем, спасибо).

Но для визуализации (частичных) корреляционных матриц вы можете использовать пакет qgraph, который предназначен для этого. Например:

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

Вычисление частичных корреляций строятся с аргументом:

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

дает тот же результат.