X1 <- seq(0, 1, len=2000)
Pa3 <- c(0, 1, 0.446, 1.739)
Y3 <- SHE(P=Pa3, x=X1)
dev.new()
plot( X1, Y3, cex.lab=1.5, cex.axis=1.5, type="l", asp=1, xaxs="i",
yaxs="i", xlim=c(0, 1), ylim=c(0, 1),
xlab="Cumulative proportion of the number of infructescences",
ylab="Cumulative proportion of the infructescence length" )
graphics.off()
Run the code above in your browser using DataLab