[R] matrix^(-1/2)

2009-12-04 Thread Greg Dropkin
re [R] matrix^(-1/2) re the discussion in November on this thread. I don't know about expm but the problem must be equivalent to solve(B^(1/2)) and a solution will exist iff B is invertible and has a square root A with A%*%A = B. For 2x2 matrices necessary and sufficient conditions for B to have

Re: [R] matrix^(-1/2)

2009-11-01 Thread David Winsemius
On Oct 31, 2009, at 9:33 PM, David Winsemius wrote: On Oct 31, 2009, at 4:39 PM, Kajan Saied wrote: Dear R-Help Team, as a R novice I have a (maybe for you very simple question), how do I get the following solved in R: Let R be a n x n matrix: \mid R\mid^{-\frac{1}{2}} solve(A) gives

Re: [R] matrix^(-1/2)

2009-11-01 Thread spencerg
A question, a comment, and an alternative answer to matrix^(-1/2): QUESTION: What's the status of the expm package, mentioned in the email you cited from Martin Maechler, dated Apr 5 19:52:09 CEST 2008? I tried both install.packages('expm') and

Re: [R] matrix^(-1/2)

2009-11-01 Thread Charles C. Berry
On Sun, 1 Nov 2009, spencerg wrote: A question, a comment, and an alternative answer to matrix^(-1/2): QUESTION: What's the status of the expm package, mentioned in the email you cited from Martin Maechler, dated Apr 5 19:52:09 CEST 2008? I tried both install.packages('expm') and

Re: [R] matrix^(-1/2)

2009-11-01 Thread spencerg
Hi, Chuck: Thanks very much, but why do I get package 'expm' is not available from install.packages(expm,repos=http://R-Forge.R-project.org;)? Best Wishes, Spencer Graves Charles C. Berry wrote: On Sun, 1 Nov 2009, spencerg wrote: A question, a comment, and an

Re: [R] matrix^(-1/2)

2009-11-01 Thread David Winsemius
On Nov 1, 2009, at 1:46 PM, spencerg wrote: Hi, Chuck: Thanks very much, but why do I get package 'expm' is not available from install.packages(expm,repos=http://R-Forge.R-project.org )? In my case I think it was it is because there is no 2.10 branch to either the:

Re: [R] matrix^(-1/2)

2009-11-01 Thread Charles C. Berry
On Sun, 1 Nov 2009, David Winsemius wrote: On Nov 1, 2009, at 1:46 PM, spencerg wrote: Hi, Chuck: Thanks very much, but why do I get package 'expm' is not available from install.packages(expm,repos=http://R-Forge.R-project.org;)? In my case I think it was it is because there is no

Re: [R] matrix^(-1/2)

2009-11-01 Thread David Winsemius
On Nov 1, 2009, at 3:12 PM, Charles C. Berry wrote: On Sun, 1 Nov 2009, David Winsemius wrote: On Nov 1, 2009, at 1:46 PM, spencerg wrote: Hi, Chuck: Thanks very much, but why do I get package 'expm' is not available from install.packages(expm,repos=http://R-Forge.R-

[R] matrix^(-1/2)

2009-10-31 Thread Kajan Saied
Dear R-Help Team, as a R novice I have a (maybe for you very simple question), how do I get the following solved in R: Let R be a n x n matrix: \mid R\mid^{-\frac{1}{2}} solve(A) gives me the inverse of the matrix R, however not the ^(-1/2) of the matrix... Thank you very much in advance!

Re: [R] matrix^(-1/2)

2009-10-31 Thread David Winsemius
On Oct 31, 2009, at 4:39 PM, Kajan Saied wrote: Dear R-Help Team, as a R novice I have a (maybe for you very simple question), how do I get the following solved in R: Let R be a n x n matrix: \mid R\mid^{-\frac{1}{2}} solve(A) gives me the inverse of the matrix R, however not the