请问老师，我有一组资料，是心律失常患者经过两种不同的处理策略（手术和药物）后随访1月，3月，6月，12月后早搏的变化资料。 Treatment_Group 1=药物; Treatment_Group 2= 手术。Visit1=处理前的baseline随访，2=处理后1月随访，3=处理后3月随访，4=处理后6月随访，5=处理后 12月随访。因为心律失常的分布为非正态分布，审稿专家要求用Poisson GEE model来做。并要求将随访Visit 和reatment_Group 作为预测变量。
我的数据输入格式如下：

ID	Treatment_Group	Visit	arrhythmia_Burden
1	1	1	29
1	1	2	9.00
1	1	3	9.00
1	1	4	9.00
1	1	5	8.00
2	1	1	22
2	1	2	8.00
2	1	3	8.00
2	1	4	9.00
2	1	5	8.00
3	2	1	12.00
3	2	2	0.24
3	2	3	0.14
3	2	4	0.09
3	2	5	0.21
4	2	1	12.00
4	2	2	0.21
4	2	3	0.13
4	2	4	0.14
4	2	5	0.03

程序如下
xi: xtgee arrhythmia_Burden  i.Visit   Treatment_Group, family(poisson) link(log) corr(exchangeable)

结果显示
-----------------------------------------------------------------------------------------
arrthymia_burden |      Coef.      Std. Err.        z       P>|z|     [95% Conf. Interval]
-------------+--------------------------------------------------------------------------------
             Vist~2  |  -1.593604   .0227498   -70.05    0.000    -1.638192   -1.549015
             Vist~3  |  -1.577053   .0225755   -69.86    0.000    -1.621301   -1.532806
             Vist~6  |  -1.643281      .0235     -69.93    0.000     -1.68934   -1.597221
           Vist~12  |  -1.455204    .0210762   -69.05   0.000    -1.496513   -1.41389
      treatment_gr |   .1543955    .015338    10.07     0.000     .1243335    .1844575
               _cons |    2.72943    .0166416   164.01    0.000     2.696813    2.762047
------------------------------------------------------------------------------------------------
请问这样处理正确吗？审稿专家要求还要先做“ an interaction term between the predictorvariables should be examined”，请问这个怎么做呢？谢谢！