[Robast-commits] r112 - in branches/robast-0.6/pkg: ROptEst/R RobAStBase/R

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Jul 22 14:44:40 CEST 2008


Author: stamats
Date: 2008-07-22 14:44:39 +0200 (Tue, 22 Jul 2008)
New Revision: 112

Modified:
   branches/robast-0.6/pkg/ROptEst/R/getRiskIC.R
   branches/robast-0.6/pkg/RobAStBase/R/getRiskIC.R
Log:
some minor changes

Modified: branches/robast-0.6/pkg/ROptEst/R/getRiskIC.R
===================================================================
--- branches/robast-0.6/pkg/ROptEst/R/getRiskIC.R	2008-07-21 14:11:37 UTC (rev 111)
+++ branches/robast-0.6/pkg/ROptEst/R/getRiskIC.R	2008-07-22 12:44:39 UTC (rev 112)
@@ -15,7 +15,7 @@
                                  neighbor = "missing",
                                  L2Fam = "L2ParamFamily"),
     function(IC, risk, L2Fam){
-        Cov <- IC at Risks[["asCov"]]  ; print(L2Fam)      
+        Cov <- IC at Risks[["asCov"]]
         return(list(asCov = list(distribution = .getDistr(L2Fam), value = Cov)))
     })
 
@@ -26,9 +26,9 @@
 setMethod("getBiasIC", signature(IC = "HampIC",
                                  neighbor = "UncondNeighborhood"),
     function(IC, neighbor, L2Fam,...){
-        if(missing(L2Fam)) 
-           {misF <- TRUE; 
-            L2Fam <- eval(IC at CallL2Fam)}
+        if(missing(L2Fam))
+            L2Fam <- eval(IC at CallL2Fam)
+
         return(list(asBias = list(distribution = .getDistr(L2Fam), 
                     neighborhood = neighbor at type, value = IC at Risks$asBias$value)))
     })

Modified: branches/robast-0.6/pkg/RobAStBase/R/getRiskIC.R
===================================================================
--- branches/robast-0.6/pkg/RobAStBase/R/getRiskIC.R	2008-07-21 14:11:37 UTC (rev 111)
+++ branches/robast-0.6/pkg/RobAStBase/R/getRiskIC.R	2008-07-22 12:44:39 UTC (rev 112)
@@ -52,8 +52,8 @@
         if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
             stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
 
-        trCov <- getRiskIC(IC, risk = asCov(), L2Fam = L2Fam)$trCov        
-        trCov$value <- sum(diag(std %*% trCov$value))
+        trCov <- getRiskIC(IC, risk = asCov(), L2Fam = L2Fam)$asCov
+        trCov$value <- sum(diag(trCov$value))
 
         prec <- checkIC(IC, L2Fam, out = FALSE)
         if(prec > tol)



More information about the Robast-commits mailing list