美文网首页
R 语言 | 计算可达矩阵

R 语言 | 计算可达矩阵

作者: Bioconductor | 来源:发表于2017-04-30 19:09 被阅读102次

在微博有位朋友问我可达矩阵的计算,于是发了点时间用R语言写出来了。

问题如下:


问题

计算过程:


计算过程

注意:是矩阵的乘法。

代码如下:

A <- matrix(c(0,0,0,0,0,0,0,1,0,0,1,1,0,0,0,0,0,0,0,0,0,0,1,1,0),nrow=5)
A1 <- A+diag(5)

repeat{
  B <- A1 * A1;
  if(identical(B,A1)) {print(B);break} else A1 <- B
}

结果:

     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    1    0    0
[2,]    0    1    1    0    0
[3,]    0    1    1    0    1
[4,]    0    0    0    1    1
[5,]    0    0    0    0    1

重点在repeat的使用。

相关文章

网友评论

      本文标题:R 语言 | 计算可达矩阵

      本文链接:https://www.haomeiwen.com/subject/gwxjbttx.html