各位大神好,我的研究是利用风险调整模型对30个医疗机构的某个率(Y为二分类变量)进行风险调整,风险调整模型包括了各种社会学以及临床学的病人相关因素和医院层面的一些变量,最后对风险调整后的该率进行比较。风险调整率的公式是单个机构实际发生Y的人数/该机构预测发生Y的人数*总体率。在模型选择上,我试了二元logistic回归和多水平logistic回归模型。软件是stata14.0。
我发现,在做二元logistic回归模型的时候(只是用病人的社会学和临床学因素)得到的预测的人数是比较理想的,即各机构的预测人数在调整前和调整都是有明显差异的。然后做多水平logistic回归模型,病人层面为一水平,医院层面为二水平。结果发现,预测的发生Y的人数和实际发生Y的人数十分接近。因此得到每一个机构调整后的率与总体率十分接近,而且各个机构之间差别不大。这与我之前查过相关文献是有出入的,文献报道的风险调整各机构的率的分布情况应该是在总体率上下波动,但我不明白为什么会出现这样的结果。请问各位大神,这样的情况是哪个环节出现了问题,希望各位大神不吝赐教,不胜感激!!