## File Name: ginverse_sym.R
## File Version: 0.08
#--- code from Eugene Demidenko: book mixed effects models
ginverse_sym <- function(A, eps=1E-8)
{
# generalized inverse of a symmetric matrix A
PV <- eigen(A, symmetric=TRUE)
V0 <- IV <- PV$values
av0 <- abs(V0)
IV[ av0 > eps] <- 1/V0[ av0 > eps]
IV[ av0 <=eps] <- 0
Ainv <- PV$vectors %*% ( IV*( t(PV$vectors) ) )
return(Ainv)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.