如果拿到的表达谱不是原始的read counts数据而是TPM值,就不能用包(比如DESeq2)来进行差异表达基因分析。我们可以手动用wilcox.test 函数手动进行分析。我的数据为log2TPM的表达矩阵 ExM# log2TPM 表达矩阵s1# 属于类型1(如 tumor)的所有样本IDs2# 属于类型2(如 normal)的所有样本IDcat("wilcox.test\n")pvalu...