Skip to content

Commit

Permalink
mc
Browse files Browse the repository at this point in the history
  • Loading branch information
jovoni committed Dec 11, 2024
1 parent 0868836 commit 2aff22c
Show file tree
Hide file tree
Showing 12 changed files with 323 additions and 298 deletions.
2 changes: 1 addition & 1 deletion README.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -19,7 +19,7 @@ knitr::opts_chunk$set(

<!-- badges: end -->

# lineaGT <a href="caravagnalab.github.io/lineaGT"><img src="man/figures/tickTack.png" align="right" height="139"/></a>
# tickTack <a href="caravagnalab.github.io/tickTack"><img src="man/figures/logo.png" align="right" height="120" alt="" /></a>

tickTack is a package to infer kinetic parameters of an evolving population whose size is observed at discrete intervals of time. The tool is able to

Expand Down
2 changes: 1 addition & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,7 @@

<!-- badges: end -->

# lineaGT <a href="caravagnalab.github.io/lineaGT"><img src="man/figures/tickTack.png" align="right" height="139"/></a>
# tickTack <a href="caravagnalab.github.io/tickTack"><img src="man/figures/logo.png" align="right" height="120" alt="" /></a>

tickTack is a package to infer kinetic parameters of an evolving
population whose size is observed at discrete intervals of time. The
Expand Down
Binary file modified man/figures/logo.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file removed man/figures/tickTack.png
Binary file not shown.
Binary file modified pkgdown/favicon/apple-touch-icon.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified pkgdown/favicon/favicon-96x96.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified pkgdown/favicon/favicon.ico
Binary file not shown.
28 changes: 27 additions & 1 deletion pkgdown/favicon/favicon.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified pkgdown/favicon/web-app-manifest-192x192.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified pkgdown/favicon/web-app-manifest-512x512.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
18 changes: 10 additions & 8 deletions vignettes/a3_Understanding_single_segment_timing.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -23,11 +23,12 @@ The timing results are obtained in the following steps:
The clonal peaks are calculated using the `get_clonal_peaks` function. The function computes the theoretical VAF peaks based on the karyotype and tumor purity.
They are computed using the following equation:

\[
$$
\text{Peak}_i = \frac{M_i \cdot P}{(N_{\text{tot}} \cdot P) + 2 \cdot (1 - P)}
\]
$$

Where:

- \( M_i \): The multiplicity of the allele of interest (e.g., the major allele or a single copy for heterozygosity).
- \( P \): The tumor purity (a value between 0 and 1).
- \( N_{\text{tot}} \): The total number of copies in the karyotype (e.g., sum of the major and minor alleles).
Expand All @@ -39,16 +40,17 @@ Where:
- The **denominator**, \( (N_{\text{tot}} \cdot P) + 2 \cdot (1 - P) \), normalizes the peak to account for contributions from both tumor and normal diploid cells in the sample.

#### Example Calculation
For a karyotype \( k = "2:1" \) (major allele = 2, minor allele = 1) and tumor purity \( P = 0.4 \):
For a karyotype $k = \text{2:1}$ (major allele = 2, minor allele = 1) and tumor purity $P = 0.4$:

1. \( N_{\text{tot}} = 2 + 1 = 3 \).
2. Compute peaks for \( M_i = 1 \) (single allele) and \( M_i = 2 \) (major allele):
\[
2. Compute peaks for \( M_i = 1 \) (single allele) and \( M_i = 2 \) (major allele)

$$
\text{Peak}_1 = \frac{1 \cdot 0.4}{(3 \cdot 0.4) + 2 \cdot (1 - 0.4)} = \frac{0.4}{1.2 + 1.2} = 0.1667
\]
\[
$$
$$
\text{Peak}_2 = \frac{2 \cdot 0.4}{(3 \cdot 0.4) + 2 \cdot (1 - 0.4)} = \frac{0.8}{1.2 + 1.2} = 0.3333
\]
$$

### Step 2: Filtering Mutations

Expand Down
Loading

0 comments on commit 2aff22c

Please sign in to comment.