A for循環與strsplit在R錯誤

我想在這篇文章的前言說我是R和生物信息學編碼的新手，並且我非常感謝來自這個有知識的社區的一些輸入。我在下面發佈的代碼的目標是生成餅圖，顯示BLAST結果中每種蛋白質的氨基酸丰度。我從UniProt上傳了一個csv文件，將其轉換爲矩陣，並寫下了下面的代碼。我不斷收到錯誤：在AAs [i] =表中（strsplit（BLAST_AA_seqs [i]，「」，useBytes = TRUE））：要替換的項目數不是替換長度的倍數。第8欄是包含氨基酸序列的輸出欄。提前致謝！A for循環與strsplit在R錯誤

mydata=read.table("CDPKbeta_BLAST_results.csv",header=TRUE,sep=",") 
mydata=as.matrix(mydata) 

AAs=c() 
BLAST_AA_seqs=c() 
for(i in 1:nrow(mydata)){ 
    print(i) 
    BLAST_AA_seqs[i]=mydata[i,8] 
    AAs[i]=table(strsplit(BLAST_AA_seqs[i],"", useBytes=TRUE)) 
    pie(AAs, col=rainbow(length(AAs)), main="Residue abundance") 
    }

來源

2014-12-02 Courtney Merz

顯示我們的樣本數據。 'head（mydata）'，'dput（mydata）'？和預期的最終結果數據/情節？ – zx8754 2014-12-02 08:53:11

lal<-c() 
lal[1]=table(strsplit("", "", useBytes = T))

空字符串的問題（BLAST_AA_seqs [i]是空字符串）

來源

2014-12-02 09:01:47 ABA

A for循環與strsplit在R錯誤

回答

相關問題