[Robkalman-commits] r89 - pkg/robKalman/R

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

Author: ruckdeschel
Date: 2025-05-15 16:52:46 +0200 (Thu, 15 May 2025)
New Revision: 89

Modified:
   pkg/robKalman/R/Util.R
Log:
Utils minimal anders formatiert

Modified: pkg/robKalman/R/Util.R
===================================================================
--- pkg/robKalman/R/Util.R	2025-05-15 14:51:46 UTC (rev 88)
+++ pkg/robKalman/R/Util.R	2025-05-15 14:52:46 UTC (rev 89)
@@ -6,8 +6,8 @@
 
 ### huberizing a vector to length b
 
-Huberize <- function(x, b, norm=EuclideanNorm, ...)
-   x*ifelse(norm(x) < b, 1, b/norm(x, ...))
+Huberize <- function(x, b, norm = EuclideanNorm, ...)
+   x * ifelse(norm(x) < b, 1, b/norm(x, ...))
     
 
 limitS <- function(S, F, Z, Q, V, tol = 10^-4, itmax = 1000)#
@@ -15,12 +15,14 @@
      {SO0 <- S + 1
       S0  <- S
       i   <- 0
+#      print(S0)
       while( (sum( (SO0 - S0)^2 ) > tol^2) && (i < itmax) )
         {i   <- i + 1
          S1  <- .getpredCov(S0, F, Q)
          SO0 <- S0
-         K   <- .getKG(S1, Z, V)
+         K   <- .getKG(S1, Z, .getDelta(S1, Z, V))
          S0  <- .getcorrCov(S1, K, Z)
+#         print(list(S0,i))
         }
      S1
      }