Re: [R] function censReg in panel data setting

2011-09-24 Thread Arne Henningsen
On 24 September 2011 23:58, felipnunes 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 = logLikAttr, fnOrig

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 ini

Re: [R] function censReg in panel data setting

2011-09-14 Thread Arne Henningsen
On 14 September 2011 00:36, Arne Henningsen wrote: > Hi Igors > > On 13 September 2011 13:27, Igors 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 contr

Re: [R] function censReg in panel data setting

2011-09-13 Thread Arne Henningsen
Hi Igors On 13 September 2011 13:27, Igors 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 minus infinity in your mode

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 000

Re: [R] function censReg in panel data setting

2011-09-11 Thread Arne Henningsen
Dear Igors On 9 September 2011 22:00, Igors 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 try to use censReg on

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 a

Re: [R] function censReg in panel data setting

2011-09-09 Thread Arne Henningsen
On 8 September 2011 09:56, Igors 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 (argument "start"), it wo

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-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 lis

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: http://r.789695.n4.nabble.com/function-censReg-in-panel-d

Re: [R] function censReg in panel data setting

2011-09-07 Thread Arne Henningsen
Dear Igors On 7 September 2011 10:39, Igors 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 maxNRCompute(fn =

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, fnOr

Re: [R] function censReg in panel data setting

2011-09-06 Thread Arne Henningsen
On 6 September 2011 07:51, Igors 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 starting values: int

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 ind

Re: [R] function censReg in panel data setting

2011-09-05 Thread Arne Henningsen
Dear Igors On 5 September 2011 23:58, Igors 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 = logLikAttr,