Performing spatial-RLQ analysis with R

Performing spatial-RLQ analysis withStephane DrayFebruary 10, 2006Thisdocumentindicateshowtoperformspatial-RLQanalysis(Drayet al., 2002) in R. Commands are written in red and outputs arewritten in blue.Contents1 Introduction 12 Data 23 Spatial-RLQ analysis 61 IntroductionR is a statistical language which can be freely downloaded from internet (http://www.r-project.org/). The package spdep written by Roger Bivand is de-voted to the creation and to the manipulation of spatial weighting matrices.You have to install it and to load it into your R session:library(spdep)You need also the ade4 package (multivariate analyses) and the maptools pack-age (if data are stored in the shape le format (ArcView GIS):library(ade4)library(maptools)Lastly, download and source the le spatialrlq.R available in the softwaresection of my website.source("spatialrlq.R")2 DataWe illustrate the method using two data sets. The data consists in a zip archivecontaining two ArcView shape les (ESRI le format). Files are extracted anddata are loaded into R with the read.shape function of the maptools library.zip.file.extract("environment.dbf", "spatialrlqdata.zip")1[1] "/tmp/RtmpkpSFHd/environment.dbf"zip.file.extract("environment.sbn", "spatialrlqdata.zip")[1] "/tmp/RtmpkpSFHd/environment.sbn"zip.file.extract("environment.sbx", "spatialrlqdata.zip")[1] "/tmp/RtmpkpSFHd/environment.sbx"zip.file.extract("environment.shp", "spatialrlqdata.zip")[1] "/tmp/RtmpkpSFHd ...