Skip to content

Commit

Permalink
Reorder report subsections about temporal signal
Browse files Browse the repository at this point in the history
  • Loading branch information
ahmig committed Nov 8, 2023
1 parent 31fb460 commit 7a6e665
Showing 1 changed file with 19 additions and 19 deletions.
38 changes: 19 additions & 19 deletions template.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -117,10 +117,10 @@ using [Freyja](https://github.com/andersen-lab/Freyja).
![Estimated lineage admixture of each sample.
Samples in the X-axis are ordered chronologically, from more ancient to newer.](`r params$freyja`){#fig-freyja}

### Phylogeny and temporal signal
### Phylogenetic reconstruction

A maximum likelihood tree of the target and context samples has been
build using [IQTREE](http://www.iqtree.org/).
built using [IQTREE](http://www.iqtree.org/).
The target samples `r stats[["monophyly"]]` monophyletic. The clade
that contains all the target samples is supported by a **UFBoot** score of
$`r stats[["boot"]]`$% and a **SH-aLRT** score of $`r stats[["alrt"]]`$% (@fig-tree_ml).
Expand All @@ -129,22 +129,6 @@ $`r stats[["boot"]]`$% and a **SH-aLRT** score of $`r stats[["alrt"]]`$% (@fig-t
target datasets and their context samples. The clade that contains the target
samples is squared in red.](`r params$tree_ml`){#fig-tree_ml}

A neighbor-joining tree has been constructed using pairwise distances
between target samples (@fig-tree), based on the allele frequencies measured from
read mappings.

![Neighbor-joining tree based on the pairwise allele
frequency-weighted distances.](`r params$tree`){#fig-tree}

Root-to-tip distances measured on this tree have been correlated
with time, obtaining a $R^2$ of **`r correlation`** and a p-value of
$`r p_value_lm`$. The estimated substitution rate is **`r sub_rate`**
substitutions per year (@fig-tempest).

![Scatterplot depicting the relationship between root-to-tip
distances and the number of days passed since the first sample. The red
line shows the linear model fit.](`r params$tempest`){#fig-tempest}

### Nucleotide diversity comparison

Nucleotide diversity (π) has been calculated for $`r div_values[["boot.reps"]]`$ random
Expand Down Expand Up @@ -206,7 +190,7 @@ at the bottom, and the latest, at the top.](`r params$SNV_s`){#fig-SNV_s}

:::

### Time dependency for the intra-host nucleotide variants
### Temporal signal of the intra-host mutations

The correlation of the allele frequency of each NV with the time since the
initial sampling has been calculated (@fig-volcano).
Expand All @@ -223,6 +207,22 @@ nucleotide variants with a significant correlation with time and sites with more
than two possible states. Each subplot depicts the progression of the allele
frequencies in time for a given genome position.](`r params$panel`){#fig-panel}

A neighbor-joining tree has been constructed using pairwise distances
between target samples (@fig-tree), based on the allele frequencies measured from
read mappings.

![Neighbor-joining tree based on the pairwise allele
frequency-weighted distances.](`r params$tree`){#fig-tree}

To estimate the substitution rate, root-to-tip distances measured on the previous
tree (@fig-tree) have been correlated with time, obtaining a $R^2$ of
**`r correlation`** and a p-value of $`r p_value_lm`$. The estimated substitution
rate is **`r sub_rate`** substitutions per year (@fig-tempest).

![Scatterplot depicting the relationship between root-to-tip
distances and the number of days passed since the first sample. The red
line shows the linear model fit.](`r params$tempest`){#fig-tempest}

### Correlation between alternative alleles

To detect possible interactions between mutations, pairwise correlation between allele
Expand Down

0 comments on commit 7a6e665

Please sign in to comment.