s****e
发帖数: 1180 | 1 有一个R程序，是用EM算法估计generalized linear model中的参数的。程序没有任何
语法问题，但就是运行的过程中总有不能求倒数的数出现。程序的算法，数据，我都有
。不知哪位大大能帮我看一下，帮我解决这个问题，我愿付报偿。站内联系 .多谢！ |
c*****m
发帖数: 4817 | 2 不能求倒數？是有0還是oo？
你這個到底是R問題，還是算法本身問題啊？ |
S******y
发帖数: 1123 | 3 Debug in R -
1) build some fake test data, which you know well
2) run your R program line by line interactively
3) stop where the error comes up
Then post your intermediate data frame and those several lines of R code
here as well as your R error message here. |
s****e
发帖数: 1180 | 4 我修改了程序，现在run起来没任何问题，但就是估计的不精确。如我以前所说，这个R
程序是从SAS/IML翻译过来的，SAS/IML的那个程序对EM算法的大部分初始值不work，但
对某些初始值work，而且能估计出main effect，但R程序几乎对所有初始值都work，但
就是估计不精确。换句话说，这是一个shrink方法，SAS/IML程序对某些初始值，可以
shrink不显著的effect，并估计出显著的main effect。但R程序虽对任何初始值都work
，但对任何初始值，都没shrink效应。哪位牛人能解决这个问题，本人愿付报偿。多谢！
【在 S******y 的大作中提到】
: Debug in R -
: 1) build some fake test data, which you know well
: 2) run your R program line by line interactively
: 3) stop where the error comes up
: Then post your intermediate data frame and those several lines of R code
: here as well as your R error message here.
|
s****e
发帖数: 1180 | 5 Ding!
【在 s****e 的大作中提到】
: 有一个R程序，是用EM算法估计generalized linear model中的参数的。程序没有任何
: 语法问题，但就是运行的过程中总有不能求倒数的数出现。程序的算法，数据，我都有
: 。不知哪位大大能帮我看一下，帮我解决这个问题，我愿付报偿。站内联系 .多谢！
|
p********a
发帖数: 5352 | |
s*r
发帖数: 2757 | 7 这个我做不了
step by step R debug本来就是很头疼的，更何况我对EM shrinkage不熟，sas里面似
乎没有什么step by step debug的工具
【在 p********a 的大作中提到】
: SIR同学呢？挣外快的机会来了！
|
D******n
发帖数: 2836 | 8 use print...hehe,debuging in R is painful.
【在 s*r 的大作中提到】
: 这个我做不了
: step by step R debug本来就是很头疼的，更何况我对EM shrinkage不熟，sas里面似
: 乎没有什么step by step debug的工具
|
s****e
发帖数: 1180 | 9 我昨天又调整了程序，得到还不错的结果。这下，我更晕，为什么同样的程序，用R写
，能得到结果，用SAS/IML写就不行。当然我还有其它数据，有不会的地方再上来请教
大家。谢谢大家的关心！ |
