• QGIS es un Sistema de Información Geográfica (SIG) de Código Abierto licenciado bajo GNU - General Public License.
  
• QGIS es un proyecto oficial de Open Source Geospatial Foundation (OSGeo).
  
• QGIS corre sobre Linux, Unix, Mac OSX, Windows y Android.
  
• QGIS soporta numerosos formatos y funcionalidades de datos vector, datos ráster y bases de datos.
  
• QGIS es el SIG Líder de Código Abierto para Escritorio.
  
• La versión actual es QGIS 2.18.13 'Las Palmas' y fue lanzada el 15.09.2017.

http://www.qgis.org

https://github.com/qgis/QGIS-Processing

• R instalado con al menos las siguientes librerías:
  sp, rgdal, maptools, raster

• QGIS instalado con el Complemento Processing activado
  

Cabecera de Script con Parámetros.
Las líneas de cabecera deben comenzar con ##. La primera hace referencia al grupo en que se encuadra el script

##Geoestadística=group

La sintaxis de los parámetros de entrada es:

<input_param>::='##'<param>=[optional] <input_type> [<default>]<EOL>
<input_type>::= (boolean | number | [long]string | file | folder |
        extent | crs | [multiple] vector [point|line|polygon] |
        [multiple] raster | table | field <parent_param> |
        selection <option>(";"<option>)+)
Ejemplos:
##Capa=vector polygon
##Area=field Capa 
##Modelo=optional selection exacto;aproximado 

La sintaxis de los parámetros de salida es:

<output_param>::='##'<param>=(output) <output_type> <EOL>
<output_type>::= (vector | raster | table)
Por ejemplo: ##output raster

Los parámetros de configuración son:

##showplots     # Muestra salida de plots
##usereadgdal   # Usa rgdal en vez de brick al abrir vectoriales
##passfilenames # Lee desde el fichero en disco, no la capa QGIS

Cuerpo de Script
Es el código en R haciendo uso de los parámetros, con una única salvedad: Una línea que se inicia con > indica que lo que le sigue, se debe enviar al visor de resultados de QGIS.

##geostatistics=group
##CapaVec=vector
##Campo=field CapaVec
##SalidaPred=output raster
##SalidaDesv=output raster
##showplots
require("automap")
require("sp")
require("raster")
require("gstat")

table=as.data.frame(CapaVec)
coordinates(table)= ~coords.x1+coords.x2
c = CapaVec[[Campo]]
# set.seed(12345)
# x=rnorm(100,100,1)
# y=rnorm(100,100,1)
# c=rnorm(100,1000,500)
# table=data.frame(x,y,c)
# coordinates(table)=~x+y

kg_result = autoKrige(c~1, table)
Prediction = raster(kg_result$krige_output,1)
stdDev = raster(kg_result$krige_output,3)
SalidaPred<-Prediction
SalidaDesv<-stdDev

#kg_result$var_model
>kg_result$var_model

modelo = tail(kg_result$var_model,1)
m = fit.variogram(kg_result$exp_var, 
                  vgm(modelo$psill, modelo$model,
                      modelo$range, modelo$kappa))
plot(kg_result$exp_var, m, main=as.character(modelo$model))

• User Manual - processing

• Calling the SpatialPosition package from within QGIS.

• ANÁLISIS ESTADÍSTICO DE DATOS ESPACIALES CON QGIS Y R

• Image Classification with RandomForests in R (and QGIS)