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
あり、cl
はclasslabel
です。次のエラーが表示されます
エラー: C スタックの使用量が限界に近すぎます
どこで間違いを犯しているのかわかりません。