Variogramas en software R

Variogramas en el software estadístico R

Para lograr hacer un variograma en el software estadístico R, siga la siguiente guía:

1) Instalación del software estadístico R.

Link para descargar R:

https://cran.r-project.org/bin/windows/base/

En los siguientes videos podrás ver como se instala el software R:

2)  Instalación del package Gstat.

Dentro de los packeage usados con el software R para geoestadística, los más conocidos son Gstat y también Rgeo. Los siguientes videos tutoriales tratan de como instalar el peckeage Gstat:

3) Manuales del packege Gstat.

El siguiente video tutorial muestra como realizar un Kriging ordinario con Gstat:

El siguiente video tutorial muestra los comandos básicos de Gstat:

4) Ejemplo de variograma y kriging ordinario con “meuse”.

A continuación una línea de comandos que se puede usar con el archivo "meuse" que Gstat trae de ejemplo:

• En rojo las líneas de comando ingresadas por el usuario.

>library(gstat)

>library(geoR)

>library(sp)

>datos<-read.delim(“clipboard”) ó también el comando >datos<-read.table(“clipboard”) ó en “clipboard” dirección directo del archivo.

>datos

"Se obtiene la siguiente salida"

> coordinates(datos)=~x+y

> summary(datos)

"Se obtiene la siguiente salida"

Object of class SpatialPointsDataFrame

Coordinates:

  min max

x  39 546

y  43 344

Is projected: NA

proj4string : [NA]

Number of points: 112

Data attributes:

      ppm              dist   

 0      :30   116,1034022:  1 

 0,69   :21   124,6515142:  1 

 1,1    :16   134,7664647:  1 

 1,4    :10   145,2067492:  1 

 1,6    : 6   150,7481343:  1 

 2,2    : 5   161,8177988:  1 

 (Other):24   (Other)    :106 

> variogram((ppm)~1,datos)

"Se obtiene la siguiente salida"

    np      dist     gamma dir.hor dir.ver   id

1  104  39.00000  6.600962       0       0 var1

2   98  43.00000  8.183673       0       0 var1

3  182  58.05170  9.760989       0       0 var1

4   96  78.00000  5.484375       0       0 var1

5  252  88.04493 10.676587       0       0 var1

6  156  94.42987 12.448718       0       0 var1

7  232 116.44349 10.566810       0       0 var1

8  224 126.01042 11.892857       0       0 var1

9  130 134.76646 15.457692       0       0 var1

10 332 149.81044 12.953313       0       0 var1

11 140 161.81780 11.603571       0       0 var1

12 390 175.65977 15.753846       0       0 var1

13 168 191.49128 13.169643       0       0 var1

> variogram=variogram((ppm)~1,datos)

> variofit=fit.variogram(variogram,model=vgm(1,"Sph",200,1))

> variofit

"Se obtiene la siguiente salida"

  model     psill   range

1   Nug  4.818052   0.000

2   Sph 10.283634 233.568

> plot(variogram,variofit)

 "Se obtiene la siguiente salida"

5) Manuales de Gstat y Meuse.

A continuación links para descargar los manuales de Gstat y uso de meuse:

Link para Gstat:

 http://www.slideshare.net/MauricioTics2016/manual-gstat

Link para Meuse:

http://www.slideshare.net/MauricioTics2016/manual-the-meuse-data-set