diff options
| author | Suraj Yerramilli | 2015-09-06 21:39:19 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2015-09-06 21:39:19 +0530 |
| commit | 0315a4a9a15030dccaec66f364196dd3a6faafd9 (patch) | |
| tree | 495dee92ee7784d3390a3630ad4cb637113abc2a | |
| parent | e307b4fc308590993437839eea90ea2d6fca0259 (diff) | |
| download | SysID-R-code-0315a4a9a15030dccaec66f364196dd3a6faafd9.tar.gz SysID-R-code-0315a4a9a15030dccaec66f364196dd3a6faafd9.tar.bz2 SysID-R-code-0315a4a9a15030dccaec66f364196dd3a6faafd9.zip | |
Small change
| -rw-r--r-- | R/estpoly.R | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/R/estpoly.R b/R/estpoly.R index 4bb1ac9..f35eceb 100644 --- a/R/estpoly.R +++ b/R/estpoly.R @@ -163,11 +163,10 @@ arx <- function(x,order=c(0,1,0)){ vcov <- sigma2 * chol2inv(qx$qr) - model <- arx(A = c(1,coef[1:na]),B = coef[na+1:(nb+1)], + model <- idpoly(A = c(1,coef[1:na]),B = coef[na+1:(nb+1)], ioDelay = nk,Ts=deltat(x)) - est <- estPoly(coefficients = model,vcov = vcov, sigma = sqrt(sigma2), + estPoly(coefficients = model,vcov = vcov, sigma = sqrt(sigma2), df = df,fitted.values=(X%*%coef)[1:N,], residuals=(Y-X%*%coef)[1:N,],call=match.call(),input=u) - est }
\ No newline at end of file |