如何在R语言中从pcalg包的GES拟合结果获取邻接矩阵?
获取pcalg包GES结果中的邻接矩阵
你可以通过以下两种方式从ges.fit$essgraph(一个EssGraph类对象)中提取邻接矩阵:
方法一:类型转换提取
将EssGraph对象直接转换为矩阵类型:adj_matrix <- as(ges.fit$essgraph, "matrix")矩阵中,
0代表变量间无边,1代表无向边,2代表有向边(方向可结合绘图结果对应)。方法二:使用内置函数提取
调用pcalg包专门的邻接矩阵提取函数:adj_matrix <- getAdjacency(ges.fit$essgraph)
提取完成后,可通过print(adj_matrix)查看结果,矩阵的行/列顺序与原始数据gmG8$x的变量顺序一致,可搭配colnames(gmG8$x)匹配对应变量名称。
内容的提问来源于stack exchange,提问作者Jackie
相关产品推荐
相关产品推荐

