我正在嘗試在 R 中創建一個簡單的腳本來模擬環境對基因的壓力，所以我有以下 ATCG 隨機生成器：
list <- c("A","T","C","G")
samp <- sample(list, 1000, replace = T)
display <- function(a,b,c) {
result <- sample(list, 1000, replace = T)
return(result)
}
a <- display(a)
b <- display(b)
c <- display(c)
代碼作業正常，但在概念上是錯誤的，因為密碼子需要作為一組字符而不是單個字符加入。所以我的問題是，如何加入這個元素 (1:4,5:8,...107:110) ？
我知道 R 中有很多關于遺傳學的軟體包，但我試圖讓它非常簡單，只是作為我遺傳學研究小組的練習。
提前致謝！
uj5u.com熱心網友回復：
也許這就是你的想法：
dat <- setNames( data.frame(
replicate(3, replicate( 1000,
paste0(sample(c("A","C","T","G"),4, replace=T),collapse="") ) ) ),
c("a","b","c") )
結果
dat
a b c
1 CATG CTAG CGAT
2 AGCT GACT TGAC
3 TCAG TAGC TCAG
4 CGAT GACT GTAC
5 TCGA TAGC CTGA
...etc
你應該能夠做的比較dat$a，dat$b和dat$c
轉載請註明出處，本文鏈接：https://www.uj5u.com/gongcheng/356533.html
下一篇：按組匯總值并以其他兩列為條件