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")


-- 
Enzo Cocca (PhD Candidate)
Research Fellow
Università di Napoli "L'Orientale"
mail: enzo.ccc@gmail.com
cell: +393495087014

	[[alternative HTML version deleted]]

