以下是我使用该multtest
包查找差异表达基因的矩阵。
ID_REF GSM362168 GSM362169 GSM362170 GSM362171 GSM362172 GSM362173 GSM362174
244901_at 5.171072 5.207896 5.191145 5.067809 5.010239 5.556884 4.879528
244902_at 5.296012 5.460796 5.419633 5.440318 5.234789 7.567894 6.908795
我使用的代码是
cl<-c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,
1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
stat <- mt.teststat(conthyp,classlabel=cl,test="t.equalvar")
conthyp
matic在哪里,cl
是classlabel
. 我收到以下错误
错误:C 堆栈使用量太接近限制
我不确定我在哪里犯了错误。