Skip to content

Commit

Permalink
brief introduction added to README
Browse files Browse the repository at this point in the history
  • Loading branch information
NBMueller committed Feb 13, 2023
1 parent a55f196 commit 8a21813
Showing 1 changed file with 9 additions and 4 deletions.
13 changes: 9 additions & 4 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,7 +1,12 @@
# Single-Cell SOMatic MolEculaR Clock testing

This repository contains scripts for running
- the **Poisson Tree (PT) Test**
This repo contains a **Poisson Tree (PT) Test** for the existence of a somatic clock in single-cell phylogenies.
In short, it tests if different cell lineages evolve at a similar rate, accumulating mutations according to a molecular clock.
As input the test requires a mutation matrix, a phylogeny of contemporaneously sampled cells, and error rates.


This repo contains scripts for running
- the **PT Test**
- and, in a subfolder (AnalysisPipelines), scripts for
- the processing of real scDNA-seq data
- the analysis of real scDNA-seq data
Expand Down Expand Up @@ -39,7 +44,7 @@ The **PT test** requires two input files:
- `-o <str>`, Output file. Default = <VCF_FILE>.poissonTree_LRT.tsv.
- `-excl <str>`, Regex pattern for samples/cells to exclude. Default = none.
- `-incl <str>`, Regex pattern for samples/cells to include. If set, only these samples/cells are included. Default = all cells.
- `-w <list of int>`, Maximum weight values. Defaut = 100, 200, ..., 1000'.
- `-w <list of int>`, Maximum weight values. Default = 100, 200, ..., 1000'.
- `-FN <float>`, Estimated FN rate (for CellPhy and infSCITE: inferred from .log/stdout file).
- `-FP <float>`, Estimated FP rate (for CellPhy and infSCITE: inferred from .log/stdout file).

Expand All @@ -55,4 +60,4 @@ python run_PT_test.py example_data example_data/data_simulated_noclock.vcf.gz ex

The former data is simulated under a molecular clock, the later with a deviation from the clock (evolutionary rate amplified by 5x in a subtree)
> ## Note
> FN and FP rate are inferred from the `.raxml.log` file
> FN and FP rate are inferred from the `.raxml.log` file

0 comments on commit 8a21813

Please sign in to comment.