corrected errors

1 files changed, 6 insertions, 6 deletions

diff --git a/R/estpoly.R b/R/estpoly.R
index b3f1f46..f0fa549 100644
--- a/R/estpoly.R
+++ b/R/estpoly.R
@@ -1,7 +1,7 @@
 #' @export
-estpoly <- function(model,fitted.values,residuals,options=NULL,
-                    call,stats,termination=NULL,datainfo){
-  out <- list(model=model,fitted.values=fitted.values,
+estpoly <- function(sys,fitted.values,residuals,options=NULL,
+                    call,stats,termination=NULL,input){
+  out <- list(sys=sys,fitted.values=fitted.values,
               residuals=residuals,input=input,call=call,
               stats=stats,options=options,termination=termination)
   class(out) <- "estpoly"
@@ -17,7 +17,7 @@ print.estpoly <- function(est,...){
 #' @export
 summary.estpoly <- function(object)
 {
-  model <- estpoly$sys
+  model <- object$sys
   if(model$type=="arx"||model$type=="armax"){
     coefs <- c(model$A[-1],model$B)
     na <- length(model$A) - 1; nk <- model$ioDelay; 
@@ -36,7 +36,7 @@ summary.estpoly <- function(object)
   params <- data.frame(Estimated=coefs,se=se)
   
   ek <- as.matrix(resid(object))
-  N <- nrow(ek); np <- nrow(TAB)
+  N <- nrow(ek); np <- nrow(params)
   mse <- t(ek)%*%ek/N
   fpe <- det(mse)*(1+np/N)/(1-np/N)
   
@@ -57,7 +57,7 @@ predict.estpoly <- function(model,newdata=NULL){
   if(is.null(newdata)){
     return(fitted(model))
   } else{
-    mod <- coef(model)
+    mod <- model$sys
     y <- outputData(newdata); u <- inputData(newdata)
     if(mod$type=="arx"){
       f1 <- Ma(c(rep(0,mod$ioDelay),mod$B))