R/freqdom.trace.R

In kidzik/freqdom: Frequency Domain Based Analysis: Dynamic PCA

#' @noRd
# @export
freqdom.trace = function(x){
	if(!is.freqdom(x)) stop("x must be an object of class freqdom")
  lags = freqdom.lags(x)
  res = list()
  res$values = rep(0,length(lags))
  for (i in 1:length(lags))
    res$values[i] = sum(diag(x$operators[,,i]))
  res$freq = x$freq
  res$lags = x$lags
  res
}

#' #' @export
#' tr.freqdom = freqdom.trace
#' 
#' #' @export
#' tr.timedom = freqdom.trace