Wie erhält man die durchschnittliche Anzahl paarweiser Unterschiede zwischen Populationen?

Ich habe genetische Daten im .structure- und .vcf-Format (und kann andere Formate mit PGDSpider leicht erreichen ). Die interessierende Population ist strukturiert, und ich möchte die durchschnittliche Anzahl paarweiser Unterschiede (über ein bestimmtes Fenster) zwischen einem beliebigen Paar von Populationen berechnen.

Gibt es ein vorhandenes Tool, das dies tun würde?

Antworten (1)

In Bezug auf die andere Bioleiterfrage möchte ich Sie auf diese Pakete hinweisen;

  1. VariantAnnotation : Bietet Funktionen zum Lesen und Bearbeiten von vcf-Dateien.
  2. GenomicRanges : Stellt Funktionen zum Erstellen von Genomintervallen bereit. Bietet auch Überlappungsfunktionen zum Berechnen von Überlappungen zwischen verschiedenen Bereichen.
  3. rtracklayer : Ermöglicht den einfachen Formatwechsel zwischen verschiedenen biologischen Dateiformaten.
  4. AnnotationHub : Stellt viele gängige und ungewöhnliche Anmerkungsdateien im R-Terminal in verschiedenen Formaten wie TxDb, OrgDb usw. usw. zur Verfügung.

Hier ist ein Tutorial, das dem, was Sie fragen, sehr ähnlich ist. Es verwendet einige der oben genannten Pakete. Hauptsächlich AnnotationHub, GenomicRanges und VariantAnnotation.

Wenn Sie das Handbuch dieser Pakete wünschen, sind sie auf den verlinkten Seiten verfügbar.