compute_snv_clonality {CLONETv2} | R Documentation |
Function to compute clonality of SNVs
Description
This function takes as input the genomic position of a SNVs and computes the percentage of genomic homogeneus cells harboring the mutation.
Usage
compute_snv_clonality(sample_id, snv_read_count, beta_table, ploidy_table,
admixture_table, error_tb = error_table, error_rate = 0.05,
n_digits = 3, n_cores = 1, annotation_style = "VEP", debug = F)
Arguments
sample_id |
the id of the analyzed sample. It must be the same value reported in column sample of tables beta_table, ploidy_table, and admixture_table |
snv_read_count |
data.frame reporting in each row the genomic coordinates of an SNV together with number of reference and alternative reads covering the position in columns rc_ref_tumor and rc_alt_tumor, respectively. See parameter annotation_style for details about column names |
beta_table |
data.frame formatted as the output of function
|
ploidy_table |
data.frame formatted as the output of function
|
admixture_table |
data.frame formatted as the output of function
|
error_tb |
data.frame that reports for each combination of coverage and number informative SNPs the expected estimation error around beta. The data.frame error_tb must contains 3 columns:
Package CLONETv2 have built in error_tb named error_table (default=error_table) |
error_rate |
expected fraction of SNV positions with outlier variant allelic fraction (default=0.05) |
n_digits |
number of digits in the output table (default=3) |
n_cores |
number of cores (default=1) |
annotation_style |
a string that corresponds to the format of the columns that describe the genomic coordinates of a SNV. Accepted values are VEP and MAF. VEP annotation describes genomic coordinates with a single column named Location. MAF format has columns Chromosome, Start_position, and End_position for each aberrant position |
debug |
return extra columns for debugging (default=F) |
Value
A data.frame that extends input table snv_read_count with columns sample, cnA, cnB, t_af, t_af_corr, SNV.clonality, and SNV.clonality.status. Columns cnA and cnB report the allele specific copy number of the genomic segment containing the SNV position. Columns t_af and t_af_corr are respectively raw and ploidy/purity adjusted tumor varian allelic fractions. SNV.clonality reports the percentage of tumor cells harboring the SNV and with allele specific copy number cnA and cnB. SNV.clonality.status column lists dicretized SNV.clonality values. Discrete states are clonal, uncertain.clonal, uncertain.subclonal, and subclonal based in threshold automatically computed on the SNV.clonality values. Empty SNV.clonality.status of an SNV indicates that clonality cannot be assessed.
Author(s)
Davide Prandi, Tarcisio Fedrizzi
Examples
## Compute SNVs clonality
snv_clonality_table_toy <- compute_snv_clonality("toy_sample",
snv_reads_toy, bt_toy, pl_table_toy, adm_table_toy)