Skip to content

vcf input with population info #651

@Ida921

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?

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