diff options
| author | Suraj Yerramilli | 2015-12-30 13:53:35 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2015-12-30 13:53:35 +0530 |
| commit | b6cef1d9634b735b4424c4d09488ea4cf3e8ad6f (patch) | |
| tree | 9847e61c4f74f05ae241733b20300b6878425163 | |
| parent | 06e9e9969f5da6b66a1acdfb62bfa3b01cfa76e0 (diff) | |
| download | SysID-R-code-b6cef1d9634b735b4424c4d09488ea4cf3e8ad6f.tar.gz SysID-R-code-b6cef1d9634b735b4424c4d09488ea4cf3e8ad6f.tar.bz2 SysID-R-code-b6cef1d9634b735b4424c4d09488ea4cf3e8ad6f.zip | |
modifying oe structure
| -rw-r--r-- | R/estpoly.R | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/R/estpoly.R b/R/estpoly.R index 97175ee..79b597f 100644 --- a/R/estpoly.R +++ b/R/estpoly.R @@ -346,7 +346,7 @@ armax <- function(x,order=c(0,1,1,0),options=optimOptions()){ #' plot(mod_oe) # plot the predicted and actual responses #' #' @export -oe <- function(x,order=c(1,1,0)){ +oe <- function(x,order=c(1,1,0),options=optimOptions()){ require(signal) y <- outputData(x); u <- inputData(x); N <- dim(y)[1] nb <- order[1];nf <- order[2]; nk <- order[3]; @@ -369,7 +369,7 @@ oe <- function(x,order=c(1,1,0)){ e <- ts(l$residuals,start = start(y),deltat = deltat(y)) model <- idpoly(B = theta[1:nb],F1 = c(1,theta[nb+1:nf]), - - ioDelay = nk,Ts=deltat(x)) + ioDelay = nk,Ts=deltat(x)) estpoly(sys = model,stats=list(vcov = l$vcov, sigma = l$sigma), fitted.values=y-e,residuals=e,call=match.call(),input=u, |