之前没用过stata，现在想用stata做多水平的poisson回归，包括固定效应和随机效应，各因素的大小用发生率比IRR表示，命令为 xi: xtmepoisson case sex i.age i.marriage i.education i.occupation || sheq :

    xtmepoisson,variance irr

①做多水平的poisson回归，我这个命令是正确的吗？

②模型的拟合优度检验要怎么看呢，或者要怎么做拟合优度检验？

③在做多水平的poisson回归之前，需要做各因素的单因素检验，同样用发生率比IRR衡量各因素的大小，要怎么做呢？