vickey_chang 发表于 2014-10-11 19:00
似然函数不应该是这个吗？
n*alpha*log(beta)-n*log(gamma(alpha))+(alpha-1)*sum(log(dat))-beta*sum(d ...

nuomin 发表于 2014-10-11 21:36
sum函数有处理NaN类型的功能吗？

qoiqpwqr 发表于 2014-10-13 13:00
sum(x, na.rm=T)

vickey_chang 发表于 2014-10-14 10:34
dat

qoiqpwqr 发表于 2014-10-14 18:56
我们的程序师一样的，用我的程序也会出现警告信息。不过这应该是数值方面的问题，原因可能是去搜索负值的 ...

qoiqpwqr 发表于 2014-10-14 18:56
我们的程序师一样的，用我的程序也会出现警告信息。不过这应该是数值方面的问题，原因可能是去搜索负值的 ...

qoiqpwqr 发表于 2014-10-14 18:56
我们的程序师一样的，用我的程序也会出现警告信息。不过这应该是数值方面的问题，原因可能是去搜索负值的 ...

vickey_chang 发表于 2014-10-15 16:57
非常感谢帮助，最近在看统计模拟与R，很迷茫，编程真的好难啊

vickey_chang 发表于 2014-10-11 17:22
试着用optim求，但是总是提示“log(gamma(alpha))”错误，恳请各位大侠赐教！！！