if (require(spdep, quiet=TRUE)){
data(vegtf)
coa1 <- dudi.coa(vegtf$veg,scannf=FALSE)
ms.coa1 <- multispati(coa1,listw=nb2listw(vegtf$nb),nfposi=2,nfnega=0,scannf=FALSE)
summary(ms.coa1)
plot(ms.coa1)
par(mfrow=c(2,2))
s.value(vegtf$xy,coa1$li[,1],area=vegtf$area,include.origin=FALSE)
s.value(vegtf$xy,ms.coa1$li[,1],area=vegtf$area,include.origin=FALSE)
s.label(coa1$c1)
s.label(ms.coa1$c1)
}
Run the code above in your browser using DataLab