[Robast-commits] r49 - pkg/ROptEst/R

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

Author: stamats
Date: 2008-02-19 13:37:59 +0100 (Tue, 19 Feb 2008)
New Revision: 49

Modified:
   pkg/ROptEst/R/getInfRobIC_asGRisk.R
Log:
check "lower" and "upper"

Modified: pkg/ROptEst/R/getInfRobIC_asGRisk.R
===================================================================
--- pkg/ROptEst/R/getInfRobIC_asGRisk.R	2008-02-19 08:10:33 UTC (rev 48)
+++ pkg/ROptEst/R/getInfRobIC_asGRisk.R	2008-02-19 12:37:59 UTC (rev 49)
@@ -33,10 +33,10 @@
             c0.old <- c0
             c0 <- try(uniroot(getInfClip, 
 ## new
-#lower = getL1normL2deriv(L2deriv = L2deriv, cent = z)/ (1 + neighbor at radius^2), 
-                        lower = .Machine$double.eps^0.75, 
-#upper = sqrt( ( Finfo + z^2 )/(( 1 + neighbor at radius^2)^2 - 1) ), 
-                        upper = upper, 
+lower = getL1normL2deriv(L2deriv = L2deriv, cent = z)/ (1 + neighbor at radius^2), 
+upper = sqrt( ( Finfo + z^2 )/(( 1 + neighbor at radius^2)^2 - 1) ), 
+#                        lower = .Machine$double.eps^0.75, 
+#                        upper = upper, 
 ##
                         tol = tol, L2deriv = L2deriv, risk = risk, 
                         neighbor = neighbor,  biastype = biastype,
@@ -131,10 +131,10 @@
             A.old <- A
             b <- try(uniroot(getInfClip, 
 ## new
-#lower = getL1normL2deriv(L2deriv = L2deriv, cent = z, stand = A, Distr = Distr)/(1+neighbor at radius^2),
-#upper = sqrt( sum( diag(A%*%Finfo%*%t(A)) + (A%*%z)^2) / ((1 + neighbor at radius^2)^2-1)),
-                         lower = .Machine$double.eps^0.75, 
-                         upper = upper, 
+lower = getL1normL2deriv(L2deriv = L2deriv, cent = z, stand = A, Distr = Distr)/(1+neighbor at radius^2),
+upper = sqrt( sum( diag(A%*%Finfo%*%t(A)) + (A%*%z)^2) / ((1 + neighbor at radius^2)^2-1)),
+#                         lower = .Machine$double.eps^0.75, 
+#                         upper = upper, 
 ##
                          tol = tol, L2deriv = L2deriv, risk = risk, 
                          biastype = biastype, Distr = Distr, neighbor = neighbor,