2012-08-06 17 views
7

, bu çıkışı dosya şuna benzer:R paketini kullanarak fasta dosyalarını yazma seqinr? Ben seqinr içinde write.fasta kullandığınızda

başka deyişle
>Sequence name 1 

>Sequence name 2 

>Sequence name 3 
...etc 

Sequence 1 Sequence 2 Sequence 3 ...etc 

, dizi adları dosyasının başında hepsi ve daha sonra diziler çıkış birlikteyiz dosyanın sonunda.

>Sequence name 1 
Sequence 1 
>Sequence name 2 
Sequence 2 
>Sequence name 3 
Sequence 3 
...etc 

write.fasta ile olması mümkün mü:

Ne yapmak istiyorum bu nedir?

+3

bir [tekrarlanabilir örnek] yayınlamak misiniz (http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)? Örneğin, 'write.fasta' olarak adlandırmak için kullandığınız kodu yayınlayın ve ona ne geçtiğini göstermek için 'dput' kullanın. –

cevap

2

Aslında, her zaman doğru yolu buluyorum ve bana hiç benzemeyen bir problem yaşamadım. Bunu dene. Bu metin aşağıda

Kopya: "test.fasta" olarak kaydet

>seq1 
agctgtagtc 
>seq2 
agtctctctt 
>seq3 
atgtataaaa 

. Sonra R siz "write.my.dna.fasta" açarsanız

my.dna<-read.fasta("test.fasta") 
write.fasta(sequences=my.dna,names=names(my.dna),file.out="write.my.dna.fasta") 

Aşağıda sonra aşağıdaki elde edecek yapın:

>seq1 
agctgtagtc 
>seq2 
agtctctctt 
>seq3 
atgtataaaa 
14

Ben benzer bir sorun yaşıyordu. Yaptığım şey, dizileri içeren vektörü bir listeye dönüştürmek ve iyi çalıştı. örneğin

, write.fasta(as.list(seq),names,file)

İlgili konular