Displaying 1 result from an estimated 1 matches for "th3_b".
Did you mean:
th3_n
2012 Apr 19
1
non-numeric argument in mle2
Hi all,
I have some problems with the mle2 function
> RogersIIbinom <- function(N0,attackR3_B,Th3_B)
{N0-lambertW(attackR3_B*Th3_B*N0*exp(-attackR3_B*(24-Th3_B*N0)))/(attackR3_B*Th3_B)}
> RogersII_B <-
mle2(FR~dbinom(size=N0,prob=RogersIIbinom(N0,attackR3_B,Th3_B)/N0),start=list(attackR3_B=1.5,Th3_B=0.04),method="Nelder-Mead",data=dat)
Error in dbinom(x, size, prob, log) :
No...