[Robkalman-commits] r72 - branches/robKalman_2012/pkg/robKalman/R

[noreply at r-forge.r-project.org]

Author: ruckdeschel
Date: 2013-09-13 16:04:43 +0200 (Fri, 13 Sep 2013)
New Revision: 72

Modified:
   branches/robKalman_2012/pkg/robKalman/R/recFilter4.R
Log:
?\195?\132nderungen Sitzung committed

Modified: branches/robKalman_2012/pkg/robKalman/R/recFilter4.R
===================================================================
--- branches/robKalman_2012/pkg/robKalman/R/recFilter4.R	2013-09-12 20:12:56 UTC (rev 71)
+++ branches/robKalman_2012/pkg/robKalman/R/recFilter4.R	2013-09-13 14:04:43 UTC (rev 72)
@@ -6,6 +6,21 @@
 ##
 #######################################################
 
+
+     initPsRet <- function(SSM,tt, exos){
+
+        if modell has uexo uexos = matrix(...) else uexos = NULL
+        psret0 <- new("SSPredictedRet",
+                              values=matrix(0,Model at pdim,length(tt)+1),
+                              call = vector("list",length(tt)+1 ),
+                              variance = array(0,dim=c(Model at pdim,Model at pdim,length(tt)+1)),
+                              uexo = uexos,...)
+        return(psret0)
+      }
+
+## analog initPrepRet initCsRet
+
+
 recFilter <- function (Model,
                        Obs,
                        times,
@@ -36,6 +51,20 @@
      cs <- vector("list", nrSteps)
      psRet <- vector("list", nrSteps)
      csRet <- vector("list", nrSteps)
+
+     if (prep) prepret0 <- initPrepRet(SSM,tt,exox...)
+     psret0 <- initPsRet(SSM,tt,exox...)
+     csret0 <- initCsRet(SSM,tt,exox...)
+
+     for(i in nrSteps){
+           # if exists prepRet ..
+           #     prepRet0 <- new("SSPreparedRet",values=matrix)
+
+          ### aus diesem Code eine Funktion machen
+          if (prep) prepRet[[i]] <- prepret0
+          psRet[[i]] <- psret0
+          csRet[[i]] <- csret0
+     }
      for(iStep in 1:nrStep){
          psRet[[iStep]] <- initSSPredOrFiltRet(pdim = Model at pdim,
                                    qdim = Model at qdim,
@@ -103,4 +132,8 @@
             cs <- ps
          }
     }
+    zusammenhängen von list(initRet, prepRet, psRet, csRet) ## SSOutput ist jetzt
+       eine Liste mit 4 Elementen; jedes einzelne der 4 ist wieder Liste
+       und enthält zB ideal robust -> unterschied zu Output.dia
+    
  }
\ No newline at end of file