似然函数为∑∑hik[(1-uik) log(1-pik)+uiklog(pik)],其中\[h_{ik} , \mu_{ik}\]已知.p满足\[logit(p_{ik})=e^{x_{i}^{'}\beta+m_{k}}\],\[x_{i}\]已知,要求\[\beta,m_{k}\]的极大似然估计。
楼主尝试了R中的ucminf,optim,maxLik这些直接优化的函数,但都出现了warning,说是在迭代过程中出现了NA值。不知道大家有没有什么更好的办法?不胜感激~~~~
|
楼主: ustc_challenge
|
1596
0
[统计软件] 用R做极大似然估计的问题 |
jg-xs1京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明


