Greetings, R helpers,
I have been using the Tps() function in the fields package to model
response surfaces for some residuals data.
But i don´t know how i can get more resolution and fill all the shape area?
I hope someone can help me.
Thank you,
Dominic
outresi <- Tps(x=as.matrix(ana[,3:4]),Y=resi)
preresi <- predict.surface(outresi)
image(preresi)
resi <- structure(c(1.18650257765518, 0.922870034874496, 1.40456497838008,
-1.98869637540279, 0.652835483390163, 3.38161411979703, 1.5301074186747,
-0.409663707199911, 1.34434770758852, -0.416787026362009, 4.18757416663512,
3.28025478124903, -5.17448326255257, 2.95828686941439, -0.782530617455976,
1.82121460192511, 0.203539209775968, 0.508322787328424, -1.94498277323241,
-4.89815558191742, 3.36600081231205, -3.43524939015204, -3.05944133924843,
1.19512325854811, -3.04309148543589, 1.40717743356163, -3.35495846881043,
2.72864346101486, 1.37230147197744, -0.106841569047528, -0.119344066154028,
-0.608350506308074, 1.44743280454379, 3.30716673292547, 0.874190165066352,
-0.916075316748891, 1.37181996156702, -1.2930355920363, -4.54134505448977,
1.80012571483517, -0.721990604703457, -2.65811792871444, 2.14103809342353,
-2.41809548989305, -1.81952424936354, 3.31004331724789, -3.30108275773548,
2.9757633493691, 0.295994951412489, 2.30145808678062, 0.422002402195169,
0.817785741134315, 3.49852828825309, 0.409779315387873, 2.30965751639276,
-1.00806803466553, 2.84451533239655, -2.35467760469766, 4.42619863347384,
2.70778104139562, 2.16518126265409, -2.88134006056208, -1.09216809797566,
-1.52993750262298, 2.65269273615996, 3.15979996437503, -0.363350001435549,
-6.39582581459343, -1.04031709170907, 1.33328272103405, -5.49726774890521,
4.08987284624148, 2.50912004072798, 1.8274159447864, -2.22385808283135,
-4.86917684419856, -0.564375776024271, 1.8708963744, -3.78120181090081,
-0.558708890799885, 4.05539617977911, -0.948307806738612,
-1.55729490567548,
1.96565015968813, 2.91601715036197, -1.25289442884952, 2.88976270172768,
-1.9916070689498, -1.48771804146596, -2.1835149830146, -2.15786396484737,
-7.13845827165779, 2.54353180775172, -3.42663298459986), .Names = c("1",
"2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",
"14", "15", "16", "17", "18", "19", "20", "21", "23", "25", "26",
"27", "28", "29", "30", "31", "32", "34", "35", "36", "37", "38",
"39", "40", "41", "42", "43", "44", "46", "47", "48", "49", "50",
"51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61",
"62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72",
"73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83",
"84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94",
"95", "96", "98", "99"))
ana <- structure(list(Long = c(521328, 643449, 547859, 495486, 510314,
570154, 651957, 591244, 653608, 616292, 562796, 669916, 648619,
670351, 607180, 519179, 610567, 625025, 631194, 640477, 531377,
618769, 684214, 516223, 623222, 524431, 511201, 560578, 641809,
590324, 497771, 550783, 502146, 648423, 574209, 510859, 549546,
639190, 625827, 585343, 598945, 682612, 489814, 485023, 527749,
560019, 513379, 589733, 623277, 583167, 526634, 613230, 570876,
518427, 660057, 636959, 558736, 505372, 594942, 594370, 532925,
658781, 654995, 615988, 522575, 555930, 639736, 522594, 612930,
576464, 538116, 547280, 566252, 540986, 525267, 498997, 598421,
536101, 516677, 640275, 544158, 597892, 664123, 587674, 518653,
540339, 670165, 632785, 629480, 657933, 523762, 520553, 534026,
585393), Lat = c(4702005, 4734442, 4801512, 4773043, 4700787,
4823868, 4772302, 4724987, 4730789, 4684126, 4696309, 4743224,
4800328, 4703493, 4644830, 4706070, 4832008, 4723028, 4833932,
4680282, 4716165, 4761493, 4690230, 4704557, 4779463, 4701030,
4649389, 4815885, 4706108, 4765782, 4711638, 4801765, 4764251,
4717609, 4644255, 4757557, 4685177, 4824470, 4812691, 4669787,
4786917, 4694030, 4738281, 4747735, 4695489, 4788103, 4798361,
4799728, 4703623, 4751251, 4658771, 4676811, 4718574, 4732806,
4730358, 4664595, 4775147, 4692421, 4688775, 4688943, 4685679,
4782554, 4822868, 4835419, 4706352, 4720490, 4780523, 4676863,
4741204, 4839683, 4799702, 4675170, 4753162, 4702027, 4771652,
4701610, 4777923, 4747354, 4694770, 4695132, 4741500, 4827467,
4691376, 4714920, 4698461, 4797835, 4758256, 4648384, 4645955,
4669863, 4675121, 4675090, 4725092, 4639803)), .Names = c("Long",
"Lat"), class = "data.frame", row.names = c("1", "2", "3", "4",
"5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15",
"16", "17", "18", "19", "20", "21", "23", "25", "26", "27", "28",
"29", "30", "31", "32", "34", "35", "36", "37", "38", "39", "40",
"41", "42", "43", "44", "46", "47", "48", "49", "50", "51", "52",
"53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63",
"64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74",
"75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85",
"86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96",
"98", "99"))
--
Dominic Royé
*
*
*Departamento de Xeografía*
*Facultade de Xeografía e Historia*
*Praza da Universidade 1*
*15782 Santiago de Compostela*
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
