Set maturity parameter bound to 1000 for large fish

quang-huynh · quang-huynh · commit ea7bb29d25ed · 2023-09-07T14:16:33.000-07:00
diff --git a/R/maturity.R b/R/maturity.R
@@ -633,9 +633,9 @@ extract_maturity_perc_re <- function(betas, random_intercepts, model) {
   out
 }
 
-binomial_perc <- function(x, a, b, perc = 0.5, linkinv, ...) {
+binomial_perc <- function(a, b, perc = 0.5, linkinv, ...) {
   f <- function(x) linkinv(a + b * x) - perc
-  uniroot(f, interval = c(-100, 100))$root
+  uniroot(f, interval = c(-100, 1000))$root
 }
 
 mat_par_delta_method <- function(model, perc = 0.5) {

Original file line number	Diff line number	Diff line change
`@@ -633,9 +633,9 @@ extract_maturity_perc_re <- function(betas, random_intercepts, model) {`
`633`	`633`	`out`
`634`	`634`	`}`
`635`	`635`
`636`		`-binomial_perc <- function(x, a, b, perc = 0.5, linkinv, ...) {`
	`636`	`+binomial_perc <- function(a, b, perc = 0.5, linkinv, ...) {`
`637`	`637`	`f <- function(x) linkinv(a + b * x) - perc`
`638`		`- uniroot(f, interval = c(-100, 100))$root`
	`638`	`+ uniroot(f, interval = c(-100, 1000))$root`
`639`	`639`	`}`
`640`	`640`
`641`	`641`	`mat_par_delta_method <- function(model, perc = 0.5) {`