Re: [R] function censReg in panel data setting

2011-09-24 Thread felipnunes
Hi guys, did you figure out a way to solve your problem, Igor? I'm still trying to understand what's happening to my estimations. My data set is also big (50,000) and I get the same error: Error in maxNRCompute(fn = logLikAttr, fnOrig = fn, gradOrig = grad, hessOrig = hess, : NA in the

Re: [R] function censReg in panel data setting

2011-09-24 Thread Arne Henningsen
On 24 September 2011 23:58, felipnunes felipnu...@gmail.com wrote: Hi guys, did you figure out a way to solve your problem, Igor? I'm still trying to understand what's happening to my estimations. My data set is also big (50,000) and I get the same error: Error in maxNRCompute(fn =

Re: [R] function censReg in panel data setting

2011-09-14 Thread Arne Henningsen
On 14 September 2011 00:36, Arne Henningsen arne.henning...@googlemail.com wrote: Hi Igors On 13 September 2011 13:27, Igors igors.lahanc...@gmail.com wrote: Any success in finding possible solutions for my problem? Somewhat. The calculation of the log-likelihood values is numerically much

Re: [R] function censReg in panel data setting

2011-09-13 Thread Igors
Hi Arne, Any success in finding possible solutions for my problem? I have tried to experiment with size of sample and I get really bad picture. I can't get it work even if sample is ~ 1000 obs. And it is way less than I would like to see working, taking into account my full sample size ~ 540

Re: [R] function censReg in panel data setting

2011-09-13 Thread Arne Henningsen
Hi Igors On 13 September 2011 13:27, Igors igors.lahanc...@gmail.com wrote: Any success in finding possible solutions for my problem? Somewhat. The calculation of the log-likelihood values is numerically much more robust/stable now. The log-likelihood contributions of some individuals became

Re: [R] function censReg in panel data setting

2011-09-11 Thread Arne Henningsen
Dear Igors On 9 September 2011 22:00, Igors igors.lahanc...@gmail.com wrote: Thank you for your answer. However it doesn't solve my problem fully. The problem is that I have much bigger data set than I sent to you (it was only a small part : 874 obs.). My full data set is 546718 obs. If I

Re: [R] function censReg in panel data setting

2011-09-09 Thread Igors
Dear Arne, Does censReg support unbalanced panel data? /Igors -- View this message in context: http://r.789695.n4.nabble.com/function-censReg-in-panel-data-setting-tp3792227p3801030.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] function censReg in panel data setting

2011-09-09 Thread Arne Henningsen
On 8 September 2011 09:56, Igors igors.lahanc...@gmail.com wrote: Does censReg expect from panel data to be balanced? No. censReg() can estimate models with unbalanced panel data. The estimation in the code that you sent me indeed does not work but if I remove the user-specified starting values

Re: [R] function censReg in panel data setting

2011-09-09 Thread Igors
Dear Arne, Thank you for your answer. However it doesn't solve my problem fully. The problem is that I have much bigger data set than I sent to you (it was only a small part : 874 obs.). My full data set is 546718 obs. If I try to use censReg on full data set, then it still gives me the same

Re: [R] function censReg in panel data setting

2011-09-08 Thread Igors
Does censReg expect from panel data to be balanced? Because in my case it is unbalanced. Could this be a reason for errors? Best, Igors -- View this message in context: http://r.789695.n4.nabble.com/function-censReg-in-panel-data-setting-tp3792227p3798113.html Sent from the R help mailing

Re: [R] function censReg in panel data setting

2011-09-07 Thread Igors
Dear Arne, Thank you for fixing the package. However I am still struggling to obtain model estmates. The same code: UpC - censReg(Power ~ Windspeed, left = -Inf, right = 2000,data=PData_In,method=BHHH,nGHQ = 4,start=c(-691.18,186.79,3.9,3.9)) Error in maxNRCompute(fn = logLikAttr, fnOrig =

Re: [R] function censReg in panel data setting

2011-09-07 Thread Arne Henningsen
Dear Igors On 7 September 2011 10:39, Igors igors.lahanc...@gmail.com wrote: However I am still struggling to obtain model estmates. The same code: UpC - censReg(Power ~ Windspeed, left = -Inf, right = 2000,data=PData_In,method=BHHH,nGHQ = 4,start=c(-691.18,186.79,3.9,3.9)) Error in

Re: [R] function censReg in panel data setting

2011-09-07 Thread Igors
Dear Arne, I have sent you an e-mail to your e-mail at gmail.com. There I have attached my data set and the code for this particular problem. I hope to hear from you soon. Many thanks in advance! Igors -- View this message in context:

Re: [R] function censReg in panel data setting

2011-09-06 Thread Igors
You chose a suitable solution for the first problem (NA in initial gradient). Unfortunately, the documentation of censReg is not very clear regarding starting values of panel tobit models. Please note that you have to specify 4 starting values: intercept, slope parameter, variance of the

Re: [R] function censReg in panel data setting

2011-09-06 Thread Arne Henningsen
On 6 September 2011 07:51, Igors igors.lahanc...@gmail.com wrote: You chose a suitable solution for the first problem (NA in initial gradient). Unfortunately, the documentation of censReg is not very clear regarding starting values of panel tobit models. Please note that you have to specify 4

[R] function censReg in panel data setting

2011-09-05 Thread Igors
Hello all, I have a problem estimating Random Effects model using censReg function. small part of code: UpC - censReg(Power ~ Windspeed, left = -Inf, right = 2000,data=PData_In,method=BHHH,nGHQ = 4) Error in maxNRCompute(fn = logLikAttr, fnOrig = fn, gradOrig = grad, hessOrig = hess, : NA

Re: [R] function censReg in panel data setting

2011-09-05 Thread Arne Henningsen
Dear Igors On 5 September 2011 23:58, Igors igors.lahanc...@gmail.com wrote: I have a problem estimating Random Effects model using censReg function. small part of code: UpC - censReg(Power ~ Windspeed, left = -Inf, right = 2000,data=PData_In,method=BHHH,nGHQ = 4) Error in maxNRCompute(fn