diff options
| author | Suraj Yerramilli | 2015-06-05 17:16:48 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2015-06-05 17:16:48 +0530 |
| commit | bcb4780f702ed6fe1b24c9f9027f7038db23b950 (patch) | |
| tree | f64fc8040de4e63777fe55a58478b5834f91d6e7 /R | |
| parent | 8d61f4692f2f0c6183c649ca19f796d9e5abf22a (diff) | |
| download | SysID-R-code-bcb4780f702ed6fe1b24c9f9027f7038db23b950.tar.gz SysID-R-code-bcb4780f702ed6fe1b24c9f9027f7038db23b950.tar.bz2 SysID-R-code-bcb4780f702ed6fe1b24c9f9027f7038db23b950.zip | |
Corrected the degrees of freedom
Diffstat (limited to 'R')
| -rw-r--r-- | R/estpoly.R | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/R/estpoly.R b/R/estpoly.R index 9fc5964..a895525 100644 --- a/R/estpoly.R +++ b/R/estpoly.R @@ -5,7 +5,7 @@ estARX <- function(data,order=c(0,1,0)){ y <- as.matrix(data$output) u <- as.matrix(data$input); N <- dim(y)[1] na <- order[1];nb <- order[2]; nk <- order[3] - nb1 <- nb+nk; n <- max(na,nb1); df <- N - n + nb1 <- nb+nk; n <- max(na,nb1); df <- N - na - nb -nk padZeros <- function(x,n) c(rep(0,n),x,rep(0,n)) yout <- apply(y,2,padZeros,n=n); |