我正在使用Matlab中的Cobra Toolbox进行双基因敲除研究,并且生长比率的输出是100乘100的矩阵,称为grRatioDble。我需要找到该矩阵的元素的行和列索引,其<0.001,排除了单基因敲除所必需的行。我有一个我要排除的行索引的列矩阵。有一个简单的方法吗?
(nb。我不能只删除矩阵中不需要的行,因为其余行的行,列索引会发生变化)
谢谢, 萨拉。
答案 0 :(得分:0)
这段代码应该做的工作:
1 /获取所有行/列索引grRatioDble<0.001
:
[row,col]=find(grRatioDble<0.001);
2 /排除不需要的行(例如包含不需要的行的向量为rows2exclude
):
row=row(~ismember(row,rows2exclude));
col=col(~ismember(row,rows2exclude));