1

文字のベクトルがある場合：

> all <- letters
> all
[1] "a" "b" "c" "d" "e" "f" "g" "h" "i" "j" "k" "l" "m" "n" "o" "p" "q" "r" "s" "t" "u" "v" "w" "x" "y" "z"

次に、文字から参照サンプルを次のように定義します。

> refSample <- c("j","l","m","s")

要素間の間隔が2（1番目から2番目）、1（2番目から3番目）、および6（3番目から4番目）である場合、要素間に同じ非ラップアラウンド間隔を持つn個のサンプルを選択するにはどうすればよいですか？？たとえば、とは有効なサンプルですが、そうで はありません。前者は3番目の要素と最後の要素の間で（6ではなく）7のインデックスの違いがありますが、後者は正しい間隔ですがラップアラウンドします。allrefSample"a","c","d","j""q" "s" "t" "z""a","c","d","k""r" "t" "u" "a"

次に、これをパラメータ化して、使用するものが何であれ、refSampleその間隔をテンプレートとして使用できるようにするにはどうすればよいですか？

4

1 に答える 1

3

これが簡単な方法です-

all <- letters                                                                                                                                                                                                                                                                
refSample <- c("j","l","m","s")                                                                                                                                                                                                                                               


pick_matches <- function(n, ref, full) {                                                                                                                                                                                                                                      
  iref <- match(ref,full)                                                                                                                                                                                                                                                     
  spaces <- diff(iref)                                                                                                                                                                                                                                                        
  tot_space <- sum(spaces)                                                                                                                                                                                                                                                    
  max_start <- length(full)  - tot_space                                                                                                                                                                                                                                      
  starts <- sample(1:max_start, n, replace = TRUE)                                                                                                                                                                                                                            
  return( sapply( starts, function(s) full[ cumsum(c(s, spaces)) ] ) )                                                                                                                                                                                                        
}                                                                                                                                                                                                                                                                             

> set.seed(1)                                                                                                                                                                                                                                                                
> pick_matches(5, refSample, all) # each COLUMN is a desired sample vector                                                                                                                                                                                                                                         
      [,1] [,2] [,3] [,4] [,5]                                                                                                                                                                                                                                                
 [1,] "e"  "g"  "j"  "p"  "d"                                                                                                                                                                                                                                                 
 [2,] "g"  "i"  "l"  "r"  "f"                                                                                                                                                                                                                                                 
 [3,] "h"  "j"  "m"  "s"  "g"                                                                                                                                                                                                                                                 
 [4,] "n"  "p"  "s"  "y"  "m"
于 2012-05-03T20:33:37.070 に答える