Estimation of the covariance structure from SNP allele frequencies

01/22/2022
by   Jan van Waaij, et al.
0

We propose two new statistics, V and S, to disentangle the population history of related populations from SNP frequency data. If the populations are related by a tree, we show by theoretical means as well as by simulation that the new statistics are able to identify the root of a tree correctly, in contrast to standard statistics, such as the observed matrix of F2-statistics (distances between pairs of populations). The statistic V is obtained by averaging over all SNPs (similar to standard statistics). Its expectation is the true covariance matrix of the observed population SNP frequencies, offset by a matrix with identical entries. In contrast, the statistic S is put in a Bayesian context and is obtained by averaging over pairs of SNPs, such that each SNP is only used once. It thus makes use of the joint distribution of pairs of SNPs. In addition, we provide a number of novel mathematical results about old and new statistics, and their mutual relationship.

READ FULL TEXT

page 20

page 21

research
06/15/2023

A CLT for the difference of eigenvalue statistics of sample covariance matrices

In the case where the dimension of the data grows at the same rate as th...
research
11/01/2022

Fluctuations of the diagonal entries of a large sample precision matrix

For a given p× n data matrix X_n with i.i.d. centered entries and a popu...
research
07/17/2020

Scalable Methods for Calculating Term Co-Occurrence Frequencies

Search techniques make use of elementary information such as term freque...
research
03/21/2018

High-dimensional covariance matrices in elliptical distributions with application to spherical test

This paper discusses fluctuations of linear spectral statistics of high-...
research
12/11/2019

Graph quilting: graphical model selection from partially observed covariances

We investigate the problem of conditional dependence graph estimation wh...
research
08/05/2021

On rank statistics of PageRank and MarkovRank

An important statistic in analyzing some (finite) network data, called P...
research
10/11/2012

Inferring clonal evolution of tumors from single nucleotide somatic mutations

High-throughput sequencing allows the detection and quantification of fr...

Please sign up or login with your details

Forgot password? Click here to reset