Genotype_value: Normalize coverage using identified/ specified normal cells...

Description Usage Arguments Value

View source: R/Genotype_value.R

Description

rho_hat: Relative coverage change for each cell in a region theta_hat: Major haplotype proportion fir each cell in a region

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
Genotype_value(
  Obj_filtered = NULL,
  type = "tumor",
  raw_counts = NULL,
  ref_counts = NULL,
  cov_adj = 1,
  ref_gtv = NULL,
  mincell = NULL,
  qt_filter = TRUE,
  cell_filter = TRUE,
  refr = TRUE
)

Arguments

Obj_filtered

An Alleloscope object with theta_hat info in the rds_list and identified/ specified normal cells and a normal region

type

Specify whether the sample is a "tumor" or "cellline". If "type" is a "cellline", param "ref_counts" needs to be specified for normal sample.

raw_counts

(required) A large binned coverage matrix (m1 bin by n1 cell) with values being read counts for all chromosomal regions of tumor sample.

ref_counts

(required only when type = "cellline") A binned coverage matrix (m2 bin by n2 cell) with values being read counts for all chromosomal regions of normal sample. n2 can be 1 for bulk sample.

cov_adj

An integer for coverage adjustment for tumor cells.

ref_gtv

A reference "genotype_values" (from scDNA-seq) to help with rho_i estimation.

mincell

An integer to filter out regions with minimum number of cells.

cell_filter

Logical (TRUE/ FALSE). Whether or not to exclude low quality cells in the output matrix.

refr

Logical (TRUE/ FALSE). Whether or not to use diplid region for normalization (otherwise, cell size is used).

Value

(rho_hat, theta_hat) of each cell for all region in the "genotype_values". Every 2 columns in the genotype_table are (rho_hat, theta_hat) of each region. Each row is a cell.


seasoncloud/Alleloscope documentation built on June 18, 2021, 7:24 p.m.