X1 <- seq(0, 1, len=2000)
Pa1 <- c(0.295, 101.485, 0.705, 0.003762)
Y1 <- SarabiaE(P=Pa1, x=X1)
dev.new()
plot( X1, Y1, 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