[R] Problems in programming a simple likelihood

2007-04-18 Thread Deepankar Basu
As part of carrying out a complicated maximum likelihood estimation, I am trying to learn to program likelihoods in R. I started with a simple probit model but am unable to get the code to work. Any help or suggestions are most welcome. I give my code below: ml

Re: [R] Problems in programming a simple likelihood

2007-04-18 Thread Stephen Weigand
Deepankar, Some general advice from a non-expert: Write your likelihoods without a for loop. This is important because the likelihood is evaluated multiple times in the maximization process and you don't want to be looping looping looping ... Always try multiple starting values Sometimes it hel

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Dimitris Rizopoulos
ium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://med.kuleuven.be/biostat/ http://www.student.kuleuven.be/~m0390867/dimitris.htm - Original Message - From: "Deepankar Basu" <[EMAIL PROTECTED]> To: Sent: Thursday, April 19, 2007 12:38 AM Subject: [R]

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Deepankar Basu
t; > Dimitris Rizopoulos > Ph.D. Student > Biostatistical Centre > School of Public Health > Catholic University of Leuven > > Address: Kapucijnenvoer 35, Leuven, Belgium > Tel: +32/(0)16/336899 > Fax: +32/(0)16/337015 > Web: http://med.kuleuven.be/biostat/ > http:

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Ravi Varadhan
> I hope it helps. > > Best, > Dimitris > > > Dimitris Rizopoulos > Ph.D. Student > Biostatistical Centre > School of Public Health > Catholic University of Leuven > > Address: Kapucijnenvoer 35, Leuven, Belgium > Tel: +32/(0)16/336899 > Fax: +32

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Deepankar Basu
> > Ph: (410) 502-2619 > > Fax: (410) 614-9625 > > Email: [EMAIL PROTECTED] > > Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html > > > > > > > -Original Message----- > From: [EMAIL PROTECTED] > [mailt