我正在使用MATLAB中的Cobra Toolbox进行双基因敲除研究,生长比率的输出是一个100×100的矩阵,称为grRatioDble。我需要找到该矩阵中元素的行和列索引,这些元素<0.001,不包括单基因敲除所必需的行。我有一个要排除的行索引的单列矩阵。有没有简单的方法可以做到这一点?
(注意:我不能只是从矩阵中删除不需要的行,然后其余单元格的行,列索引更改)
这段代码应该完成工作:
-
获取所有行/列索引,其中
grRatioDble<0.001
:[row,col] = find(grRatioDble<0.001);
排除不需要的行 (假设包含不需要的行的向量是
rows2exclude
):row = row(~ismember(row, rows2exclude)); col = col(~ismember(row, rows2exclude));