一直在使用Stata在复制中运行负二项式回归。不确定Stata是如何做到这一点的,但想知道是否有一个R功能/包做同样的事情? R将让我更好地了解它是如何工作的,因为我可以看到代码。
答案 0 :(得分:4)
查看glm.nb
包中的MASS
函数。如果您对“幕后”发生的事情感兴趣,只需在命令提示符下输入其名称即可查看该函数的源代码。
答案 1 :(得分:1)
如果你更习惯使用R,那么这可能就是你要走的路;但是,如果你对Stata中的“引擎盖”感兴趣,你可以通过使用
来看到与R中发生的事情大致相同的事情。set trace on
查看正在运行的代码(或从SSC下载tr
)或使用
viewsource xtnbreg.ado
查看xtnbreg
运行的实际代码。
如果您对Stata如何计算xtnbreg
中的结果感兴趣,可以在[XT]手册第367-370页中详细讨论似然函数,其中包含参考文献。