diff options
| author | Suraj Yerramilli | 2016-03-18 18:16:03 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2016-03-18 18:16:03 +0530 |
| commit | 9a894b78aea2bfeb77042fd3f9e4598aee6223b2 (patch) | |
| tree | e973cd4f24f222ae67acf61198422abbef3e40a1 | |
| parent | e73ac91c966dc20e34a67659804faa4a39efbddd (diff) | |
| download | SysID-R-code-9a894b78aea2bfeb77042fd3f9e4598aee6223b2.tar.gz SysID-R-code-9a894b78aea2bfeb77042fd3f9e4598aee6223b2.tar.bz2 SysID-R-code-9a894b78aea2bfeb77042fd3f9e4598aee6223b2.zip | |
fixing integrated noise part
| -rw-r--r-- | R/estpoly.R | 7 |
1 files changed, 4 insertions, 3 deletions
diff --git a/R/estpoly.R b/R/estpoly.R index 814b225..d1747f3 100644 --- a/R/estpoly.R +++ b/R/estpoly.R @@ -163,11 +163,12 @@ residplot <- function(model,newdata=NULL){ #' #' @export arx <- function(x,order=c(1,1,1),lambda=0.1,intNoise=FALSE){ - y <- outputData(x); u <- inputData(x); N <- dim(y)[1] + y <- outputData(x); u <- inputData(x) if(intNoise){ - y <- apply(y,2,integfilter) - u <- apply(u,2,integfilter) + y <- apply(y,2,diff) + u <- apply(u,2,diff) } + N <- dim(y)[1] na <- order[1];nb <- order[2]; nk <- order[3] nb1 <- nb+nk-1 ; n <- max(na,nb1); df <- N-na-nb |