-
Notifications
You must be signed in to change notification settings - Fork 15
/
Copy pathquantify.Rd
215 lines (177 loc) · 5.93 KB
/
quantify.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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/methods_quantify.R
\name{quantify}
\alias{quantify}
\alias{quantify.default}
\alias{quantify.APA_discovery}
\alias{quantify.PESCAn_discovery}
\alias{quantify.CSCAn_discovery}
\alias{quantify.ATA_discovery}
\alias{quantify.ARA_discovery}
\alias{quantify.saddle_discovery}
\alias{quantify.IIT_discovery}
\title{Quantification of results}
\usage{
quantify(discovery, ...)
\method{quantify}{APA_discovery}(
discovery,
size = 3,
metric = "median",
shape = "center_vs_quadrants",
IDX = NULL,
...
)
\method{quantify}{PESCAn_discovery}(
discovery,
size = 5,
metric = "median",
shape = "circle",
IDX = NULL,
...
)
\method{quantify}{CSCAn_discovery}(
discovery,
size = 5,
metric = "median",
shape = "circle",
IDX = NULL,
...
)
\method{quantify}{ATA_discovery}(
discovery,
size = 3,
metric = "median",
shape = "insulation",
IDX = NULL,
...
)
\method{quantify}{ARA_discovery}(discovery, size = 3, shape = "ARA", ...)
\method{quantify}{saddle_discovery}(discovery, ...)
\method{quantify}{IIT_discovery}(discovery, ...)
}
\arguments{
\item{discovery}{A \code{discovery} object as returned by GENOVA analysis functions.}
\item{...}{further arguments passed to or from other methods.
take the middle 3x3 matrix of the APA).}
\item{size}{An \code{integer} of length one to determine the size of features
of interest in bins.}
\item{metric}{Either \code{"median"} or \code{"mean"} to summarise features.}
\item{shape}{A character of length 1 specifying what shape to use. See the
section shapes.}
\item{IDX}{The \code{IDX} part of a contacts object. Used only in converting
features expressed in bins back to genomic space. This is rarely needed,
but is useful for APAs ran with extended loops where features aren't 1:1
traceable to the input.}
}
\description{
A good amount of \code{discovery} objects can be quantified.
What exactly is to be quantified differs per \code{discovery} type.
\describe{
\item{saddle disoveries}{can be used to compute compartment strengths.}
\item{ARMLA discoveries}{such as APA, PE-SCAn, ATA and ARA compare
different regions of their outputs.}
\item{IIT discoveries}{summarise their values by neighbours.}
}
}
\section{Shapes}{
The quantification of ARMLA discoveries require a shape to distinguish
regions to quantify.
\subsection{ARA, PESCAn and CSCAn}{
APA, PESCAn and CSCAn require one of the following:
\itemize{
\item{\code{"center_vs_quadrants"}}
\item{\code{"center_vs_rest"}}
\item{\code{"circle"}}
}
The \code{size} parameter determines the number of bins of the central
foreground.
In the illustrations below, red is considered foreground and blue is
considered background.
The \code{"center_vs_quadrants"} option does not include region directly
horizontal or vertical of the centre as background.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-center-vs-quandrants.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-center-vs-quandrants.png}\out{\end{center}}
}
The \code{"center_vs_rest"} option sees everything but the centre as background.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-center-vs-rest.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-center-vs-rest.png}\out{\end{center}}
}
The \code{"circle"} option is like \code{"center_vs_rest"} but rounds corners
of the central foreground. Note that for \code{size <= 3} these two options
are equivalent.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-circle.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-circle.png}\out{\end{center}}
}
}
\subsection{ATA}{
ATA requires one of the following:
\itemize{
\item{\code{"insulation"}}
\item{\code{"cornerpeak"}}
\item{\code{"checker"}}
}
In the illustrations below, red is considered foreground and blue is
considered background. The line indicates the diagonal.
The \code{"insulation"} option compares within-TAD contacts to between-TAD
contacts.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-insulation.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-insulation.png}\out{\end{center}}
}
The \code{"cornerpeak"} option compares the intersection of boundaries versus
within-TAD contacts.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-cornerpeak.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-cornerpeak.png}\out{\end{center}}
}
The \code{"checker"} option compares the within-TAD contacts to the
contacts between it's immediate neighbours.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-checker.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-checker.png}\out{\end{center}}
}
}
\subsection{ARA}{
ARA requires one of the following:
\itemize{
\item{\code{"ARA"}}
\item{\code{"stripes"}}
}
The \code{size} argument controls the width of the stripes.
The \code{"ARA"} option reports about both 3' and 5' stripes and regions as
well as the bins that span the locus. They are indicated in different colours
below.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-ARA.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-ARA.png}\out{\end{center}}
}
The \code{"stripes"} options reports the values and distances of the stripes.
The 5' distances are encoded as negative, whereas 3' distances are positive.
\if{html}{
\out{<div style="text-align: center">}\figure{quant-stripes.png}{options: style="width:62px;max-width:100\%;"}\out{</div>}
}
\if{latex}{
\out{\begin{center}}\figure{quant-stripes.png}\out{\end{center}}
}
}
}
\examples{
NULL
}