Check derivations of real77 and real77r
dpritchard opened this issue · comments
Daniel Pritchard commented
We need to check the derivations / behaviour of real77 and real77r.
Some files that might help are attached.
Daniel Pritchard commented
I think the best solutions might be:
frair_real77_alt3 <- function(X, b, q, h, T) {
a <- (b*X^q)
return((a*X*T)/(1+a*X*h))
}
frair_real77r_alt3 <- function(X, b, q, h, T) {
a <- (b*X^q)
return(X-lamW::lambertW0(a*h*X*exp(-a*(T-h*X)))/(a*h))
}
Daniel Pritchard commented
Daniel Pritchard commented
The solutions mentioned above have now been implemented as flexp
and flexpnr
, respectively.