function gstat r with pyhton
Hi, As you probably noticed, it is r-help list, not python-help list. Regards, Pascal
On 05/02/2013 09:38 PM, Enzo Cocca wrote:
hi,
can someone help me to put this code R in python ?
library(gstat)
VGM_PARAM_A3 <- gstat(id="bos_bison",
formula=combusto~1,locations=~coord_x+coord_y, data=archezoology_table,
nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "calcinati", strie~1,
locations=~coord_x+coord_y, archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "camoscio", cervo~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "capriolo", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "cervo", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "combusto", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "coni", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "pdi", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "stambecco", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "strie", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, "canidi", pdi~1,
locations=~coord_x+coord_y,archezoology_table, nmax = 10)
VGM_PARAM_A3 <- gstat(VGM_PARAM_A3, model=vgm(1, "Sph", 5, 0),
fill.all=TRUE)
ESV_A3 <- variogram(VGM_PARAM_A3, cutoff=9)
VARMODEL_A3 = fit.lmc(ESV_A3, VGM_PARAM_A3)
plot(ESV_A3, model = VARMODEL_A3,
pch=20, cex=0.7, col="red",
main="Linear Model of Coregionalization for A3 variables")
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.