[R] gmp: Error in solve.bigz(B) : System is singular
Bertrand Haas
bhaas at fas.harvard.edu
Wed Dec 21 01:59:44 CET 2011
With a matrix such as C I do not have any problem:
>library(gmp)
> C
V1 V2 V3 V4 V5 V6 V7
[1,] 1 0 0 0 1 0 0
[2,] 0 1 0 0 0 1 0
[3,] 0 0 1 0 0 0 1
[4,] 0 0 0 1 0 0 0
[5,] 0 0 0 0 1 0 0
[6,] 0 0 0 0 0 1 0
[7,] 0 0 0 0 0 0 1
> solve.bigz(C)
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] "1" "0" "0" "0" "-1" "0" "0"
[2,] "0" "1" "0" "0" "0" "-1" "0"
[3,] "0" "0" "1" "0" "0" "0" "-1"
[4,] "0" "0" "0" "1" "0" "0" "0"
[5,] "0" "0" "0" "0" "1" "0" "0"
[6,] "0" "0" "0" "0" "0" "1" "0"
[7,] "0" "0" "0" "0" "0" "0" "1"
But with a matrix such as B I do:
> B
V1 V2 V3 V4 V5 V9 V13
[1,] 1 1 1 1 0 0 0
[2,] 0 0 0 0 1 0 0
[3,] 0 0 0 0 0 1 0
[4,] 0 0 0 0 0 0 1
[5,] 1 0 0 0 1 1 1
[6,] 0 1 0 0 0 0 0
[7,] 0 0 1 0 0 0 0
> det(B)
[1] 1
> solve(B)
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
V1 0 -1 -1 -1 1 0 0
V2 0 0 0 0 0 1 0
V3 0 0 0 0 0 0 1
V4 1 1 1 1 -1 -1 -1
V5 0 1 0 0 0 0 0
V9 0 0 1 0 0 0 0
V13 0 0 0 1 0 0 0
> solve.bigz(B)
Error in solve.bigz(B) : System is singular
Can someone help me understand what's going on?
More information about the R-help
mailing list