s****e
发帖数: 1180 | 1 谁能帮我调试好一个程序，我愿付给一定报偿。写了一个程序，好长时间了，但就是
run了之后有问题，不收敛，error越来越大。以前也碰到过这种情况，没放在心上，现
在又一次遇到，下决心克服，但还是无奈。不知论坛上有何牛人，能帮一下，我愿付一
定报偿。算法，程序，程序run出的结果我都有。这是一个用EM algorithm估计
generalized linear mixed model里的参数的问题。当参数个数大约在500左右的情况
时，已经验证过了，行得通。现在，加上了interaction，参数个数在10000左右。我用
SAS/IML(which is similar to R, if you can write a program in R, it is also
OK)写了个程序，error总是无限增大，无法收敛。
重新写也行。只要把问题解决了就行。我的意思是，我的程序起码能给大家一些参考。
如果你能证明这个算法在这个问题上行不通啥的，也行。也算解决了问题。不过算法是
别人发表过的。
站内联系。多谢！ | s****e
发帖数: 1180 | 2 Ding!
【在 s****e 的大作中提到】
: 谁能帮我调试好一个程序，我愿付给一定报偿。写了一个程序，好长时间了，但就是
: run了之后有问题，不收敛，error越来越大。以前也碰到过这种情况，没放在心上，现
: 在又一次遇到，下决心克服，但还是无奈。不知论坛上有何牛人，能帮一下，我愿付一
: 定报偿。算法，程序，程序run出的结果我都有。这是一个用EM algorithm估计
: generalized linear mixed model里的参数的问题。当参数个数大约在500左右的情况
: 时，已经验证过了，行得通。现在，加上了interaction，参数个数在10000左右。我用
: SAS/IML(which is similar to R, if you can write a program in R, it is also
: OK)写了个程序，error总是无限增大，无法收敛。
: 重新写也行。只要把问题解决了就行。我的意思是，我的程序起码能给大家一些参考。
: 如果你能证明这个算法在这个问题上行不通啥的，也行。也算解决了问题。不过算法是
|
|
