1 files changed, 11 insertions, 5 deletions
diff --git a/R/nonparam.R b/R/nonparam.R
index 068f9b7..7df30d9 100644
--- a/R/nonparam.R
+++ b/R/nonparam.R
@@ -125,11 +125,17 @@ plot.impulseest <- function(model,sig=0.975){
 #' 
 #' @export 
 step <- function(model){
-  title <- paste("Step Response \n From",model$x,"to",model$y)
-  stepResp <- cumsum(coef(model))
-  plot(model$lags,stepResp,type="s",xlab="Lag",ylab= model$y,
-       main = title)
-  abline(h=0)
+  par(mfrow=c(model$noutputs,model$ninputs))
+  
+  stepplot <- function(model){
+    title <- paste("Step Response \n From",model$x,"to",model$y)
+    stepResp <- cumsum(coef(model))
+    plot(model$lags,stepResp,type="s",xlab="Lag",ylab= model$y,
+         main = title)
+    abline(h=0)
+  }
+  l <- model[seq(model$noutputs*model$ninputs)]
+  p <- lapply(l,stepplot)
 }
 
 #' Estimate frequency response