Skip to content

vcf input with population info #651

Description

@Ida921

Q1: I have a VCF file containing 518 samples, which belong to 25 populations. How can I use angsd to specify the samples to be included for analysis when the input is a VCF file? The -b parameter seems to be applicable only for BAM files. Since I need to calculate FST between these populations, I'll have to compute the SFS for each population separately.
Q2 : I tried splitting the total VCF file into subsets according to population information and calculated the corresponding SFS for each population using the subset VCF files. However, some population VCF files produced the following error during the calculation: [ERROR](safreader.cpp:306)ret==band[1]*sizeof(T) && band[0]+band[1]<=saf->nChr+1. How can I figure this out?

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions