在sas中做多水平logistic模型的结果，与在stata的不一样，如图。最后改用proc nlmixed做，又出现问题，改为什么估计方法比较好呢？在proc glimmix怎么对随机效应进行检验呢？请各位帮忙解答下。