2013-05-07 72 views
2

我有一個數據幀(8個變量451個OBS),其具有兩列(6 & 7)如下所示:選擇從兩列較高的值,R

Major  Minor 
    C:726  T:2 
    A:687  G:41 
    T:3  C:725 

我想創建一個總結這一點的專欄。要做到這一點,我不關心在每個單元的字母,但我想更大的數量保持,它在任何行即我希望它看起來是這樣的:

Summary_column 
    726 
    687 
    725 

沒有必要,但對那些想知道我在做什麼的人來說,這是VCFtools程序的輸出結果;它具有計數VCF中等位基因的計數功能,但有時當它顯然更常見時,它會將等位基因命名爲「輕微」。

感謝您的幫助!

回答

3

我會做這樣的事情:

extract <- function(v) { 
    gsub("^.*:", "", v) 
} 
within(d, Summary_column <- pmax(extract(Major), extract(Minor))) 

其中給出:

Major Minor Summary_column 
1 C:726 T:2   726 
2 A:687 G:41   687 
3 T:3 C:725   725 
+0

這是太酷了。 – cianius 2013-05-07 10:24:13

+0

你會如何編輯以獲得較小的數字? – cianius 2013-05-07 10:34:52

+1

只需用'pmin'替換'pmax'即可。 – juba 2013-05-07 10:35:51