I realize function write FASTA expects a list with two items, respectively,
description and sequence.
However, just passing a list won't work (please, see code at the bottom of this
message)
I saw there is the helper function CharacterToFASTArecords(x) that presumably
generates the right input data format.
It would b very useful to get some example of CharacterToFASTArecords(x) usage.
The on-line documentation reads:
"For CharacterToFASTArecords, the (possibly named) character vector to be
converted to a list of FASTA records as one returned by readFASTA"
Since I have description and sequnce in separate variables ... I do not know
how to use it.
zz <- file (filname,"w")
write(miRNA.rec, zz, append = FALSE)
write(miRNA.seq,zz, append = TRUE)
#
geneDesc <- paste (">",gene.id, "|",
gene.map[i,"ensembl_transcript_id"], sep="")
geneSeq <- gene.seq[i,"3utr"]
gene.string <- list(desc=geneDesc, seq=geneSeq)
writeFASTA (gene.string, zz)
Thank you in advance for your help.
Maura
tutti i telefonini TIM!
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
