[GRASS-user] Cluster whit area?

Etienne DELAY (deust 2010) etienne.delay at etu.unilim.fr
Wed Jul 28 09:12:44 EDT 2010


Sorry, an error in copy / past with a capital letter here or one should
not....
But an other probleme so :
s <- stepFlexclust(y, k=2:10, nrep=20) it's ok
y.pam <- pam(y, 5, stand=TRUE)
> plot(y$easting, y$northing, col=y.pam$clustering, main="Bugsites
Spatial Clustering, 5 classes", cex=0.5, pch=16, xlab="Easting",
ylab="Northing")
> points(y.pam$medoids, pch=15, col=1:5, cex=1.25) OK

but

for(i in 1:5)
{
segments(x0=y.pam$medoids[i,][[1]], y0=y.pam$medoids[i,][[2]],
x1=y$easting[y$cluster == i], y1=y$northing[y$cluster ==i], col=i, lty=3)
}

Erreur in segments(x0 = y.pam$medoids[i, ][[1]], y0 = y.pam$medoids[i,  :
impossible to mix the coordinates of zero length with other

So what else ... :-)

-- 
Cordialement

Etienne DELAY
Deust "Webmaster et gest. d'intranet" (IDSkype : etienne.delay.tic)
Master IUP "valorisation du patrimoine et aménagement du territoire"


More information about the grass-user mailing list