以下是我用于使用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堆栈使用率太接近极限
我不确定我在哪里犯错误。