查找矩阵中元素的行索引 <0.001, excluding certain rows

Question

我正在 MATLAB 中使用 Cobra Toolbox 进行双基因敲除研究，生长率的输出是一个 100 x 100 的矩阵，称为

grRatioDble

。我需要找到该矩阵元素的行索引和列索引，这些元素是 <0.001, excluding the rows which were essential on single gene knockout. I have a one-column matrix of the row indexes that I want to exclude. Is there an easy way to do this?

（注意：我不能只从矩阵中删除不需要的行，因为其余单元格的行、列索引会发生变化）

Answer 1

这段代码应该可以完成这项工作：

1. 获取所有行/列索引，其中

   grRatioDble<0.001

   :

   [row,col] = find(grRatioDble<0.001);
   

2. 排除不需要的行（假设包含不需要的行的向量是

   rows2exclude

   ）：

   row = row(~ismember(row, rows2exclude));
   col = col(~ismember(row, rows2exclude));