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 http://r.789695.n4.nabble.com/file/n4291676/gene.txt gene.txt -- View this message in context: http://r.789695.n4.nabble.com/rho-statistics-for-dinucleotide-abundance-from-a-sequence-file-tp4291676p4291676.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.