2016-03-03 65 views
2

我有這樣伯爵出現並排除失蹤

aa<-structure(c("AABB", "AABB", NA, "AABB", "AABB", "AABB", "AABB", 
      "AABB", "AABB", "AABB", "AAAA", "AAAA", NA, "AAAA", "AAAA", "AAAA", 
      "AAAA", "AAAA", "AAAA", "AAAA", "BBBB", NA, NA, NA, "AAAA", "AAAA", 
      NA, NA, NA, NA, "AAAA", NA, NA, NA, "AAAA", "BBBB", NA, NA, NA, 
      NA, "AABB", NA, NA, NA, "AABB", "AAAA", NA, NA, NA, NA, "AAAA", 
      "AAAA", "AAAA", "BBBB", "AAAA", "BBBB", "BBBB", "BBBB", "BBBB", 
      "BBBB", "AABB", "AABB", "AABB", "AAAA", "AABB", "AAAA", "AABB", 
      "AAAA", "AAAA", "AAAB", "BBBB", "BBBB", NA, "AABB", "AABB", "AABB", 
      "AABB", "AABB", "AABB", "AABB", "AAAA", "AAAA", NA, "AAAA", "AAAA", 
      "AAAA", "AAAA", "AAAA", "AAAA", "AAAA", "BBBB", "BBBB", NA, "BBBB", 
      "BBBB", "AAAA", "AAAA", "BBBB", "BBBB", "ABBB"), .Dim = c(10L, 10L)) 

的數據集我想在每個計算「A」。我嘗試了2種方法。

str_countstringr包裝。

> apply(aa,2,str_count,"A") 
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
[1,] 2 4 0 4 2 4 2 0 4  0 
[2,] 2 4 1 1 1 4 2 0 4  0 
[3,] 1 1 1 1 1 4 2 1 1  1 
[4,] 2 4 1 1 1 0 4 2 4  0 
[5,] 2 4 4 4 2 4 2 2 4  0 
[6,] 2 4 4 0 4 0 4 2 4  4 
[7,] 2 4 1 1 1 0 2 2 4  4 
[8,] 2 4 1 1 1 0 4 2 4  0 
[9,] 2 4 1 1 1 0 4 2 4  0 
[10,] 2 4 1 1 1 0 3 2 4  1 

這裏我得到1的缺失。但我喜歡有NA

regex

dosage<-function(string,char){ 

    x<-sapply(regmatches(string, gregexpr(char, string)), length) 
    return(x) 
} 

apply(aa,2,dosage,"A") 

     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
[1,] 2 4 0 4 2 4 2 0 4  0 
[2,] 2 4 0 0 0 4 2 0 4  0 
[3,] 0 0 0 0 0 4 2 0 0  0 
[4,] 2 4 0 0 0 0 4 2 4  0 
[5,] 2 4 4 4 2 4 2 2 4  0 
[6,] 2 4 4 0 4 0 4 2 4  4 
[7,] 2 4 0 0 0 0 2 2 4  4 
[8,] 2 4 0 0 0 0 4 2 4  0 
[9,] 2 4 0 0 0 0 4 2 4  0 
[10,] 2 4 0 0 0 0 3 2 4  1 

在這裏我得到0,但我又喜歡有NA

我該怎麼做?

+2

你可以只更換'回報(X)'和'ifelse(is.na(串),NA,X)'和'的功能dosage'將工作,只要你想 – rawr

+1

@rawr是的,這個工作... – Ananta

+0

你的'apply'函數給了我預期的輸出 – Jaap

回答

1

您可以直接將str_count直接格式化爲matrix

library(stringr) 
matrix(str_count(aa, "A"), nrow=nrow(aa), ncol = ncol(aa)) 
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
[1,] 2 4 0 4 2 4 2 0 4  0 
[2,] 2 4 NA NA NA 4 2 0 4  0 
[3,] NA NA NA NA NA 4 2 NA NA NA 
[4,] 2 4 NA NA NA 0 4 2 4  0 
[5,] 2 4 4 4 2 4 2 2 4  0 
[6,] 2 4 4 0 4 0 4 2 4  4 
[7,] 2 4 NA NA NA 0 2 2 4  4 
[8,] 2 4 NA NA NA 0 4 2 4  0 
[9,] 2 4 NA NA NA 0 4 2 4  0 
[10,] 2 4 NA NA NA 0 3 2 4  1 
+0

仍然給我第一個輸出。似乎沒有認識到「NA」缺失。我試過'is.na(aa)',它給了預期的TRUE/FALSE,雖然 – Ananta

+0

@Ananta這對我有用。也許嘗試在新的R會話。 ...其實你的原始'apply'代碼也適用於我。 – Frank

+0

有趣的是,原來的和這個解決方案在Windows中工作,但不是Linux。 – Ananta