plot(geyser$duration,geyser$waiting)
clus <- kmeans(geyser,centers=2)
plot(geyser$waiting,geyser$duration,pch=19,cex=0.8,col=clus$cluster)
gey <- geyser
gey[,1] <- (geyser[,1]-mean(geyser[,1]))/sqrt(var(geyser[,1]))
gey[,2] <- (geyser[,2]-mean(geyser[,2]))/sqrt(var(geyser[,2]))
clus <- kmeans(gey,centers=2)
plot(geyser$waiting,geyser$duration,pch=19,cex=0.8,col=clus$cluster)
mean(geyser)
var(geyser)
tapply(geyser[,1],clus$cluster,mean)
tapply(geyser[,2],clus$cluster,mean)
tapply(geyser[,1],clus$cluster,var)
tapply(geyser[,2],clus$cluster,var)
which(clus$cluster == 1)
cl <- as.factor(clus$cluster)
x <- c(1,3,6,4,2,5)
order(x)
x[order(x)]
sort(x)
rank(x)