在R中设置循环

时间:2014-02-19 17:58:34

标签: r bioinformatics

我已经在下面的帖子中讨论过类似的问题

How to set a for -loop in R

每个文件内容如下:

FILE_1.FASTA

>>TTBK2_Hsap ,(CK1/TTBK)
MSGGGEQLDILSVGILVKERWKVLRKIGGGGFGEIYDALDMLTRENVALKVESAQQPKQVLKMEVAVLKKLQGKDHVCRFIGCGRNDRFNYVVMQLQGRNLADLRRSQSRGTFT

FILE_2.FASTA

>>TTBK2_Hsap ,(CK1/TTBK)
MSGGGEQLDILSVGILVKERWKVLRKIGGGGFGEIYDALDMLTRENVALKVESAQQPKQVLKMEVAVLKKLQGKDHVCRFIGCGRNDRFNYVVMQLQGRNLADLRRSQSRGTFT

然而,R中还有另一个包,其工作原理如下:

extractAPAAC(x, props = c("Hydrophobicity", "Hydrophilicity"), lambda = 30,
w = 0.05, customprops = NULL)

我尝试创建一个函数来为多个文件序列运行它,程序看起来像这样

read_and_extract <- function(fasta) {
  seq <- readFASTA(fasta)[[1]]
  return(extractAPAAC(seq, props = c("Hydrophobicity", "Hydrophilicity"), lambda = 30,
             w = 0.05, customprops = NULL))

}

setwd("H:\\CC")
fasta_files <- dir(pattern = "[.]fasta$")
aa_comp <- vapply(fasta_files, read_and_extract, rep(pi, 80))
write.csv(aa_comp, file = "C:\\Users\\PAAC.csv")

此程序显示错误

错误:意外','在“w = 0.05,”

但我已经给出w = 0.05的默认值,有人能告诉我实际问题在哪里吗?

0 个答案:

没有答案