[Robast-commits] r414 - branches/robast-0.8/pkg/ROptEst/R

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

Author: ruckdeschel
Date: 2010-09-02 22:35:14 +0200 (Thu, 02 Sep 2010)
New Revision: 414

Modified:
   branches/robast-0.8/pkg/ROptEst/R/getAsRisk.R
Log:
ROptEst: fixed a buglet in asL1 asL4 risks

Modified: branches/robast-0.8/pkg/ROptEst/R/getAsRisk.R
===================================================================
--- branches/robast-0.8/pkg/ROptEst/R/getAsRisk.R	2010-09-02 20:32:21 UTC (rev 413)
+++ branches/robast-0.8/pkg/ROptEst/R/getAsRisk.R	2010-09-02 20:35:14 UTC (rev 414)
@@ -365,8 +365,8 @@
         else{
              s <- getInfV(L2deriv, neighbor, biastype, clip, cent, stand)
              r <- neighbor at radius
-             b <- clip
-             w <- r*b/s^.5
+             b <- r*clip
+             w <- b/s^.5
              pp <- 2*pnorm(w)-1
              dp <- 2*dnorm(w)
              L1 <- b*pp+s*dp
@@ -387,7 +387,7 @@
         else{
               s <- getInfV(L2deriv, neighbor, biastype, clip, cent, stand)
               r <- neighbor at radius
-              b <- clip
+              b <- r*clip
               L4 <- 3*s^2+6*s*b^2+b^4
         }
         return(list(asL4 = L4))