[Robast-commits] r321 - branches/robast-0.7/pkg/ROptEst/R

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

Author: stamats
Date: 2009-07-15 05:55:59 +0200 (Wed, 15 Jul 2009)
New Revision: 321

Modified:
   branches/robast-0.7/pkg/ROptEst/R/getAsRisk.R
Log:
small bug in getAsRisk. Now radiusMinimaxIC works for models with trafo.

Modified: branches/robast-0.7/pkg/ROptEst/R/getAsRisk.R
===================================================================
--- branches/robast-0.7/pkg/ROptEst/R/getAsRisk.R	2009-07-02 18:07:24 UTC (rev 320)
+++ branches/robast-0.7/pkg/ROptEst/R/getAsRisk.R	2009-07-15 03:55:59 UTC (rev 321)
@@ -75,11 +75,12 @@
         comp <- .getComp(L2deriv, DistrSymm, L2derivSymm, L2derivDistrSymm)
         z.comp <- comp$"z.comp"
         A.comp <- comp$"A.comp"
+        DA.comp <- abs(trafo) %*% A.comp != 0
         
         eerg <- .LowerCaseMultivariate(L2deriv = L2deriv, neighbor = neighbor, 
              biastype = biastype, normtype = normtype, Distr = Distr, 
              trafo = trafo, z.start = z.start, A.start, z.comp = z.comp, 
-             A.comp = A.comp,  maxiter = maxiter, tol = tol)
+             A.comp = DA.comp,  maxiter = maxiter, tol = tol)
         erg <- eerg$erg
         bias <- 1/erg$value