-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathCNV_analysis.Rd
61 lines (46 loc) · 2.03 KB
/
CNV_analysis.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
49
50
51
52
53
54
55
56
57
58
59
60
61
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/CNV_analysis.R
\name{CNV_analysis}
\alias{CNV_analysis}
\title{CNV analysis}
\usage{
CNV_analysis(
scMuffinList = NULL,
mc.cores = 1,
reference = NULL,
min_cells = 100,
min_genes = 200,
wnd_size = 100,
center_cells = TRUE,
center_genes = FALSE,
expr_lim = FALSE,
method = "mean",
na.rm = FALSE,
z.score = FALSE,
eps = NULL,
n_comp = 10,
gene_ann = NULL,
...
)
}
\arguments{
\item{scMuffinList}{scMuffinList object}
\item{mc.cores}{number of cores;}
\item{reference}{one-column data.frame with a reference expression profile; rownames must match those of genes_by_cell;}
\item{min_cells}{minimum numbver of cells in which a gene must be expressed;}
\item{min_genes}{minimun number of genes expressed in a cell;}
\item{wnd_size}{number of adjacent genes considered;}
\item{center_cells}{whether to center cells. Default is TRUE}
\item{center_genes}{whether to center genes or not. Default is FALSE. When using a reference it is useful to set this to TRUE}
\item{expr_lim}{min and max values of relative expression; by default, lower and higher whiskers returned by grDevices::boxplot.stats will be used. Set to NA or anything other value such that length(expr_lim) != 2 to disable the removal of outliers.}
\item{method}{"mean": subtract the average profile of the reference cluster to every cell (default); "min_max" (from Tirosh et al., DOI: 10.1126/science.aad0501) can be used but is not tested yet}
\item{na.rm}{whether to remove 0 values in CNV estimation}
\item{z.score}{whether to use z-scores of the cluster median CNV profile, instead of the median itself.}
\item{eps}{absolute threshold to call CNV regions.}
\item{n_comp}{number of PCA components used for clustering}
\item{gene_ann}{optional data.frame with gene annotation with mandatory columns "Chromosome", "symbol" and "pos" (genomic location). If NULL gene locations will be collected from org.Hs.eg.db}
\item{...}{arguments passed to Seurat::FindClusters()}
}
\description{
CNV analysis
}