ランダムな文字列を生成して特定の文字列の割合を計算する方法は？

Question

私は初めてのRコースをやっています。演習の1つは、1000塩基のランダムなDNAストリングを作成することです。そのためには、GCパーセント （GC％）を計算します。

私はDNA塩基ベクトルを作成し、シーケンスを作成しようとしました、しかし、結果は

DNA <- c("A","G","T","C") 
seq <- strrep(DNA, 250) 

このnoobのための任意の提案は正しくないのですか？

Answer 1

基本的にはRの非常に基本的な、すなわちsum、sample、pasteの機能をチェックアウトする必要があります。解決策は非常に簡単です：ちょうどそれを一緒に貼り、

DNA <- c("A","G","T","C") 
chain <- sample(DNA, 1000, replace = T) 
GCcontent <- sum(chain %in% c('G','C'))/length(chain) 
# chain as a string 
paste(chain, collapse="")