我的数据是上市公司十年4921条研发投入数据,时间有gap,但能保证连续7年及以上有数据,想要做的是对每家公司进行分组回归,构建garch模型,求出研发投入变动趋势,出现错误:flat log likelihood encountered, cannot find uphill direction,我猜测可能原因是数据不够,但是已经发表的文章上用类似的数据也做出了结果。所以我想请问各位大佬有什么解决办法吗?这是我的代码:
- use c:\users\administrator\desktop\7year.dta,clear
- tsset company year
- egen g=group(company)
- sum g
- local n=r(max)
- gen f=.
- forv i=1/`n'{
- arch rd L.rd if g==`i',arch(1) garch(1)
- predict e if e(sample),r
- replace f=e if e(sample)
- drop e
- }
- drop g
复制代码
部分结果
- Iteration 344: log likelihood = -100.04287
- (switching optimization to BHHH)
- Iteration 345: log likelihood =-100.0428(not concave)
- Iteration 346: log likelihood = -100.04267(not concave)
- Iteration 347: log likelihood = -100.04267(not concave)
- Iteration 348: log likelihood = -100.04267(not concave)
- Iteration 349: log likelihood = -100.04267(not concave)
- (switching optimization to BFGS)
- BFGS stepping has contracted, resetting BFGS Hessian (87)
- Iteration 350: log likelihood = -100.04267
- Iteration 351: log likelihood = -100.04267(backed up)
- Iteration 352: log likelihood = -100.04266(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (88)
- Iteration 353: log likelihood =-100.0425
- Iteration 354: log likelihood = -100.04223(backed up)
- Iteration 355: log likelihood = -100.04223(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (89)
- Iteration 356: log likelihood = -100.04142
- Iteration 357: log likelihood = -100.04134(backed up)
- Iteration 358: log likelihood = -100.04133(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (90)
- Iteration 359: log likelihood = -100.04061
- (switching optimization to BHHH)
- Iteration 360: log likelihood = -100.04058(not concave)
- Iteration 361: log likelihood = -100.04051(not concave)
- Iteration 362: log likelihood = -100.04051(not concave)
- Iteration 363: log likelihood = -100.04051(not concave)
- Iteration 364: log likelihood = -100.04051(not concave)
- (switching optimization to BFGS)
- BFGS stepping has contracted, resetting BFGS Hessian (91)
- Iteration 365: log likelihood = -100.04051
- Iteration 366: log likelihood = -100.04051(backed up)
- Iteration 367: log likelihood = -100.04051(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (92)
- Iteration 368: log likelihood =-100.0403
- Iteration 369: log likelihood = -100.04014(backed up)
- Iteration 370: log likelihood = -100.04014(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (93)
- Iteration 371: log likelihood = -100.03937
- Iteration 372: log likelihood = -100.03935(backed up)
- Iteration 373: log likelihood = -100.03934(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (94)
- Iteration 374: log likelihood =-100.0388
- (switching optimization to BHHH)
- Iteration 375: log likelihood = -100.03878(not concave)
- Iteration 376: log likelihood = -100.03869(not concave)
- Iteration 377: log likelihood = -100.03869(not concave)
- Iteration 378: log likelihood = -100.03869(not concave)
- Iteration 379: log likelihood = -100.03869(not concave)
- (switching optimization to BFGS)
- BFGS stepping has contracted, resetting BFGS Hessian (95)
- Iteration 380: log likelihood = -100.03869
- Iteration 381: log likelihood = -100.03869(backed up)
- Iteration 382: log likelihood = -100.03869(backed up)
- Iteration 383: log likelihood = -100.03864
- BFGS stepping has contracted, resetting BFGS Hessian (96)
- Iteration 384: log likelihood = -100.03859
- Iteration 385: log likelihood = -100.03859(backed up)
- Iteration 386: log likelihood = -100.03857(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (97)
- Iteration 387: log likelihood = -100.03836
- Iteration 388: log likelihood = -100.03826(backed up)
- Iteration 389: log likelihood = -100.03826(backed up)
- (switching optimization to BHHH)
- Iteration 390: log likelihood = -100.03755(not concave)
- Iteration 391: log likelihood = -100.03748(not concave)
- Iteration 392: log likelihood = -100.03748(not concave)
- Iteration 393: log likelihood = -100.03748(not concave)
- Iteration 394: log likelihood = -100.03748(not concave)
- (switching optimization to BFGS)
- BFGS stepping has contracted, resetting BFGS Hessian (98)
- Iteration 395: log likelihood = -100.03748
- Iteration 396: log likelihood = -100.03748(backed up)
- Iteration 397: log likelihood = -100.03746(backed up)
- BFGS stepping has contracted, resetting BFGS Hessian (99)
- Iteration 398: log likelihood =-100.0372
- Iteration 399: log likelihood = -100.03715(backed up)
- Iteration 400: log likelihood = -100.03715(backed up)
- flat log likelihood encountered, cannot find uphill directi
- > on
- r(430);
复制代码