I need to compute two equations related with trace and inverse of a around
30000 x 30000 density matrices. The equations are
-trace( W_i %**% C) and -trace(W_i %**% C %*% W_j C)
I know W_i, W_j and inverse of C. These equations are related with Pearson
estimating functions. I am trying to use R and package Matrix, but I
couldn't compute the C matrix, using solve() or chol() and chol2inv(). I do
not know with is possible using solve() to solve a system of equation and
after compute the trace. It is common to use solve function to compute
something like C^{-1} W = solve(C, W), but my equation is a little bit
different. Any help is welcome. I am thinking about to use RcppArmadillo,
but I am not sure that it is able to compute my equations.
Thank you everyone.
--
Wagner Hugo Bonat
LEG - Laboratório de Estatística e Geoinformação
UFPR - Universidade Federal do Paraná
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
