# Load the dataset
data(iris)
# Summary statistics
summary(iris)
## Sepal.Length Sepal.Width Petal.Length Petal.Width
## Min. :4.300 Min. :2.000 Min. :1.000 Min. :0.100
## 1st Qu.:5.100 1st Qu.:2.800 1st Qu.:1.600 1st Qu.:0.300
## Median :5.800 Median :3.000 Median :4.350 Median :1.300
## Mean :5.843 Mean :3.057 Mean :3.758 Mean :1.199
## 3rd Qu.:6.400 3rd Qu.:3.300 3rd Qu.:5.100 3rd Qu.:1.800
## Max. :7.900 Max. :4.400 Max. :6.900 Max. :2.500
## Species
## setosa :50
## versicolor:50
## virginica :50
##
##
##
# Histogram of each variable
hist(iris$Sepal.Length, main="Sepal Length Distribution", xlab="Sepal Length", col="blue")

# Scatter plot between Sepal Length and Petal Length
plot(iris$Sepal.Length, iris$Petal.Length, main="Sepal vs Petal Length", xlab="Sepal Length", ylab="Petal Length", col="red")

# Correlation matrix and heatmap
cor_matrix <- cor(iris[,1:4])
heatmap(cor_matrix)
