# NOT RUN {
data(matFish)
gonad_mat = gonad_mature(matFish, varNames = c("total_length", "stage_mat"), inmName = "I",
matName = c("II", "III", "IV"), method = "fq", niter = 50)
plot(gonad_mat, xlab = "Total length (cm.)", ylab = "Proportion mature", col = c("blue", "red"))
# }
Run the code above in your browser using DataLab