R/ews.trans.R

Defines functions ews.trans

Documented in ews.trans

ews.trans <-
function(x,scales=NULL){
  N=length(x)
  if (length(scales)==0) {
    J=floor(log(N,2))
  } else {
    J=length(scales)
    }
  res = modwt(x,filter="haar")
  out = matrix(0,N,J)
  for (j in 1:J){
    out[,j]=res@W[[j]]
  }
  for (i in 1:(J)) out[,i]=c(out[(2^i):N,i],out[1:(2^i-1),i])
  return(out^2)
}

Try the wbsts package in your browser

Any scripts or data that you put into this service are public.

wbsts documentation built on July 1, 2020, 5:23 p.m.