[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