Simple Density plots

#create datasets
exceldata<- read.csv("D:/Jeremiah/Count of species per village.csv")
data<- data.frame(exceldata)
funestus_prop <-data$funestus_prop
gambiae_prop <-data$gambiae_prop

#plot first kernel density plot
kd1 <- density(funestus_prop)
plot(kd1, col='blue', lwd=2)

#plot second kernel density plot
kd2 <- density(gambiae_prop)
plot(kd2, col='red', lwd=2)

#plot combined kernel density plot
kd1 <- density(funestus_prop)
plot(kd1, col='blue', lwd=2, main='Combined Kernel Density Plot')
kd2 <- density(gambiae_prop)
lines(kd2, col='red', lwd=2)