Re: [R] Maximum likelihood estimation

2008-09-04 Thread toh
. > > I think you intended to optimize over c(a,b,p, lambda), so you need to > specify them as a single vector. > > You may be making life unnecessarily hard for yourself: see function mle() > in package stats4. > > Showing your code without a verbal description

[R] Maximum likelihood estimation

2008-09-04 Thread toh
Hi R-experts, I'm new to R in mle. I tried to do the following but just couldn't get it right. Hope someone can point out the mistakes. thanks a lot. t <- c(1:90) y <- c(5,10,15,20,26,34,36,43,47,49,80,84,108,157,171,183,191,200,204,211,217,226,230, 234,236,240,243,252,254,259,263,264,268,271,27

[R] How to write the log function in R

2008-09-01 Thread toh
Hi, I need help to write the following log-likelihood function in R: log(L) = sum [(y_k - y_k-1).log(m(t_k) - m(t_k-1)) - (m(t_k) - m(t_k-1)) - log((y_k -y_k-1)!)] I need to write this log-likelihood function in order to find the parameters by MLE method. Appreciate the help. -- View this m