| 所在主题: | |
| 文件名: luetkepohl(2007) .pdf | |
| 资料下载链接地址: https://bbs.pinggu.org/a-1544879.html | |
| 附件大小: | |
|
我想用matlab中的fminunc命令做最大似然估计,请问我这样写的错误在哪里?
mylog_lik.m里的code是: function val=mylog_lik(B,omega) t = 200; T = 280; cov1 = [5.2486 1.9729 0.7363;1.9729 3.0687 0.9202;0.7363 0.9202 5.3797]; cov2 = [37.7688 8.0768 8.5434;8.0768 6.4612 6.9920;8.5434 6.9920 38.9579]; val=(t-1)*(log(det(B*B'))+trace(inv(cov1*(B*B'))))/2+(T-t+1)*(log(det(B*omega*B'))+trace(cov2*(B*omega*B')^(-1)))/2; maxlike.m里是 B0=[0.1 0.1 0.1;0.1 0.1 0.1;0.1 0.1 0.1]; omega0 = [0.1 0 0 ; 0 0.1 0; 0 0 0.1]; x0 = [B0,omega0]; [x,fval,exitflag,output,grad,hessian]=fminunc('mylog_lik',x0); 初始值应该根据什么设置呢?初始值应该是向量吗? |
|
熟悉论坛请点击新手指南
|
|
| 下载说明 | |
|
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。 2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。 3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明