[R-sig-Geo] Error in SpatioTemporal package: system is exactly singular!!!
Mahbubeh Parsaeian
mparsa1050 at gmail.com
Mon Nov 10 09:36:54 CET 2014
Dear all,
I am using the package SpatioTemporal and trying to analyze PM
data.When I use estimate function I receive this message:Lapack
routine dgesv: system is exactly singular: U[1,1] = 0
> est.EBOD.model <- estimate(EBOD.model, x.init,
+ type="p", hessian.all=FALSE)
Optimisation using starting value 1/2
N = 9, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 31227 |proj g|= 15
At iterate 10 f = 3036.8 |proj g|= 2.4429
iterations 17
function evaluations 31
segments explored during Cauchy searches 21
BFGS updates skipped 0
active bounds at final generalized Cauchy point 0
norm of the final projected gradient 0.0555717
final function value 3036.35
F = 3036.35
final value 3036.345750
converged
Error in diag(solve(res[[i]]$hessian)) :
error in evaluating the argument 'x' in selecting a method for
function 'diag': Error in solve.default(res[[i]]$hessian) :
Lapack routine dgesv: system is exactly singular: U[1,1] = 0
I think this error occur due to high collinearity between stations but I have
checked the correlations and they are not high:
| PM31 PM32 PM33 PM34 PM35 PM36
-------------+------------------------------------------------------
PM31 | 1.0000
PM32 | 0.3547 1.0000
PM33 | 0.6414 -0.1596 1.0000
PM34 | 0.7102 0.1937 0.6225 1.0000
PM35 | 0.2578 0.6152 -0.2569 0.3733 1.0000
PM36 | 0.4386 0.6826 0.2227 0.2510 0.1930 1.0000
I have attached R code with necessary files.
If it is possible please help me to understand the problem and find a
solution for this problem.
Thanks for your help in advance.
Kindest regards,
Mahboubeh
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20141110/b038f66e/attachment.html>
-------------- next part --------------
ID X Y ST_TYPE long lat
PM31 535737 3944857 DOE 5818.10382 35495.82329
PM32 539717 3958807 DOE 5861.32662 35621.34539
PM33 539235 3962529 DOE 5856.0921 35654.83594
PM34 572740 3955767 DOE 6219.9564 35593.99147
PM35 530506 3951224 DOE 5761.29516 35553.11355
PM36 552708 3952968 DOE 6002.40888 35568.80606
-------------- next part --------------
date PM31 PM32 PM33 PM34 PM35 PM36
"2005-01-03" 77.5699 95.2212 147.689 37.1137
"2005-01-17" 110.342 122.781 182.969 36.9129
"2005-01-31" 95.5716 72.2073 129.292 58.9851 50.1795
"2005-02-14" 114.085 130.161 149.608 87.9972 143.677
"2005-02-28" 83.8428 96.9447 123.782 67.7066 61.4356
"2005-03-14" 57.8186 65.538 69.4267 36.1276 38.3897
"2005-03-28" 71.6618 97.3639 121.942 67.8738 87.4313
"2005-04-11" 75.1343 97.8515 111.284 76.6714 74.5232
"2005-04-25" 79.1766 95.7183 112.91 71.5106 76.2023
"2005-05-09" 78.6504 103.109 111.662 74.7137 88.088
"2005-05-23" 65.184 87.7625 105.992 76.1015 85.5665
"2005-06-06" 82.9721 86.7729 109.757 77.7298 132.142
"2005-06-20" 135.263 136.418 164.165 140.193 163.272
"2005-07-04" 121.158 121.749 146.32 140.625 108.66
"2005-07-18" 98.4822 97.8361 126.504 160.444 79.5385
"2005-08-01" 94.9587 98.7085 127.415 120.845 77.7862
"2005-08-15" 96.4996 100.166 122.235 151.058 70.8648
"2005-08-29" 105.197 96.5641 127.23 140.768 72.1686
"2005-09-12" 103.536 93.1749 129.97 104.036 70.5773
"2005-09-26" 117.163 100.519 132.652 81.0315 76.647
"2005-10-10" 126.582 102.314 140.32 104.718 79.8498
"2005-10-24" 102.995 86.8288 119.014 89.4026
"2005-11-07" 74.2086 76.8305 102.563
"2005-11-21" 116.235 119.553 159.828 166.642
"2005-12-05" 144.588 139.648 183.305 147.069 144.578
"2005-12-19" 76.5928 87.4176 120.522 82.1613 91.7381 53.8913
"2006-01-02" 70.5746 97.2782 134.748 81.0197 87.4118 76.9064
"2006-01-16" 88.5054 99.1791 123.285 83.6371 89.0786 62.7051
"2006-01-30" 88.1744 110.816 130.031 78.7251 123.367 53.1701
"2006-02-13" 76.8055 84.9524 108.151 69.103 114.934 41.3044
"2006-02-27" 95.1755 113.905 120.664 86.6606 200.587 57.4599
"2006-03-13" 85.6476 101.04 109.335 82.0466 155.265 57.8504
"2006-03-27" 49.6982 70.9892 67.1267 47.7369 57.445 40.9695
"2006-04-10" 72.4127 98.3976 90.8763 69.684 55.7429
"2006-04-24" 80.4402 104.118 88.7109 82.0519 59.355
"2006-05-08" 89.9933 102.121 74.2535 96.4501 63.3243
"2006-05-22" 77.4586 90.492 78.855 60.9779 56.184
"2006-06-05" 84.5304 94.6069 90.9458 68.8621 60.14
"2006-06-19" 136.636 93.9548 115.662 93.8705 90.3535
"2006-07-03" 107.806 106.083 102.609 95.2917 83.2577
"2006-07-17" 122.409 104.043 106.983 95.5096 87.0191
"2006-07-31" 137.549 115.832 133.276 111.135 100.955
"2006-08-14" 103.747 100.364 107.231 85.0293 83.3055
"2006-08-28" 118.556 123.106 100.238 83.7716
"2006-09-11" 116.378 147.033 102.196 82.2395
"2006-09-25" 92.1979 75.6208
"2006-10-09" 70.1287
"2006-10-23"
"2006-11-06"
"2006-11-20" 80.1089
"2006-12-04" 64.9309
"2006-12-18" 87.7313
"2007-01-01" 101.774 66.2574 150.156 90.5265 71.6703 64.5003
"2007-01-15" 100.495 66.5318 170.52 102.094 134.453 68.2614
"2007-01-29" 77.9259 46.6097 98.261 54.0592 81.4183
"2007-02-12" 90.4378 53.5467 102.269 62.4913 92.2482
"2007-02-26" 76.6092 56.8426 111.855 68.8661 139.618
"2007-03-12" 51.9963 44.1585 67.795 44.5355 86.5422
"2007-03-26" 81.7353 42.128 62.4435 45.8908 47.1416
"2007-04-09" 113.552 54.0614 133.382 73.9804 80.1158
"2007-04-23" 83.2845 57.1326 111.957 74.3273 98.0926
"2007-05-07" 86.816 55.0209 116.222 77.7053 128.1
"2007-05-21" 108.149 58.9473 117.634 87.2891 134.753 20.09
"2007-06-04" 110.039 63.2283 126.817 91.4376 80.3093
"2007-06-18" 123.926 67.4053 146.125 87.9957
"2007-07-02" 104.214 63.9684 111.194 81.6195
"2007-07-16" 100.962 57.1797 72.8825 65.4598
"2007-07-30" 84.7595 57.835 71.8973
"2007-08-13" 123.358 61.1251 76.8427
"2007-08-27" 102.787 68.2816 68.9308
"2007-09-10" 105.254 61.2981 78.0429
"2007-09-24" 93.1118 54.1961 65.5058
"2007-10-08" 104.324 57.8489 76.9441
"2007-10-22" 99.9095 60.1559 76.1873
"2007-11-05" 124.083 33.9783 68.5597
"2007-11-19" 104.772 35.3028
"2007-12-03" 133.454 48.2154
"2007-12-17" 85.3491 96.2912
"2007-12-31" 122.797 112.788 134.537 127.701
"2008-01-14" 165.618 172.479 144.245 189.575
"2008-01-28" 153.621 164.707 144.473 196.324
"2008-02-11" 130.148 195.437 178.815 211.807
"2008-02-25" 109.305 169.605 119.031 231.803
"2008-03-10" 137.78 186.657 117.635 229.716
"2008-03-24" 152.794 243.888 164.935 222.352
"2008-04-07" 112.418 173.446 114.935 158.093
"2008-04-21" 106.74 159.047 113.188 164.39
"2008-05-05" 115.549 182.117 91.3219 157.709
"2008-05-19" 140.279 182.894 118.829 187.56
"2008-06-02" 161.378 207.491 133.703 168.7
"2008-06-16" 146.705 185.675 102.121 105.322
"2008-06-30" 178.716 205.289 117.286 134.265
"2008-07-14" 131.585 172.276 113.916 219.31
"2008-07-28" 125.637 163.964 100.458 214.321
"2008-08-11" 165.192 74.8103 188.636
"2008-08-25" 183.972 90.7737 207.991
"2008-09-08" 209.087 118.911 225.755
"2008-09-22" 205.156 168.917 194.58
"2008-10-06" 207.747 153.293 177.335
"2008-10-20" 187.881 123.448 133.399
"2008-11-03" 166.398 120.447 132.203
"2008-11-17" 169.18 120.247 135.522
"2008-12-01" 165.836 168.043 173.11
"2008-12-15" 204.601 186.894 150.29
"2008-12-29" 166.444 151.099 119.4 61.462
"2009-01-12" 82.4194 132.752 136.816 111.246 49.502
"2009-01-26" 74.1327 141.747 120.307 96.469 50.4087
"2009-02-09" 103.907 232.775 122.074 136.124 63.651
"2009-02-23" 87.3664 123.911 110.347 100.379 51.9734
"2009-03-09" 96.4177 187.952 92.669 114.558 66.6497
"2009-03-23" 77.8626 63.1432 90.4968 41.1579
"2009-04-06" 88.9518 129.781 83.1349 105.339 58.4464
"2009-04-20" 82.8313 130.906 76.5668 105.523 56.0572
"2009-05-04" 132.652 187.073 92.93 137.377 76.3742
"2009-05-18" 101.43 71.8219 89.7404 115.27 71.2514
"2009-06-01" 99.9719 65.8528 71.7692 107.674 82.9182
"2009-06-15" 78.1662 57.1018 54.2652 101.295 56.4237
"2009-06-29" 246.554 283.637 245.939 286.863 331.682
"2009-07-13" 204.205 103.325 141.136 118.785
"2009-07-27" 141.658 170.048 124.913 168.5 153.33
"2009-08-10" 140.164 201.447 116.624 163.202 96.5306
"2009-08-24" 121.941 183.783 98.78 150.312 93.8478
"2009-09-07" 119.194 192.189 102.802 137.66 92.8439
"2009-09-21" 114.985 178.612 84.7091 130.553 67.0556
"2009-10-05" 146.171 71.128
"2009-10-19" 114.56 52.3936
"2009-11-02" 122.596 58.8266
"2009-11-16" 80.6761 35.0138
"2009-11-30" 95.2526 52.6186
"2009-12-14" 107.174 44.0446
"2009-12-28" 78.1087 104.706 84.9622 114.828 57.801
"2010-01-11" 76.156 103.155 100.474 131.615 79.964
"2010-01-25" 63.2156 69.8904 69.3502 54.0596 31.6249
"2010-02-08" 68.6425 75.8325 107.074 46.2049 61.3008
"2010-02-22" 93.3217 94.5992 109.205 65.4482 76.6534
"2010-03-08" 119.216 130.381 137.978 59.3783 104.736 101.718
"2010-03-22" 58.6934 55.5835 45.6692 45.1942 38.7566
"2010-04-05" 85.1385 86.8483 73.8888 70.4449 64.2642
"2010-04-19" 54.4261 79.7597 78.2344 51.3323 101.886 41.4635
"2010-05-03" 63.0048 93.7391 88.9033 60.8356 109.598 86.1476
"2010-05-17" 124.055 86.5296 121.696 169.753 91.751
"2010-05-31" 118.921 144.453 109.033 91.9822 186.868 84.3003
"2010-06-14" 139.933 139.271 109.009 104.299 172.948 94.7655
"2010-06-28" 85.4189 118.788 77.3933 82.6777 191.335 71.393
"2010-07-12" 169.364 106.434 107.12 227.61 101.903
"2010-07-26" 59.8349 111.741 71.8738 73.7122 163.302 58.2537
"2010-08-09" 63.0457 93.9072 63.3552 55.6514 136.338 49.1748
"2010-08-23" 61.4978 93.0257 66.3509 60.6845 150.422 54.8798
"2010-09-06" 64.5507 102.051 69.0402 64.2097 150.851 66.2128
"2010-09-20" 60.9397 126.155 67.6608 79.8722 175.262 67.7369
"2010-10-04" 117.955 76.311 173.67 57.2803
"2010-10-18" 110.546 68.4844 164.308 59.302
"2010-11-01" 125.413 75.7955 185.27 54.9296
"2010-11-15" 147.096 104.434 242.387 76.4001
"2010-11-29" 151.028 105.019 250.323 76.3402
"2010-12-13" 109.341 71.9809 162.971 64.9604
"2010-12-27"
"2011-01-10"
"2011-01-24"
"2011-02-07"
"2011-02-21" 53.8375
"2011-03-07" 70.6412
"2011-03-21" 48.6999 50.2322 39.4448 81.8781
"2011-04-04" 137.931 131.922 114.944 181.915
"2011-04-18" 113.979 80.4552 83.0858 142.973
"2011-05-02" 119.312 99.2889 91.9776 147.397
"2011-05-16" 93.2647 88.5526 68.9065 142.218
"2011-05-30" 126.686 134.845 135.341 217.304
"2011-06-13" 95.1443 100.762 90.9106 171.738
"2011-06-27" 92.9993 118.493 95.064 179.306
"2011-07-11" 67.6446 83.6045 69.277 137.388
"2011-07-25" 73.9157 88.6678 67.0895 132.536
"2011-08-08" 79.7905 114.03 79.4362 146.848
"2011-08-22" 37.2528 108.899 69.4501
"2011-09-05" 97.3273 62.5054
"2011-09-19" 109.855 68.7606
"2011-10-03" 108.434 63.4556
"2011-10-17" 77.3514
"2011-10-31"
"2011-11-14" 68.8996
"2011-11-28" 81.0618
"2011-12-12" 101.936 57.0327 179.451
"2011-12-26" 84.9502 49.6994 154.534
"2012-01-09" 81.5517 67.814 151.578
"2012-01-23" 62.1386 53.7422 126.119
"2012-02-06" 76.3158 73.2991 144.364
"2012-02-20" 72.8109
"2012-03-05" 79.2844
"2012-03-19" 53.1878 48.6094 104.024 35.8646
"2012-04-02" 75.7904 69.301 135.773 50.7425
"2012-04-16" 74.2118 78.0367 174.378 34.7276
"2012-04-30" 57.3518 72.0551 109.647
"2012-05-14" 82.6797 132.924
"2012-05-28" 58.4923 104.919
"2012-06-11"
"2012-06-25" 73.2088
"2012-07-09" 94.2812 98.7566 60.091
"2012-07-23" 106.812 78.8137 76.1181
"2012-08-06" 98.8931 68.028 67.0282
"2012-08-20" 81.2889 64.0033 58.9723
"2012-09-03" 97.6001 73.4868 64.2404
"2012-09-17" 79.8471 68.9743
"2012-10-01" 75.6195 66.6569
"2012-10-15" 60.4909 59.4116
"2012-10-29" 76.2034 60.7136
"2012-11-12" 51.8174 46.0817
"2012-11-26" 91.786 87.2726
"2012-12-10" 69.3033 54.5673
"2012-12-24" 95.2767 94.535
"2013-01-07" 66.0259
"2013-01-21" 46.1013
"2013-02-04" 51.2419
"2013-02-18" 44.9151
"2013-03-04" 65.2624
"2013-03-18" 53.0524
"2013-04-01"
"2013-04-15"
"2013-04-29"
"2013-05-13"
"2013-05-27"
"2013-06-10"
"2013-06-24"
"2013-07-08"
"2013-07-22"
"2013-08-05"
"2013-08-19"
"2013-09-02"
"2013-09-16"
"2013-09-30"
"2013-10-14"
"2013-10-28"
"2013-11-11"
"2013-11-25"
"2013-12-09"
"2013-12-23"
-------------- next part --------------
A non-text attachment was scrubbed...
Name: R-Script.R
Type: application/octet-stream
Size: 1180 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20141110/b038f66e/attachment.obj>
More information about the R-sig-Geo
mailing list