Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
mg14 authored Oct 14, 2019
1 parent d65c5e7 commit dd316e1
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,11 +1,11 @@
# MutationTime.R
MutationTime.R is an R script to time somatic mutations relative to clonal and subclonal copy number states and calculate
# MutationTimeR
MutationTimeR is an R package to time somatic mutations relative to clonal and subclonal copy number states and calculate
the relative timing of copy number gains.
It has been used by the PCAWG consortium to calculate mutation times of 2,778 whole genome sequencing samples. http://dx.doi.org/10.1101/161562

MutationTime.R needs the following input objects
```{R}
source("MutationTime.R")
library("MutationTimeR")
vcf <- readVcf("myvcf.vcf") # Point mutations, needs info columns t_alt_count t_ref_count
bb <- GRanges(, major_cn= , minor_cn=, clonal_frequency=purity) # Copy number segments, needs columns major_cn, minor_cn and clonal_frequency of each segment
```
Expand All @@ -31,4 +31,4 @@ bb$timing_param <- mt$P
bbToTime(bb)
```

<img src="MutationTime.png" alt="MutationTime.R output" width="50%"/>
<img src="MutationTime.png" alt="MutationTime.R output" width="50%"/>

0 comments on commit dd316e1

Please sign in to comment.