Hi Blue Jay,
If your sequences are small (<10Kb) and you have a few samples (~100)
the seqinr package from CRAN itself is a very straightforward way of
handling DNA sequences, but if you plan to do more sophisticated things
that can drain your RAM and processor time, I would definitely recommend
the more versatile but complex Biostrings package from Bioconductor.
www.bioconductor.org/packages/2.2/bioc/html/*Biostrings*.html
Happy DNARing!
Keo.
El 23/08/10 05:30 a.m., blue jay escribió:
Hi,
I want to analyse DNA sequence data (mtDNA) in R as in calculate Fst,
Heterozygosity and such summary statistics. Package Adagenet converts the
DNA sequence into retaining only retaining the polymorphic sites and then
calcuates Fst.. but is there any other way to do this? I mean analyse the
DNA sequence as it is.. and calculate the statistics?
Thanks!
[[alternative HTML version deleted]]
______________________________________________
[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.
______________________________________________
[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.
