data("green")
plot(OUTPUT~COST, data=green)
x <- seq(min(green$COST), max(green$COST), length.out=1001)
# dea
lines(x, dea_est(green$COST, green$OUTPUT, x, type="dea"),
lty=1, col="red")
# fdh
lines(x, dea_est(green$COST, green$OUTPUT, x, type="fdh"),
lty=2, col="blue")
# lfdh
lines(x, dea_est(green$COST, green$OUTPUT, x, type="lfdh"),
lty=3, col="green")
legend("topleft",legend=c("dea","fdh","lfdh"), col=c("red","blue","green"),
lty=1:3)Run the code above in your browser using DataLab