2
$\begingroup$

I have genetic data in .structure and .vcf format (and can easily reach other formats with PGDSpider). The population of interest is structured and I would like to calculate the average number of pairwise differences (over a given window) between any pair of populations.

Is there an existing tool that would do it?

Improve this question
$\endgroup$

1 Answer 1

Reset to default
2
$\begingroup$

With relation to the other bioconductor question I would point you towards these packages;

  1. VariantAnnotation: Provides functionality for reading and manipulating vcf files.
  2. GenomicRanges: Provides functionality for creating genome intervals. Also provides overlap functions for calculating overlaps between various ranges.
  3. rtracklayer: Allows you to easily change formats between different biological file formats.
  4. AnnotationHub: Makes many common and uncommon annotation files available in the R terminal in various formats such as TxDb, OrgDb etc. etc.

Here's a tutorial very similar to what you are asking for. It utilises some of the packages mentioned above. Principally, AnnotationHub, GenomicRanges and VariantAnnotation.

If you want the manual of these packages, they are available on the linked pages.

Improve this answer
$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .