[GRASS-user] Cluster whit area?

[Dylan Beaudette]

On Wednesday 28 July 2010, Etienne DELAY (deust 2010) wrote:
> 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 ... :-)

Sorry about that, there were 2 typos on that page... Please try the updated 
version:

http://casoilresource.lawr.ucdavis.edu/drupal/node/340

Cheers,
Dylan

-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341