Change tol param in ginv

R/varExp.R

@@ -163,7 +163,7 @@ calculateVarExp <- function(std_bG, std_bI, matcor, varY, N, type) {
         stop("type must be in c(\"G\", \"I\", \"J\")", call. = FALSE)
     }
     q <- qr(matcor)$rank
-    inv <- ginv(matcor, tol = 0.1)
+    inv <- ginv(matcor, tol = 0.0001)
     return((N * (crossprod(t(crossprod(std_bG, inv)), std_bG) + 
                      crossprod(t(crossprod(std_bI, inv)), std_bI)) -
                 q) / ((N - q) * varY))