data(mafragh)
par(mfrow = c(3,2))
s.label(mafragh$xy, inc = FALSE, neig = mafragh$neig,
sub = "Samples & Neighbourhood graph")
coa1 <- dudi.coa(mafragh$flo, scan = FALSE)
s.value(mafragh$xy, coa1$li[,1], sub = "Axis 1 - COA")
pca1 <- dudi.pca(mafragh$xy, scan = FALSE)
s.value(mafragh$xy, pca1$li[,1], sub = "Axis 1 - PCA")
s.class(pca1$li, mafragh$partition, sub = "Plane 1-2 - PCA")
s.class(coa1$li, mafragh$partition, sub = "Plane 1-2 - COA")
s.chull(mafragh$xy, mafragh$partition, optchull = 1)
par(mfrow=c(1,1))
link1 <- area2link(mafragh$area)
neig1 <- neig(mat01 = 1*(link1>0))
nb1 <- neig2nb(neig1)
par(mfrow = c(2,1))
area.plot(mafragh$area,center = mafragh$xy,clab=0.75)
area.plot(mafragh$area,center = mafragh$xy,graph=neig1)
if (require(maptools, quiet = TRUE) & require(spdep, quiet = TRUE)) {
lw1 <- apply(link1,1,function(x) x[x>0])
listw1 <- nb2listw(nb1,lw1)
coa1 <- dudi.coa(mafragh$flo, scan = FALSE, nf = 4)
ms1 <- multispati(coa1, listw1, scan = FALSE, nfp = 2, nfn = 0)
summary(ms1)
par(mfrow = c(2,2))
barplot(coa1$eig)
barplot(ms1$eig)
s.corcircle(ms1$as)
plot(coa1$li[,1], ms1$li[,1])
}
par(mfrow = c(1,1))
Run the code above in your browser using DataLab