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) : 
  Non-numeric argument to mathematical function

Can somenone explain met what this error means? All my parameters and data 
are defined, so I don't understand what the non-numeric argument is??? 
With different equations the script does work...........

For this function to run you need packages: bbmle and emdbook.
Testdata I pasted here

N0      FR
5       4
3       2
5       3
5       4
6       4
5       4
4       2
5       4
6       5
5       3
10      5
14      7
12      6
17      9
10      4
16      8
15      8
14      5
15      6
15      5
19      10
20      11
22      10
21      12
25      12
26      11
25      10
25      12
25      11
24      9
35      12
33      12
35      13
36      12
31      12
37      12
39      13
36      14
35      12
35      10
46      11
44      10
41      12
43      14
45      12
48      11
45      13
45      9
49      12
45      11
56      12
59      15
52      11
54      10
51      12
55      12
61      13
55      12
54      10
55      12

Kind regards,
Met vriendelijke groeten,
Joachim

Don't waste paper! Think about the environment before printing this e-mail

______________________________________

Joachim Audenaert
Adviesdienst Gewasbescherming
Proefcentrum voor Sierteelt
Schaessestraat 18
B-9070 Destelbergen
Belgium
Tel. +32 9 353 94 71
Fax +32 9 353 94 95
E-mail: joachim.audena...@pcsierteelt.be
www.pcsierteelt.be
______________________________________
        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to