Displaying 1 result from an estimated 1 matches for "dinucleotide".
Did you mean:
dinucleotides
2012 Jan 16
1
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...