Search Archives

By the way, in my last post I forgot to return negative of llfn, hence the llfn will be as follows: llfn <- function(param) { bx <- param[1] b1 <- param[2] b2 <- param[3] b3 <- param[4] b <- param[5] lL1...

Of course, that was the trick! It works now. Thank you very much Rui, I am very grateful. I hope this thread will help others as well. Best, Rui Barradas wrote > > Once again, sorry. > I had a different llfn in...

Hello again, You are absolutely right about probabilities.. Thanks for reminding me about that. I did exactly how you said but in the end I receive the error : "objective function in optim evaluates to length 12 not 1". I checked...

Hello John, Thanks for your response. Can you explain why it is better to use NM or nmk instead of optim? By the way, my original model consists 24 parameters. In this case is it still better to use NM...

Hello Rui, First of all, thanks a lot! 1. I changed par to param, 2. t is a variable too, a binary one, b is the parameter associated to it, 4. Yes, this is where I am stuck actually. I...

Hello, I am new to R and I am trying to estimate a discrete model with three choices. I am stuck at a point and cannot find a solution. I have probability functions for occurrence of these choices, and then...

Hello again, I changed the name to tt. and for a and tt actually I was getting them from data, I didn't put them here in the question. Now I restructured my code and below I copy the full...