diff options
| author | Suraj Yerramilli | 2016-02-15 14:01:19 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2016-02-15 14:01:19 +0530 |
| commit | 4c167a17836565cffd41ce27bf378eaf4a42130d (patch) | |
| tree | 9cfc431b51ecfc28ee1b2b290be326d278eac8b7 | |
| parent | f6a75d66c0486600ec67b066d27621486fc9a913 (diff) | |
| download | SysID-R-code-4c167a17836565cffd41ce27bf378eaf4a42130d.tar.gz SysID-R-code-4c167a17836565cffd41ce27bf378eaf4a42130d.tar.bz2 SysID-R-code-4c167a17836565cffd41ce27bf378eaf4a42130d.zip | |
minor corrections
| -rw-r--r-- | R/estpoly.R | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/R/estpoly.R b/R/estpoly.R index ddfaf93..e39c796 100644 --- a/R/estpoly.R +++ b/R/estpoly.R @@ -261,11 +261,11 @@ armax <- function(x,order=c(0,1,1,0),options=optimOptions()){ uout <- apply(u,2,padZeros,n=n) # Initial Parameter Estimates - mod_arx <- iv4(x,c(na,nb,nk)) # fitting ARX model + mod_arx <- iv(x,c(na,nb,nk)) # fitting ARX model eps_init <- matrix(resid(mod_arx)) - mod_ma <- arima(eps_init,order=c(0,0,1),include.mean = F) - e_init <- matrix(mod_ar$residuals); e_init[is.na(e_init)] <- 0 - theta0 <- matrix(c(mod_arx$sys$A[-1],mod_arx$sys$B,mod_ar$coef)) + mod_ma <- arima(eps_init,order=c(0,0,nc),include.mean = F) + e_init <- matrix(mod_ma$residuals); e_init[is.na(e_init)] <- 0 + theta0 <- matrix(c(mod_arx$sys$A[-1],mod_arx$sys$B,mod_ma$coef)) l <- levbmqdt(yout,uout,order,e_init,obj=armaxGrad, theta0=theta0,N=N,opt=options) |