[Robast-commits] r264 - branches/robast-0.7/pkg/RandVar/inst/doc branches/robast-0.7/pkg/RobLox/R pkg/RobLox/R

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

Author: stamats
Date: 2009-03-11 08:10:32 +0100 (Wed, 11 Mar 2009)
New Revision: 264

Modified:
   branches/robast-0.7/pkg/RandVar/inst/doc/RandVar.pdf
   branches/robast-0.7/pkg/RobLox/R/rlsOptIC_TuMad.R
   pkg/RobLox/R/rlsOptIC_TuMad.R
Log:
corrected small bug in computation of bias in case of TuMad IC ...

Modified: branches/robast-0.7/pkg/RandVar/inst/doc/RandVar.pdf
===================================================================
(Binary files differ)

Modified: branches/robast-0.7/pkg/RobLox/R/rlsOptIC_TuMad.R
===================================================================
--- branches/robast-0.7/pkg/RobLox/R/rlsOptIC_TuMad.R	2009-03-07 19:17:49 UTC (rev 263)
+++ branches/robast-0.7/pkg/RobLox/R/rlsOptIC_TuMad.R	2009-03-11 07:10:32 UTC (rev 264)
@@ -36,7 +36,7 @@
     a.mad <- qnorm(0.75)
     b.mad <- 1/(4*a.mad*dnorm(a.mad))
 
-    bias <- A.loc^2*(16/25/sqrt(5)*a^5)^2 + b.mad^2
+    bias <- sqrt(A.loc^2*(16/25/sqrt(5)*a^5)^2 + b.mad^2)
 
     fct1 <- function(x){ A.loc*x*(a^2 - x^2)^2*(abs(x) < a) }
     body(fct1) <- substitute({ A.loc*x*(a^2 - x^2)^2*(abs(x) < a) },

Modified: pkg/RobLox/R/rlsOptIC_TuMad.R
===================================================================
--- pkg/RobLox/R/rlsOptIC_TuMad.R	2009-03-07 19:17:49 UTC (rev 263)
+++ pkg/RobLox/R/rlsOptIC_TuMad.R	2009-03-11 07:10:32 UTC (rev 264)
@@ -19,7 +19,7 @@
     a.mad <- qnorm(0.75)
     b.mad <- 1/(4*a.mad*dnorm(a.mad))
 
-    b.2 <- sqrt(A.loc^2*(16/25/sqrt(5)*a^5)^2 + b.mad^2)
+    b.2 <- sqrt(A.loc^2*(16/25/sqrt(5)*a^5)^2 + b.mad^2)^2
 
     return(.TuMadrlsGetvar(a = a) + r^2*b.2)
 }
@@ -36,7 +36,7 @@
     a.mad <- qnorm(0.75)
     b.mad <- 1/(4*a.mad*dnorm(a.mad))
 
-    bias <- A.loc^2*(16/25/sqrt(5)*a^5)^2 + b.mad^2
+    bias <- sqrt(A.loc^2*(16/25/sqrt(5)*a^5)^2 + b.mad^2)
 
     fct1 <- function(x){ A.loc*x*(a^2 - x^2)^2*(abs(x) < a) }
     body(fct1) <- substitute({ A.loc*x*(a^2 - x^2)^2*(abs(x) < a) },