# The 1000 Genomes Project Data

The outlier detection method (Lindskou, Svante Eriksen, and Tvedebrink 2019) arose from a problem in the forensic science community where it is of great interest to make statements about the geographical origin of a DNA sample. This is in general a very complicated matter. More over, when DNA markers are in linkage disequilibrium things get even more complicated. The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations (The 1000 Genomes Project Consortium 2015). In the molic package we include the final data from the project which includes $$2504$$ DNA profiles coming from five different continental regions; Europe (EUR), Africa (AFR), America (AMR), East Asia (EAS), South Asia (SAS). Each DNA profile is the compound of two rows in the data tgp_dat, one for each chromosome copy. This makes the outlier procedure slightly more complicated since we must fit a model to each copy and aggregate the information (but we shall see in a moment that it is not hard to do using the molic package).

The data includes $$276$$ SNP markers grouped in $$97$$ microhaplotypes; in other words, the SNPs within a microhaplotype are so close that they cannot be assumed to be in linkage disequilibrium and we must take into account their mutual dependencies. We do this with the use of decomposable graphical models. See the outlier_intro for more details on the model. The microhaplotype grouping is provided is the list tgp_haps.

## Setup and Data

We peak at data:

## # A tibble: 1,006 x 304
##    sample_name pop_meta rs28529526 rs10134526 rs260694 rs11123719 rs11691107
##    <chr>       <chr>    <chr>      <chr>      <chr>    <chr>      <chr>
##  1 HG00096a    EUR      G          C          T        T          C
##  2 HG00097a    EUR      G          T          T        C          C
##  3 HG00099a    EUR      G          C          T        T          C
##  4 HG00100a    EUR      G          C          T        T          C
##  5 HG00101a    EUR      G          C          T        T          C
##  6 HG00102a    EUR      G          C          T        T          C
##  7 HG00103a    EUR      G          C          T        C          C
##  8 HG00105a    EUR      G          T          T        C          C
##  9 HG00106a    EUR      G          C          T        C          C
## 10 HG00107a    EUR      G          C          T        T          C
## # … with 996 more rows

We now fit the interaction graphs for each chromosome copy. We use fit_components since we know, that SNPs between haplotypes are independent. This also makes the procedure much faster.

## Example 1 - Testing which observations within data are outliers

Fitting the model:

Model information:

##
##  --------------------------------
##   Simulations: 500
##   Variables: 302
##   Observations: 503
##   Estimated mean: 488.65
##   Estimated variance: 319.11
##  --------------------------------
##   Critical value: 516.5493
##   Alpha: 0.05
##   <mixed_outlier, outlier, outlier_model, list>
##  --------------------------------

The distribution of the test statistic:

Locating the outliers:

## # A tibble: 19 x 302
##    rs28529526 rs10134526 rs260694 rs11123719 rs11691107 rs12469721 rs3101043
##    <chr>      <chr>      <chr>    <chr>      <chr>      <chr>      <chr>
##  1 G          C          T        T          C          A          C
##  2 G          C          T        C          C          T          T
##  3 G          C          T        T          C          T          T
##  4 G          C          T        C          C          A          T
##  5 G          T          T        T          C          A          T
##  6 G          C          T        C          C          T          T
##  7 G          T          T        T          C          A          C
##  8 G          C          T        C          C          A          C
##  9 G          C          T        T          C          A          C
## 10 G          C          T        C          C          T          C
## 11 G          C          T        T          C          T          C
## 12 G          C          T        C          C          A          T
## 13 G          C          T        T          C          T          T
## 14 G          C          T        T          T          T          T
## 15 G          C          T        T          C          A          C
## 16 G          C          T        T          C          T          T
## 17 G          C          T        T          C          T          C
## 18 G          T          T        T          C          T          C
## 19 G          T          T        T          C          A          T

Retrieving the deviance test statistic for individual observations:

## [1] 537.8832
## [1] 496.7722
## [1] 0.004
## [1] 0.328

## Example 2 - Testing if a new observation is an outlier

Testing if an American is an outlier in Europe:

##
##  --------------------------------
##   Simulations: 500
##   Variables: 302
##   Observations: 504
##   Estimated mean: 482.83
##   Estimated variance: 300.25
##  --------------------------------
##   Critical value: 512.2617
##   Deviance: 508.0878
##   P-value: 0.078
##   Alpha: 0.05
##   <mixed_outlier, novelty, outlier_model, list>
##  --------------------------------

Thus, z does not deviate significantly from the Europe database. The red area depicts the $$0.05$$ significance level and the dotted line is the observed deviance of z (508.0877826). This is not a complete surprise, since Americans and Europeans have many similarities.

# References

Lindskou, Mads, Poul Svante Eriksen, and Torben Tvedebrink. 2019. “Outlier Detection in Contingency Tables Using Decomposable Graphical Models.” Scandinavian Journal of Statistics. https://doi.org/10.1111/sjos.12407.

The 1000 Genomes Project Consortium. 2015. “A Global Reference for Human Genetic Variation.” Nature 526 (7571): 68. https://doi.org/10.1038/nature15393.