-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathcallAlleleSpecificCNfromHSCN.Rd
48 lines (41 loc) · 1.61 KB
/
callAlleleSpecificCNfromHSCN.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/callASCN.R
\name{callAlleleSpecificCNfromHSCN}
\alias{callAlleleSpecificCNfromHSCN}
\title{Call allele specific copy number in single cell datasets}
\usage{
callAlleleSpecificCNfromHSCN(
hscn,
eps = 1e-12,
maxCN = NULL,
selftransitionprob = 0.95,
progressbar = TRUE,
ncores = 1,
fillmissing = TRUE
)
}
\arguments{
\item{hscn}{hscn object from callHaplotypeSpecificCN}
\item{eps}{default 1e-12}
\item{maxCN}{maximum copy number to infer allele specific states, default=NULL which will use the maximum state from CNbins}
\item{selftransitionprob}{probability to stay in the same state in the HMM, default = 0.999, set to 0.0 for an IID model}
\item{progressbar}{Boolean to display progressbar or not, default = TRUE, will only show if ncores == 1}
\item{ncores}{Number of cores to use, default = 1}
\item{fillmissing}{For bins with missing counts fill in values based on neighbouring bins}
}
\value{
allele specific copy number object which includes dataframe similar to input with additional columns which include
* `A` (Major allele copy number)
* `B` (Minor allele copy number)
* `state_AS_phased` (phased state of the form A|B )
* `state_AS` (mirrored state of the form A|B)
* `LOH` (is bin LOH or not)
* `state_phase` (state describing which is the dominant allele and whether it is LOH or not)
* `state_BAF` (binned discretized BAF value calculated as B / (A + B))
}
\description{
Call allele specific copy number in single cell datasets
}
\details{
In the allele specific copy number inference A is always > B and state_AS_phased == state_AS
}