Computes the root of a symmetric and positive semidefinite matrix.
root.matrix(X)
a symmetric matrix of same dimensions as X
a symmetric and positive semidefinite matrix
X <- matrix(c(1,2,2,8), ncol=2)
test <- root.matrix(X)
## control results
X
test %*% test
Run the code above in your browser using DataLab