Hi Rouyer,
You can redefine dpss.taper as follows
dpss.taper.2 <- function (n, k, nw = 4, nmax = 2^(ceiling(log(n, 2))))
{
if (n > nmax)
stop("length of taper is greater than nmax")
w <- nw/n
if (w > 0.5)
stop("half-bandwidth parameter (w) is greater than 1/2")
if (k <= 0)
stop("positive dpss order (k) required")
v <- matrix(0, nrow = nmax, ncol = (k + 1))
storage.mode(v) <- "double"
out <- .Fortran("dpss", nmax = as.integer(nmax), kmax = as.integer(k),
n = as.integer(n), w = as.double(w), v = v, sig = double(k +
1), totit = integer(1), sines = double(n), vold = double(n),
u = double(n), scr1 = double(n), ifault = integer(1),
PACKAGE = "waveslim")
return(list( v=out$v[1:n, 1:k],
eigen=1+out$sig[1:k],
iter=out$totit,
n=n,
w=w,
ifault=out$ifault) );
}
or you can calculate the eigenvalues from the tapers as done in
the tridiagonal dpss calculation at
http://lib.stat.cmu.edu/sapaclisp/multitaper.lisp
Karim
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
