Description Usage Arguments Details Value Author(s) Examples
This function analyzes ordered data series to identify regional biases using an moving (running) approximated binomial test.
1 |
v |
data vector |
span |
numeric vector. Each element is used to define
the number of points to include when the approximated binomial test
is applied to |
summarize |
function that is used to summarize the results from
multiple spans. if NULL, a matrix with |
movbin
applies a moving binomial test to sequential windows of
elements of v
. Within each span a z-score from an approximated
binomial is computed such that z=(2*r - n)/sqrt(n)
where
r
is the number of positive relative gene expression values and
n
is the number of non-zero values within each window.
For convenience, this function allows for the specification of multiple
window sizes using the span
argument. The result of a
movbin
call will generate a matrix with length(span)
rows
and length(v)
columns. Each row of the matrix represents the
data generated from each span. This matrix can be returned or the matrix
from can be condensed to a single vector of length v by
applying a summary function summarize
to the matrix columns.
Either a matrix or a vector containing the summarized z-scores from the applied binomial test.
Kyle A. Furge, Ph.D., kyle.furge@vai.org and Karl J. Dykema, karl.dykema@vai.org
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | x <- c(rnorm(50,mean=1),rnorm(50,mean=-1),rnorm(100))
layout(1:2)
plot(x,type="h",ylim=c(-5,5))
## apply the approximated binomial with a single span
mb <- movbin(x,span=25,summarize=NULL)
lines(mb[1,])
## try a few different span ranges
mb <- movbin(x,span=c(10,25,50),summarize=NULL)
lines(mb[1,]) ## span of 10
lines(mb[2,]) ## span of 25
lines(mb[3,]) ## span of 50
## average the results from the different spans
plot(x,type="h",ylim=c(-5,5))
mb <- movbin(x,span=c(10,25,50),summarize=mean)
lines(mb,col="blue")
mb <- movbin(x,span=c(10,25,50),summarize=median)
lines(mb,col="red")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.