On 10 January 2013 01:04, Yao He <yao.h.1...@gmail.com> wrote: > In fact I want to calculate the gene frequency of each SNP.
Why don't you use bioconductor for your analysis instead of trying to develop by your own? For example: http://www.bioconductor.org/help/course-materials/2008/MGED08/BiostringsMGED2008.pdf For example alphabetFrequency function in the ShortRead package. I am sure you can handle I/O somehow to interface with available tools in bioconductor. ______________________________________________ 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.
