请教各位大神。我现在要用stata做面板数据的负二项回归，patent2是因变量，role是自变量，而且是个取值为1,2,3,4的虚拟变量，还有几个控制变量，我想研究role对patent2的影响，要怎么做负二项回归分析？我自己琢磨了一下，进行了如下的尝试，不知道对不对？请各位多多指教！

1. 我先声明了截面变量和时间变量
. tsset firmid year
2. 然后定虚拟变量
. xi i.role
出来结果  i.role            _Irole_1-4          (naturally coded; _Irole_1 omitted)
3.假定使用随机效应模型
. xi: xtnbreg patent2 i.role density reach clustering partner_presample,re
出来结果

不知道这一过程对不对？
我想分析role对patent2的影响，那这一结果该如何分析？
急！请大家多给意见，谢谢啦~