ANGSD: Analysis of next generation Sequencing Data

Latest tar.gz version is (0.938/0.939 on github), see Change_log for changes, and download it here.

Association

From angsd
Jump to navigation Jump to search

Association

Association can be performed using two approaches. One based on testing differences in allele frequencies between cases and control while the other is based on a generalized linear framework which also allows for quantitative traits and for including additional covariates. Both methods takes the uncertaincy of the genotypes into account.

Case control association using allele frequencies

To test for differences in the allele frequency genotype likelihood need to be provided or estimated. The test is an implimentation of the likelihoods ratio test for differences between cases and controls Kim2011.If alignment files are your input then -GL must be invoked.

-doAsso [int]

1: The test is performed assuming the minor allele is known
3: The test is performed summing over all possible minor alleles

-yBin [file]

a file containing the case control status. 0 being the controls, 1 being the cases and -999 being missing phenotypes. The file should contain a single phenotype entry per line. Example of cases control phenotype file

Score statistic

To perform the test in a generalized linear framework posterior genotype probabilities must be provided or estimated. The approach is published here citation. If alignment files are your input then -doLike, -doMaf, -doPost must be invoked. If input files are genotype likelihoods then -doMaf, -doPost must be invoked. Beagle output files can be used directly.

-doAsso [int]

2: The test is based on a score statistic from a generalized linear framework

-yBin [file]

a file containing the case control status. 0 being the controls, 1 being the cases and -999 being missing phenotypes. The file should contain a single phenotype entry per line. Example of cases control phenotype file

-yQuant [file]

a file containing the phenotype values.-999 being missing phenotypes. The file should contain a single phenotype entry per line. Example of a quantitative phenotype file

-cov [file]

a files containing additional covariates in the analysis. Each lines should contain the additional covariates for a single individuals. Thus the number of lines should match the number of individuals and the number of coloums should match the number of additional covariates. Example of a covariate file

-minHigh [int]

default = 10
This approach needs a certain amount of variability in the genotype probabilities. minHigh filters out sites that does not have at least [int] number of heterozygoes and homogoes genotype with at least 0.9 probability. This filter avoids the scenario where all individuals are heterozygoes with a high probability.

-minCount [int]

default = 10
The minimum expected minor alleles in the sample. This is the frequency multiplied by to times the number of individuals. Performing association on extremely low minor allele frequencies does not make sence.

output

Association

Score statistic (prefix lrt*)

Chromosome Position major minor Frequency N LRT highHe highHo
  • Chromosome

The Chromosome

  • Position

The physical Position

  • major

Often the most common allele. If beagle input is used it might be the minor (see frequency)

  • major

Often the minorallele. If beagle input is used it might be the major (see frequency)

  • Frequency

The frequency estimate. The choice of estimation is determined by the *doMaf option.

  • N

The number of individuals with non-missing data. That is the individuals who have both some sequencing data for the given site and have phenotype data

  • LRT

The likelihood ratio statistic. This statistic is chi square distributed with one degree of freedom. Sites that fails one of the filters are given the value -999.000000

  • highHe

Number of sites with a heterozygoes posterior probability above 0.9

  • highHe

Number of sites with a homozygote posterior probability above 0.9


example:

chr     position        major   minor   frequency       Nind    LRT     highHe  highHo
1       13999950        2       3       0.208311        330     -999.000000     0       20
1       14000023        1       0       0.032045        330     4.429480        20      310
1       14000072        2       3       0.015167        330     0.262688        10      320
1       14000113        0       2       0.015194        330     0.404170        10      320
1       14000202        2       0       0.245345        330     0.562326        90      60
1       14000375        3       1       0.015211        330     0.263012        10      320
1       14000851        3       1       0.015177        330     0.478610        10      320
1       14000873        2       0       0.303030        330     1.032727        100     230
1       14001008        3       1       0.015158        330     3.831819        10      320
1       14001018        3       1       0.300947        330     0.859651        80      240