reb (version 1.50.0)

movt: movt

Description

This function analyzes ordered data series to identify regional biases using an moving (running) approximated t-test.

Usage

movt(v,span=NULL,summarize=mean)

Arguments

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 v. While mixed for the defaults, the span can be specified as fraction of the observation or actual sizes, but not a mixture - defaults to: seq(25,length(v)*.3,by=5)

summarize

function that is used to summarize the results from multiple spans. if NULL, a matrix with length(span) rows and length(v) columns is returned.

Value

Either a matrix or a vector containing the summarized z-scores from the applied t-test.

Details

movt acts very similar to movbin

See Also

movbin

Examples

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 <- movt(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 <- movt(x,span=c(10,25,50),summarize=mean)
lines(mb,col="blue")

mb <- movt(x,span=c(10,25,50),summarize=median)
lines(mb,col="red")  

mb <- movt(x,span=c(10,25,50),summarize=max)
lines(mb,col="green")