[R] Error on easy way for JoSAE Package
Ana-Maria Dobre
ana24maria at gmail.com
Fri May 11 09:01:24 CEST 2012
Dear Mr/Mrs,
I have got the the data sets in the attachement (txt).
The fact is that on applying the easy method for JoSAE package on these
data sets I got the following error:
> result <- eblup.mse.f.wrap(domain.data = d.data, lme.obj = fit.lme1)
Error in solve.default(asympt.var.covar) :
Lapack routine dgesv: system is exactly singular
Could you please explain me how to solve this error or what kind of changes
i have to do on my data sets?
Thank you.
Kind regars,
Ana Dobre
-------------- next part --------------
>head(amigo)
siruta med migr y_m y_f ydirect mediu_um mediu_rm mediu_uf mediu_rf
1 1 1017 1 3 3 0 301 301 0 0 0
2 1 1071 3 0 0 0 0 0 0 0 0
3 1 1151 1 0 0 0 0 0 0 0 0
4 1 1213 1 2 1 1 164 118 0 46 0
5 1 1348 1 14 5 9 0 0 0 0 0
6 1 1455 1 0 0 0 0 0 0 0 0
female male emig_coef c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 c12 c13 c14 c15 c16
1 0 301 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 46 118 46 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
c17 c18 c19 c20 c21 c22 c23 c24 c25 c26 c27 c28 c29 c30 c31 c32 c33 c34 c35
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 46 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 0 0 0 0 1 0 0 0 0 0 1 0 3 1 0 0 0 0 2
6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
c36 c37 c38 c39 c40 c41 c42 c43 c44 c45 c46 c47 c48 c49 c50 c51 c52 c53 c54
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 2 0 1 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0
6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
c55 c56 c57 c58 c59 c60 c61 c62 c63 c64 c65 c66 c67 c68 c69 c70 c71 c72 c73
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
c74 c75 c76 c77 c78 c79 c80 c81 c82 c83 c84 c85 c86 c87 c88 c89 c90 c91 c92
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
c93 c94 c95 c96 c97 c98 c99 c100 c101 gr_1f gr_2f gr_3f gr_4f gr_5f eduinf_f
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
4 0 0 0 0 0 0 0 0 0 0 46 0 0 0 46
5 0 0 0 0 0 0 0 0 0 0 1 10 2 0 7
6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
edusup_f gr_1m gr_2m gr_3m gr_4m gr_5m eduinf_m edusup_m
1 0 85 80 0 136 0 301 0
2 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0
4 0 0 0 118 0 0 118 0
5 6 1 1 3 1 0 4 2
6 0 0 0 0 0 0 0 0
>head(rpl)
siruta med pop mediu_um mediu_rm mediu_uf mediu_rf female male c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 c12 c13 c14 c15
1 1 1017 1 65679 31270 0 34409 0 34409 31270 257 258 260 286 295 259 245 284 267 307 399 412 562 570 607
2 1 1071 3 2584 0 1255 0 1329 1329 1255 9 10 18 15 15 16 14 16 10 14 7 16 18 14 21
3 1 1151 1 5983 2971 0 3012 0 3012 2971 20 19 26 28 26 37 35 27 47 42 38 36 44 44 44
4 1 1213 1 26919 13251 0 13668 0 13668 13251 110 101 119 115 121 97 84 96 118 142 147 148 195 238 229
5 1 1348 1 21343 10525 0 10818 0 10818 10525 99 110 118 105 112 98 101 103 99 117 125 147 162 197 158
6 1 1455 1 8372 4219 0 4153 0 4153 4219 38 40 53 34 32 37 43 51 55 52 41 62 81 70 57
c16 c17 c18 c19 c20 c21 c22 c23 c24 c25 c26 c27 c28 c29 c30 c31 c32 c33 c34 c35 c36 c37 c38 c39 c40 c41 c42 c43 c44 c45 c46 c47 c48
1 616 593 568 524 570 621 645 627 580 570 555 606 615 509 472 548 584 619 649 733 438 500 502 628 561 585 672 631 629 690 670 699 712
2 10 19 13 14 9 13 11 13 20 20 12 17 16 17 15 21 20 22 18 20 11 16 18 14 7 13 18 7 8 15 14 20 18
3 44 58 60 56 44 31 35 34 30 36 48 35 64 46 50 65 60 63 89 77 51 51 43 52 54 47 39 36 45 44 45 45 42
4 221 246 236 187 183 219 234 225 207 222 201 216 210 163 211 177 219 203 239 247 170 154 175 194 182 190 212 226 200 227 256 271 242
5 190 189 178 161 171 174 196 185 196 182 186 177 172 167 155 140 172 138 177 187 93 119 111 110 139 146 155 141 150 177 192 180 177
6 68 72 92 61 57 52 51 57 64 61 51 54 54 56 66 71 76 66 90 81 44 48 63 68 76 72 82 51 58 46 60 69 69
c49 c50 c51 c52 c53 c54 c55 c56 c57 c58 c59 c60 c61 c62 c63 c64 c65 c66 c67 c68 c69 c70 c71 c72 c73 c74 c75 c76 c77 c78 c79 c80 c81
1 579 615 543 553 503 419 406 351 280 279 227 253 209 214 257 229 231 267 209 210 196 194 198 207 167 199 179 179 153 148 122 126 104
2 10 26 21 18 16 13 18 15 16 20 12 9 11 9 14 14 13 18 14 19 19 25 27 21 22 25 28 24 21 14 16 12 14
3 37 43 40 36 40 34 24 32 24 34 18 13 28 23 32 20 24 33 25 18 25 25 25 19 19 26 18 20 12 13 15 13 14
4 254 224 237 208 184 158 172 174 136 137 113 127 138 105 133 126 113 126 106 151 129 114 120 119 106 98 93 89 77 95 78 73 74
5 147 173 165 130 159 140 137 138 118 112 96 102 95 82 120 97 111 85 102 87 100 102 95 91 85 76 79 80 57 75 67 62 46
6 56 53 59 58 61 41 56 51 42 35 44 35 38 35 34 51 39 48 39 37 42 43 25 29 34 22 34 35 14 12 22 14 13
c82 c83 c84 c85 c86 c87 c88 c89 c90 c91 c92 c93 c94 c95 c96 c97 c98 c99 c100 c101 gr_1f gr_2f gr_3f gr_4f gr_5f eduinf_f edusup_f
1 98 90 38 23 21 33 26 19 20 19 9 8 3 1 3 1 1 0 1 0 5268 5914 8519 11436 3272 29129 5280
2 9 3 5 4 3 5 1 5 2 2 0 2 1 0 1 0 0 0 0 0 213 142 244 368 362 1278 51
3 9 13 2 3 6 3 4 7 5 1 2 2 0 0 0 1 0 0 0 0 513 428 848 845 378 2785 227
4 52 51 23 19 14 13 22 17 13 12 9 3 4 1 1 0 1 0 0 1 2060 2180 2961 4563 1904 12255 1413
5 38 31 13 13 11 8 12 7 11 7 6 6 2 2 2 0 4 0 0 0 1851 1822 2243 3440 1462 9903 915
6 6 12 2 5 5 4 4 2 5 2 3 0 0 0 0 0 0 0 0 0 746 635 964 1295 513 3789 364
gr_1m gr_2m gr_3m gr_4m gr_5m eduinf_m edusup_m
1 5472 5856 7012 10718 2212 25059 6211
2 215 163 263 393 221 1170 85
3 632 443 716 945 235 2558 413
4 2131 2290 2970 4491 1369 11419 1832
5 1980 1896 2293 3309 1047 9269 1256
6 848 658 1001 1305 407 3761 458
More information about the R-help
mailing list