Analyzing the Iris Dataset

Introduction

We’re going to analyze Fisher’s data (again)

iris = read.csv("iris.csv")
iris$Species = factor(iris$Species)
head(iris)

##   Sepal.Length Sepal.Width Petal.Length Petal.Width Species Code
## 1          5.1         3.5          1.4         0.2  setosa    1
## 2          4.9         3.0          1.4         0.2  setosa    1
## 3          4.7         3.2          1.3         0.2  setosa    1
## 4          4.6         3.1          1.5         0.2  setosa    1
## 5          5.0         3.6          1.4         0.2  setosa    1
## 6          5.4         3.9          1.7         0.4  setosa    1

cor(iris$Petal.Length, iris$Petal.Width)

## [1] 0.9628654

Let’s add a simple plot

plot(iris$Petal.Length, iris$Petal.Width, 
     pch = 16, col = iris$Species)