[R] function gstat r with pyhton

Pascal Oettli kridox at ymail.com
Thu May 2 14:51:07 CEST 2013


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.
>



More information about the R-help mailing list