data(iris)
names(iris)
summary(iris$Petal.Length)
summary(iris$Sepal.Width)
summary(iris$Petal.Length)
summary(iris$Petal.Width)
summary(iris$Species)
var(iris$Petal.Length)
plot(iris$Petal.Length, iris$Petal.Width, pch=21,bg=c("red","green","yellow")[as.numeric(iris$Species)])
par(mfrow=c(1,2))
plot(jitter(as.numeric(iris$Species)), iris$Petal.Length, pch=19, xlab="Species")
plot(iris$Petal.Length ~ iris$Species)

dados1=iris[,1:4]
b=hclust(dist(dados1))
par(mfrow=c(1,1))
plot(b)
res=cutree(b,3)
#No grupo 1 
summary(dados1$Petal.Width[res==1])
#No grupo 2 
summary(dados1$Petal.Width[res==2])
#No grupo 3 
summary(dados1$Petal.Width[res==3])