summaryrefslogtreecommitdiff
path: root/R
diff options
context:
space:
mode:
authorSuraj Yerramilli2015-06-07 19:02:23 +0530
committerSuraj Yerramilli2015-06-07 19:02:23 +0530
commita00eafd068de773660a157794b98b4b99c5de4db (patch)
tree98973b7d28e2fd78ea72b4a0ee4d4060b3e0b014 /R
parentc7059970d5ff5c12f075d836cd7213245d791c6b (diff)
downloadSysID-R-code-a00eafd068de773660a157794b98b4b99c5de4db.tar.gz
SysID-R-code-a00eafd068de773660a157794b98b4b99c5de4db.tar.bz2
SysID-R-code-a00eafd068de773660a157794b98b4b99c5de4db.zip
Corrected minor errors
Diffstat (limited to 'R')
-rw-r--r--R/sim.R3
1 files changed, 2 insertions, 1 deletions
diff --git a/R/sim.R b/R/sim.R
index 3892958..6f802f3 100644
--- a/R/sim.R
+++ b/R/sim.R
@@ -42,7 +42,8 @@ sim.arx <- function(model,input,sigma=0){
# u <- apply(input,2,padLeftZeros)
for(i in n+1:length(input)){
- reg <- matrix(c(-(y[i-1:na]),u[i-nk:nb1]),ncol=na+nb1)
+ if(nk==0) v <- u[i-0:(nb-1)] else v <- u[i-nk:nb1]
+ reg <- matrix(c(-(y[i-1:na]),v),ncol=na+nb1)
y[i] <- reg%*%coef + rnorm(1,sd = sigma)
}
return(y[n+1:length(input)])