[Robast-commits] r202 - pkg/RobLox/R

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

Author: stamats
Date: 2008-11-03 07:31:42 +0100 (Mon, 03 Nov 2008)
New Revision: 202

Modified:
   pkg/RobLox/R/rlsOptIC_AL.R
   pkg/RobLox/R/roblox.R
Log:
small bug in computation of as. covariance

Modified: pkg/RobLox/R/rlsOptIC_AL.R
===================================================================
--- pkg/RobLox/R/rlsOptIC_AL.R	2008-11-03 06:30:31 UTC (rev 201)
+++ pkg/RobLox/R/rlsOptIC_AL.R	2008-11-03 06:31:42 UTC (rev 202)
@@ -176,7 +176,7 @@
                 a1 <- A[1, 1]/sdneu^2
                 a3 <- A[2, 2]/sdneu^2
                 a2 <- a[1]/sd/a3 + 1
-                asVar <- sd^2*.ALrlsVar(b = b/sd, a1 = a1, a2 = a2, a3 = a3)
+                asVar <- sdneu^2*.ALrlsVar(b = b/sdneu, a1 = a1, a2 = a2, a3 = a3)
                 res <- list(A = A, a = sdneu*cent(IC)/sdalt, b = b, d = NULL,
                             risk = list(asMSE = mse, asBias = b, trAsCov = mse - r^2*b^2,
                                         asCov = asVar), 

Modified: pkg/RobLox/R/roblox.R
===================================================================
--- pkg/RobLox/R/roblox.R	2008-11-03 06:30:31 UTC (rev 201)
+++ pkg/RobLox/R/roblox.R	2008-11-03 06:31:42 UTC (rev 202)
@@ -156,7 +156,7 @@
     a1 <- A1/est[2]^2
     a3 <- A2/est[2]^2
     a2 <- a[2]/est[2]/a3 + 1
-    asVar <- est[2]^2*.ALrlsVar(b = b, a1 = a1, a2 = a2, a3 = a3)
+    asVar <- est[2]^2*.ALrlsVar(b = b/est[2], a1 = a1, a2 = a2, a3 = a3)
 
     return(list(est = est, A1 = A1, A2 = A2, a = a, b = b, asvar = asVar))
 }