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....@gmail.com
cell: +393495087014
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
