annotate spp/man/get.binding.characteristics.Rd @ 6:ce08b0efa3fd draft

Uploaded
author zzhou
date Tue, 27 Nov 2012 16:11:40 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
1 \name{get.binding.characteristics}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
2 \alias{get.binding.characteristics}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
3 %- Also NEED an '\alias' for EACH other topic documented here.
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
4 \title{ Calculate characteristics of observed DNA-binding signal from
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
5 cross-correlation profiles }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
6 \description{
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
7 The methods calculates strand cross-correlation profile to determine binding
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
8 peak separation distance and approximate window size that should be used
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
9 for binding detection. If quality scores were given for the tags,
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
10 which quality bins improve the cross-correlation pattern.
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
11 }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
12 \usage{
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
13 get.binding.characteristics(data, srange = c(50, 500), bin = 5, cluster = NULL, debug = F, min.tag.count = 1000, acceptance.z.score = 3, remove.tag.anomalies = T, anomalies.z = 5,accept.all.tags=F)
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
14 }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
15 %- maybe also 'usage' for other objects documented here.
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
16 \arguments{
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
17 \item{data}{ Tag/quality data: output of \code{read.eland.tags} or similar function }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
18 \item{srange}{ A range within which the binding peak separation is
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
19 expected to fall. Should be larger than probe size to avoid artifacts. }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
20 \item{bin}{ Resolution (in basepairs) at which cross-corrrelation
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
21 should be calculated. bin=1 is ideal, but takes longer to calculate. }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
22 \item{cluster}{ optional snow cluster for parallel processing }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
23 \item{debug}{ whether to print debug messages }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
24 \item{min.tag.count}{ minimal number of tags on the chromosome to be
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
25 considered in the cross-correlation calculations }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
26 \item{acceptance.z.score}{ A Z-score used to determine if a given tag
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
27 quality bin provides significant improvement to the strand cross-correlation }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
28 \item{remove.tag.anomalies}{ Whether to remove singular tag count peaks prior to
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
29 calculation. This is recommended, since such positions may distort the
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
30 cross-correlation profile and increase the necessary computational time. }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
31 \item{anomalies.z}{ Z-score for determining if the number of tags at a
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
32 given position is significantly higher about background, and should be
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
33 considered an anomaly.}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
34 \item{accept.all.tags}{ Whether tag alignment quality calculations
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
35 should be skipped and all available tags should be accepted in the
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
36 downstream analysis.}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
37 }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
38 \value{
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
39 \item{cross.correlation }{ Cross-correlation profile as an $x/$y data.frame}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
40 \item{peak }{Position ($x) and height ($y) of automatically detected
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
41 cross-correlation peak.}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
42 \item{whs} { Optimized window half-size for binding detection (based
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
43 on the width of the cross-correlation peak) }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
44 \item{quality.bin.acceptance} { A list structure, describing the
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
45 effect of inclusion of different tag quality bins on
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
46 cross-correlation, and a resolution on which bins should be
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
47 considered.
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
48 \item{informative.bins} { A boolean vector indicating whether the
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
49 inclusion of tags from the tag quality bin specified in the name
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
50 attribute significantly increases cross-correlation profile near
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
51 the peak.}
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
52 \item{quality.cc} { A list giving the cross-correlation profile
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
53 after the inclusion of the tags from different quality bins }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
54 }
ce08b0efa3fd Uploaded
zzhou
parents:
diff changeset
55 }