全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
4430 6
2010-08-07
小弟没什么matlab经验
如果我已经写出了likelihood function,怎么做mle estimate 呢?是不是直接最大化这个写出来得likelihood function 就行了,用什么函数最大化呢?我看到有个mle function,不知道可以用么?

另外,有关estimate regime switching with TIME-VARYING TRANSITION PROBABILITES ,哪位 同学有写过程序的,请发出来参考一下,万分感谢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-8-7 15:09:16
MLE就是个最优化过程而已 别想太复杂了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-8-7 16:55:01
哦,在符号运算里边,要有正态累计函数怎么办啊?其实我就是用mle求estimate,写出来的likelihood function 中有正态累计函数,matlab的 normcdf()的参数只能是real的
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-8-7 17:13:33
试试直接把normcdf写到方程里面
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-2-15 11:22:38
1# clark1025

用fminsearch做

先在script file里面猜测一组数据(对要估计的参数猜测)eg:guess=[1 2 3](代表你要估计的三个参数,1 2 3是我自己随便写的)
然后用 estimate=fminsearch('myfct',guess);

调用函数
function[f]=myfct(x)
global y T;   %把script里的sample设成global variable
T=length(y);
yita=x(1);
alpha1=x(2);
beta1=x(3);
然后设定你的log likelihood function
lnl=lnl=-(T-1)/2*log(2*pi)-slnh/2-sz/2;(随便举的例子)
然后
f=-lnl;(这个必须要,因为fminsearch是找最小值,加了负号才变成找最大值)
end;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-2-15 11:38:37
感谢5楼的朋友
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群