# NOT RUN {
data(crabdata)
classify_data = classify_mature(crabdata, varNames = c("carapace_width", "chela_height"),
varSex = "sex_category", selectSex = NULL, method = "ld")
my_mature = morph_mature(classify_data, method = "fq", niter = 50)
plot(my_mature, xlab = "Carapace width (mm.)", ylab = "Proportion mature", col = c("blue", "red"))
# }
Run the code above in your browser using DataLab