perl -anle '
BEGIN { $is_tpm = 1 }
if ($. > 1) {
for ($i=1; $i<@F; $i++) {
$sum[$i] += $F[$i]
}
}
END {
for ($i=1; $i<@sum; $i++) {
$is_tpm = 0 if ($sum[$i] < 999000 || $sum[$i] > 1001000);
printf "Col_%d: %.0f\n", $i, $sum[$i]
}
print $is_tpm ? "\n结论: 纯TPM矩阵 ✓" : "\n结论: 非TPM或已TMM标准化 ✗"
}
' genes.TPM.matrix
```
---
## 参数说明
```
-a : 自动分割每行到@F数组(相当于awk的自动分割)
-n : 循环读取每一行,但不自动打印
-l : 自动chomp行尾换行符
-e : 执行后面的Perl代码