Description

Instructions: Please submit the pdf file on Gradescope and Rmd file on Canvas.
Part 1 (Iris)
Background: Edgar Anderson’s Iris Data
The R data description follows:
This famous (Fisher’s or Anderson’s) iris data set gives the measurements in centimeters of the variables sepal length and width and petal length and width, respectively, for 50 flowers from each of 3 species of iris. The species are Iris setosa, versicolor, and virginica.
Task
The purpose of this task is to construct a complex plot using both base R graphics and ggplot. Consider the follwoing base R plot.
# Base plot plot(iris$Sepal.Length,iris$Petal.Length,col=iris$Species,xlab=”Sepal Length”,ylab=
# loop to construct each LOBF
for (i in 1:length(levels(iris$Species))) { extract <- iris$Species==levels(iris$Species)[i]
abline(lm(iris$Petal.Length[extract]~iris$Sepal.Length[extract]),col=i) }
# Legend legend(“right”,legend=levels(iris$Species),fill = 1:length(levels(iris$Species)), cex = .
# Add points and text
points(iris$Sepal.Length[15],iris$Petal.Length[15], pch = “*”, col = “black”) text(iris$Sepal.Length[15]+.4,iris$Petal.Length[15],”(5.8,1.2)”,col=”black”) points(iris$Sepal.Length[99],iris$Petal.Length[99], pch = “*”, col = “red”) text(iris$Sepal.Length[99]+.35,iris$Petal.Length[99],”(5.1,3)”,col = “red”) points(iris$Sepal.Length[107],iris$Petal.Length[107],pch = “*”, col = “green”) text(iris$Sepal.Length[107],iris$Petal.Length[107]+.35,”(4.9,4.5)”,col = “green”)
“Petal Length”,main=
75)
Gabriel’s Plot

Sepal Length
1) Produce the exact same plot from above using ggplot as opposed to Base R graphics. That is, plot Petal Length versus Sepal Length split by Species. The colors of the points should be split according to Species. Also overlay three regression lines on the plot, one for each Species level. Make sure to include an appropriate legend and labels to the plot. Note: The function coef() extracts the intercept and the slope of an estimated line.

Part 2 (World’s Richest)
Background
Tasks
2) Open the file and make a new variable (dataframe) containing only the year, “P99”, “P99.5” and “P99.9” variables; these are the income levels which put someone at the 99th, 99.5th, and 99.9th, percentile of income. What was P99 in 1993? P99.5 in 1942? You must identify these using your code rather than looking up the values manually.
wtid <- read.csv(“wtid-report.csv”, as.is = TRUE)
### your code goes here
3) Plot the three percentile levels against time using ggplot. Make sure the axes are labeled appropriately, and in particular that the horizontal axis is labeled with years between 1913 and 2012, not just numbers from 1 to 100. Also make sure a legend is displayed that describes the multiple time series plot. Write one or two sentences describing how income inequality has changed throughout time. Remember library(ggplot2).