knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
require(gpgr) require(devtools) require(dplyr) require(knitr)
PURPLE: Purity Ploidy Estimator (https://github.com/hartwigmedical/hmftools/tree/master/purple).
PURPLE combines B-allele frequency, read depth ratios, small variants and structural variants to estimate the purity and copy number profile of a tumor sample.
It outputs several files, some of which are displayed below.
cnv_som <- system.file("extdata/purple/purple.cnv.somatic.tsv", package = "gpgr") |> gpgr::purple_cnv_som_process()
Description
cnv_som$descr |> knitr::kable(format = "html", caption = "PURPLE Somatic CNVs (per chromosome) Columns.")
cnv_som$tab |> dplyr::slice(1:10) |> knitr::kable(format = "html", caption = "PURPLE Somatic CNVs (per chromosome) Summary Table.")
umccr_key_genes <- system.file("extdata/ref/umccr_cancer_genes_2019-03-20.tsv", package = "gpgr") cnv_som_gene <- system.file("extdata/purple/purple.cnv.gene.tsv", package = "gpgr") |> gpgr::purple_cnv_som_gene_process(g = umccr_key_genes)
Description
cnv_som_gene$descr |> knitr::kable(format = "html", caption = "PURPLE Somatic CNVs (per gene) Columns.")
cnv_som_gene$tab |> dplyr::slice(1:10) |> knitr::kable(format = "html", caption = "PURPLE Somatic CNVs (per gene) Summary Table.")
cnv_germ <- system.file("extdata/purple/purple.cnv.germline.tsv", package = "gpgr") |> gpgr::purple_cnv_germ_process()
Description
cnv_germ$descr |> knitr::kable(format = "html", caption = "PURPLE Germline CNVs (per chromosome) Columns.")
cnv_germ$tab |> dplyr::slice(1:10) |> knitr::kable(format = "html", caption = "PURPLE Germline CNVs (per chromosome) Summary Table.")
purity <- system.file("extdata/purple/purple.purity.tsv", package = "gpgr") |> gpgr::purple_purity_read() purity$summary |> knitr::kable(format = "html", caption = "PURPLE Purity Summary Table.")
kat <- system.file("extdata/purple/purple.somatic.vcf.gz", package = "gpgr") |> purple_kataegis() kat$data |> knitr::kable(format = "html", caption = "PURPLE Kataegis Table.")
Description
knitr::kable(kat$description, format = "html", caption = "Kataegis column description.")
qc <- system.file("extdata/purple/purple.qc", package = "gpgr") |> gpgr::purple_qc_read() qc$summary |> knitr::kable(format = "html", caption = "PURPLE QC Summary Table.")
pkgs_of_interest <- c("base", "gpgr") si <- gpgr::session_info_kable(pkgs_of_interest) si$si_pkg si$si_pl
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.