2017-06-06 81 views
-1

我正在處理兩個數據幀都是20502行,我想要執行t.test,第一行是一個數據幀,第一行是另一個數據幀,兩個數據幀的20502行。T.test對於R中的每一行

PS:我在數據有缺失值(NA)幀

+1

你爲什麼要這樣做?要麼以錯誤的方式使用數據幀,要麼使用不正確的t檢測。 – Consistency

+0

我想查看正常樣本和腫瘤樣本的探針之間的差異甲基化。 –

回答

1

如果我理解你正確地做什麼,這是不使用data.frames的典型方式 - 通常,應該作爲變量而不是行,並且可用的data.frame工具的主機使用這個假設。話雖這麼說,對於data.frames A及以下B ...

A <- as.data.frame(matrix(1:12,nrow=3,ncol=4)) 
A 
## V1 V2 V3 V4 
## 1 1 4 7 10 
## 2 2 5 8 11 
## 3 3 6 9 12 

B <- as.data.frame(matrix(2:13,nrow=3,ncol=4)) 
B 
## V1 V2 V3 V4 
## 1 2 5 8 11 
## 2 3 6 9 12 
## 3 4 7 10 13 

...你想要做的是c(1,4,7,10)c(2,5,8,11)等的於各行的2樣本t檢驗, 一直往下。你可以到達那裏

pvals <- rep(NA, nrow(A)) 
for(i in 1:nrow(A)) pvals[i] <- t.test(A[i,],B[i,])$p.value