B.1 Inicio de una sesión y de carga de datos

Después de iniciar la sesión R, cargar geoR con el comando library (o require). Si el paquete se carga correctamente aparece un mensaje.

library(geoR)

B.1.1 Archivos de datos

Normalmente, los datos se almacenan como un objeto (una lista) de la clase geodata. Un objeto de esta clase contiene obligatoriamente dos elementos:

  • $coords: las coordenadas de las posiciones de los datos.

  • $data: los valores observados de la variables.

Opcionalmente pueden tener otros elementos, como covariables y coordenadas de las fronteras de la zona de estudio.

Hay algunos conjuntos de datos incluidos en el paquete de distribución.

# data()                    # lista todos los conjuntos de datos disponibles
# data(package = "geoR")    # lista los conjuntos de datos en el paquete geoR

data(wolfcamp)              # carga el archivo de datos wolfcamp
summary(wolfcamp)
## Number of data points: 85 
## 
## Coordinates summary
##       Coord.X   Coord.Y
## min -233.7217 -145.7884
## max  181.5314  136.4061
## 
## Distance summary
##         min         max 
##   0.3669819 436.2067085 
## 
## Data summary
##      Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
##  312.1095  471.8218  547.7156  610.2845  774.1778 1088.4209

Se pueden importar directamente un archivo de datos en formato texto:

ncep <- read.geodata('ncep.txt', header = FALSE, coords.col = 1:2, data.col = 4)
# plot(ncep)
# summary(ncep)

También se puede convertir un data.frame a un objeto geodata:

ncep.df <- read.table('ncep.txt', header = FALSE)
names(ncep.df) <- c('x', 'y', 't', 'z')
# str(ncep.df)
# Nota: los datos son espacio-temporales, pero geoR sólo admite datos 2D

datgeo <- as.geodata(ncep.df, coords.col = 1:2, data.col = 4)
# plot(datgeo)
# summary(datgeo)

O objetos de datos espaciales (entre ellos los compatibles del paquete sp), por ejemplo el siguiente código crea un objeto SpatialPointsDataFrame y lo convierte a geodata:

library(sp)
load("caballa.galicia.RData")
coordinates(caballa.galicia) <- c("x","y")
proj4string(caballa.galicia) <- CRS("+proj=longlat +ellps=WGS84")

datgeo <- as.geodata(caballa.galicia["lcpue"])
# Problemas con coordenadas duplicadas  (ver ?duplicated)
# plot(datgeo)
# summary(datgeo)

En la documentación de las funciones as.geodata y read.geodata hay más información sobre cómo importar/convertir datos.