plot(rho - normal ~ rho, data = simulation, type = "l", col = "blue",
ylim = c(0, 0.016),
ylab = expression(rho - bar(rho)[sim]),
xlab = expression(rho))
points(rho - logistic ~ rho, data = simulation, type = "l", col = "red")
points(rho - gumbel ~ rho, data = simulation, type = "l", col = "grey")
legend("topright", legend = c("Normal", "Logistic", "Gumbel"),
col = c("blue", "red", "grey"), pch = "l" )
Run the code above in your browser using DataLab