movt                   package:reb                   R Documentation

_m_o_v_t

_D_e_s_c_r_i_p_t_i_o_n:

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

_U_s_a_g_e:

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

_A_r_g_u_m_e_n_t_s:

       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.

_D_e_t_a_i_l_s:

     'movt' acts very similar to 'movbin'

_V_a_l_u_e:

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

_A_u_t_h_o_r(_s):

     Kyle A. Furge, Ph.D., kyle.furge@vai.org and Karl J. Dykema,
     karl.dykema@vai.org

_S_e_e _A_l_s_o:

     'movbin'

_E_x_a_m_p_l_e_s:

     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")  

