似然函数为∑∑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值。不知道大家有没有什么更好的办法？不胜感激~~~~