if (FALSE) setwd("..") library(spdep) load("dados/ctba.RData") names(ctba) spplot(ctba, names(ctba)[9:ncol(ctba)]) viz <- poly2nb(ctba) dpad <- scale(ctba@data[,9:10]) vizw <- nbcosts(viz, dpad) nbw <- nb2listw(viz, vizw) ### cria arvore geradora minima mst <- mstree(nbw, 1) plot(ctba) plot(mst, coordinates(ctba), asp=1, add=TRUE) sk1 <- skater(mst[,1:2], dpad, 1) table(sk1$groups) sk2 <- skater(sk1, dpad, 1) table(sk1$groups, sk2$groups) sk3 <- skater(sk2, dpad, 1) table(sk2$groups, sk3$groups) sk4 <- skater(sk3, dpad, 1) table(sk3$groups, sk4$groups) sk10 <- skater(sk4, dpad, 5) table(sk10$groups) par(mfrow=c(2,3), mar=c(0,0,0,0)) plot(ctba, border=gray(.7)) plot(mst, coordinates(ctba), asp=1, add=TRUE) plot(ctba, col=sk1$groups+1) plot(ctba, col=sk2$groups+1) plot(ctba, col=sk3$groups+1) plot(ctba, col=sk4$groups+1) par(mar=c(3,3,1,1), mgp=c(2,1,0), new=FALSE) plot(sk10$ssw)