• R语言
  • 基因表达差异比较的问题

一个老问题再拿出来问一下

因为最近想把单样本,多样本比较整理一下,关于这个问题,想再请教一下各位

我有一列基因在两个样本的表达量:

<br />
Gene	S1	S2<br />
ART4	4.92143	2.89385<br />
AZGP1	3.87462	3.86413<br />
C20orf103	3.38678	3.451<br />
C3	3.85766	3.66852<br />
C7orf16	4.01758	4.49015<br />
CDH5	6.68535	3.90784<br />
CPM	3.57492	4.26991<br />
CYP4X1	3.45229	3.59845<br />
DDX3Y	7.3952	7.60647<br />
EDNRB	6.1241	3.27362<br />
EIF1AY	8.79199	8.71892<br />
ELTD1	6.90271	4.9844<br />
GABRA2	4.28023	3.57325<br />
GABRP	6.74015	4.95259<br />
GDF3	3.1609	3.38629<br />
GLIPR1	5.21847	3.50253<br />
GPR50	3.93796	3.66198<br />
HIST1H2BB	4.10468	4.1237<br />
HPSE	3.67252	3.93885</p>
<p>
</p>

我想找出哪些基因在S1和S2之间是有显著差异的

t-test 好像只能给出两组数总的pvalue, 如果我想得到每个基因在s1和s2之间的pvalue, 应该用什么test

如果s1 和s2 有重复,我知道怎么做,但是对于这种没有重复的两个样本比较,除了找fold-change, 统计学上还有什么合适的test吗。。。

Thanks a lot!