We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 77e58ef commit 03f9709Copy full SHA for 03f9709
R/macro_fisher_null.R
@@ -69,10 +69,8 @@ macro_fisher_null <- function(X,
69
sm_denom <- as.numeric(as.matrix(1 + Matrix::crossprod(v,info_inverse)%*%v))
70
sm_half_num <- info_inverse%*%v
71
72
- #strictly speaking, this isn't the inverse info -- its the inverse of
+ #strictly speaking, info_inverse isn't the inverse info -- its the inverse of
73
# (an approximation to) the hessian of the augmented lagrangian
74
- info_inverse <- info_inverse #- Matrix::tcrossprod(sm_half_num)/sm_denom
75
- # the tcrossprod in the last line is the final step for the sherman-morrison update
76
77
#compute derivative of augmented lagrangian
78
lag_deriv <- lapply(1:J,
0 commit comments