R help - Jan 2012 - rho stat from a fasta sequence file

Hi all, 

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 

I will be grateful if anyone solve this.. I've also attached the sequence
file.. 
Thanks in advance.. 
Utpal 

--
View this message in context:
http://r.789695.n4.nabble.com/rho-stat-from-a-fasta-sequence-file-tp4298621p4298621.html
Sent from the R help mailing list archive at Nabble.com.

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 ...