我知道如何生成DAG.

library(pcalg)
dag <- randDAG(10, 1, "er") 

但现在我需要得到上面DAG的邻接矩阵.我想确实存在这样的功能,但我不知道它是什么.

聊天的人给了我错误的结果:

# Get the adjacency matrix from the DAG
adjacency_matrix <- amat(dag)

# Print or inspect the adjacency matrix
print(adjacency_matrix)

但不幸的是,amat(Dag)中的错误:找不到函数"amat"

推荐答案

你可以用showAmat

library(pcalg)

dag <- randDAG(10, 1, "er") 
amat <- showAmat(dag)
amat
#>    1         2 3 4         5 6         7         8 9 10
#> 1  0 0.0000000 0 0 0.4144231 0 0.0000000 0.0000000 0  0
#> 2  0 0.0000000 0 0 0.0000000 0 0.0000000 2.0000000 0  0
#> 3  0 0.0000000 0 0 0.0000000 0 0.7542406 0.7771955 0  0
#> 4  0 0.0000000 0 0 0.0000000 0 0.0000000 0.0000000 0  0
#> 5  2 0.0000000 0 0 0.0000000 0 0.0000000 0.0000000 0  0
#> 6  0 0.0000000 0 0 0.0000000 0 0.0000000 0.0000000 0  0
#> 7  0 0.0000000 2 0 0.0000000 0 0.0000000 0.0000000 0  0
#> 8  0 0.3098213 2 0 0.0000000 0 0.0000000 0.0000000 0  0
#> 9  0 0.0000000 0 0 0.0000000 0 0.0000000 0.0000000 0  0
#> 10 0 0.0000000 0 0 0.0000000 0 0.0000000 0.0000000 0  0

R相关问答推荐

重复组的运行计数

如何以编程方式将X轴勾号上的希腊符号合并到R图中?

使用ggplot将平滑线添加到条形图

根据收件箱中的特定值提取列名

如何在ggplot 2线性图的每个方面显示每个组的误差条?

更改Heatmap Annotation对象的名称

从gtsummary包中使用tBL_strata()和tBL_summary()时删除变量标签

如何求解arg必须为NULL或deSolve包的ode函数中的字符向量错误

如何利用模型函数在格图中添加双曲/指数曲线

为什么当用osmdata映射R时会得到相邻状态?

如何使用列表中多个列表中的第一条记录创建数据框

从外部文件读取多个值作为字符向量

有效识别长载体中的高/低命中

绘制采样开始和采样结束之间的事件

从多个线性回归模型中提取系数

在保留列表元素属性的同时替换列表元素

在列表中排列R数据框中的列顺序

在纵向数据集中创建新行

如何在R中创建条形图,使条形图在y轴上围绕0.5而不是0构建条形图?

Data.table::Shift type=允许扩展数据(&Q;LAG&Q;)