Displaying 1 result from an estimated 1 matches for "funcl".
Did you mean:
func
2012 Jan 03
0
optim - Self-Start values - growth function
...3271.703374, m= -1.870935)) ; start values
obtained: A from data, k & m substituting. But I obtain this error:
/Error en numericDeriv(form[[3L]], names(ind), env) :
Missing value or an infinity produced when evaluating the model /
b) est.o4 <- optim( c(30.3656, 2.841345, 3270.109), funcL, method="SANN",
hessian=TRUE, yr= data$Yia, x=data$Yi)
where funcL <- function( par, yr, x ) {
A = par[1]
m = par[2]
k = par[3]
y <- A*exp(-k* (1/(-k/(log(x/A)))^(1/m)+a)^m)
sum((yr-y)^2) }
Output is as follows:
/$par
[1] 46.9067228 -0.7824855 16.5397317
$value
[1] 977.2446...