From 01531784e9cb874d948c47d445c52c558a2735d2 Mon Sep 17 00:00:00 2001
From: Suraj Yerramilli
Date: Fri, 4 Sep 2015 10:24:49 +0530
Subject: corrected the phase in spa

---
 R/nonparam.R | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

(limited to 'R')

diff --git a/R/nonparam.R b/R/nonparam.R
index 7df30d9..4d8496a 100644
--- a/R/nonparam.R
+++ b/R/nonparam.R
@@ -175,9 +175,11 @@ spa <- function(data,npad=255){
   
   # Non-parametric Estimation of Spectral Densities - 
   # WOSA and Hanning window
-  gamma <- SDF(temp,method="wosa",sampling.interval = data$Ts,npad=npad)
-  freq <- seq(from=1,to=ceiling(npad/2),by=1)/ceiling(npad/2)*pi/data$Ts
-  out <- idfrd(response = gamma[,2]/gamma[,3],freq=freq,Ts= data$Ts)
+  gamma <- SDF(temp,method="wosa",sampling.interval = deltat(data),
+               npad=npad)
+  freq <- attributes(gamma)$frequency*2*pi
+  out <- idfrd(response = Conj(gamma[,2])/Mod(gamma[,3]),freq=freq,
+               Ts= deltat(data))
   return(out)
 }
 
-- 
cgit