This repository contains all computational analyses used in my 2020 paper “The evolutionary history of Neandertal and Denisovan Y chromosomes”.
You can find here:
- my complete data processing and data-generation pipeline -
processing of BAM files, genotype calling, simulations
- scripts used by the
src/) and R functions used for analyses (
- R Jupyter notebooks used to generate every result in the paper
- cached results in the rds format - useful for re-producing all
tables and figures in the paper without having to first re-process
sequencing data or run SLiM simulations (
data/rds/); for full data access see here.
Note that the repository is structured as an R package, with the
exception that it contains not only the required R package components
(R code in
DESCRIPTION file etc.), but also the data and
analyses themselves. Therefore, after cloning this repository, you can
devtools::load_all(".") from the root of the repository
in an R console and immediately start exploring the data yourself.
All software dependencies used in the project are described in detail in this notebook.
In addition to cached results in
data/rds in this repository, all
sequencing data (BAM/VCF/FASTA files), raw simulation outputs, as well
as the sequences of capture probes can be found at: