#Charger les données spatiales pour déterminer les groupes
spa <- read.csv ('http://www.davidzeleny.net/anadat-r/data-download/DoubsSpa.csv', row.names = 1)
spa <- spa[,-8]

#Visualiser la matrice de données 
View (spa)

#Ajouter les numéros de site
numbers<-(1:30)
numbers<-numbers[!numbers%in%8] 
spa$site<-numbers

#Faire des groupes en fonction de la latitudey<82=group1, 82<y<156=group2, y>156=group3
spa.group<-ddply(.data=spa, .variables=.(x, y, site), .fun= summarise, group = if(y <= 82) 1 else if (y <= 156) 2 else 3)

#Ordonner par site
spa.group<-spa.group[with(spa.group, order(site)), ]