utpal <utpalmtbi <at> gmail.com> writes:
> I have a sequence file (fasta format) and want to calculate the rho
> statistics for dinucleotide abundance value on my data.. the code which I
> use is (using seqinr library and current working directory)
>
> seq_info<-read.fasta("gene.txt")
> rho(seq_info[1],2)
> but it yields only the dinucleotides, not their rho values, i.e,
> > rho(seq_info[1],2)
> aa ac ag at ca cc cg ct ga gc gg gt ta tc tg tt
You will likely have better luck with this query on the
Bioconductor mailing list ...
______________________________________________
[email protected] 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.
