LeBlanc_Palakurthi#3
Conversation
araiho
left a comment
araiho
left a comment
There was a problem hiding this comment.
looks good! Next step would be to simplify with for loops. 4/4
|
|
||
| write.fasta(sequences=motif1, names=motif1Names,file.out="motif1.fasta") | ||
| write.fasta(sequences=motif2, names=motif2Names,file.out="motif2.fasta") | ||
| write.fasta(sequences=NoMotif, names=NoMotifNames,file.out="NoMotif.fasta") No newline at end of file |
There was a problem hiding this comment.
Interesting approach. Very creative! Could also have:
library(stringr)
inFile=scan("motifsort.fasta",sep="\n",what=character())
motif1=file("motif1.fasta",open="w")
motif2=file("motif2.fasta",open="w")
other=file("nomotif.fasta",open="w")
for(i in 1:length(inFile)){
if(str_detect(inFile[i],">")){
seqID=inFile[i]
}else{
if(str_detect(inFile[i],"AKKPRVZE")){
writeLines(seqID,motif1)
writeLines(inFile[i],motif1)
}else if(str_detect(inFile[i],"AAQWWRNYGG")){
writeLines(seqID,motif2)
writeLines(inFile[i],motif2)
}else{
writeLines(seqID,other)
writeLines(inFile[i],other)
}
}
}
close(motif1)
close(motif2)
close(other)
3 participants
No description provided.