R SPSS GENLIN生成合并数据的Exp（B）和CI

是否有任何方法可以修改语法或利用R让SPSS GENLIN函数在多次插补后为合并数据生成exp（B）和95%CI？当前语法如下所示。语法中汇总后的（指数化）对每一个其他插补进行exp（B）和CI，但不对汇总数据进行处理。输出也显示了，但我裁剪了图像，只显示了30个插补中的1个

* Generalized Estimating Equations. 
GENLIN AnyMinorComp (REFERENCE=FIRST) BY Preop_Ambulation (ORDER=DESCENDING) WITH Age 
  /MODEL Preop_Ambulation Age INTERCEPT=YES 
 DISTRIBUTION=BINOMIAL LINK=LOGIT 
  /CRITERIA METHOD=FISHER(1) SCALE=1 MAXITERATIONS=100 MAXSTEPHALVING=5 PCONVERGE=1E-006(ABSOLUTE) 
    SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 LIKELIHOOD=FULL 
  /REPEATED SUBJECT=Study_ID WITHINSUBJECT=Site_Code SORT=YES CORRTYPE=INDEPENDENT ADJUSTCORR=YES 
    COVB=ROBUST MAXITERATIONS=100 PCONVERGE=1e-006(ABSOLUTE) UPDATECORR=1 
  /MISSING CLASSMISSING=EXCLUDE 
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).

试试这个：

DATASET DECLARE  temp.
OMS
  /SELECT TABLES
  /IF COMMANDS=['Generalized Linear Models'] SUBTYPES=['Parameter Estimates']
  /DESTINATION FORMAT=SAV NUMBERED=TableNumber_ OUTFILE='temp' VIEWER=YES.
GENLIN AnyMinorComp (REFERENCE=FIRST) BY Preop_Ambulation (ORDER=DESCENDING) WITH Age
  /MODEL Preop_Ambulation Age INTERCEPT=YES DISTRIBUTION=BINOMIAL LINK=LOGIT 
  /CRITERIA METHOD=FISHER(1) SCALE=1 MAXITERATIONS=100 MAXSTEPHALVING=5 PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 LIKELIHOOD=FULL
  /REPEATED SUBJECT=Study_ID WITHINSUBJECT=Site_Code SORT=YES CORRTYPE=INDEPENDENT ADJUSTCORR=YES 
    COVB=ROBUST MAXITERATIONS=100 PCONVERGE=1e-006(ABSOLUTE) UPDATECORR=1
  /MISSING CLASSMISSING=EXCLUDE
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).
OMSEND.
DATASET ACTIVATE temp.
COMPUTE ExpB=EXP(B).
COMPUTE Lower_A=EXP(Lower).
COMPUTE Upper_A=EXP(Upper).
EXECUTE.

---

这不会改变输出，但在新创建的数据集中，它将进行计算。

我已经查找了GENLIN函数，从我所知道的情况来看，它使用广义估计方程（GEE）对结构化数据建模（即最常处理混合效应模型的数据类型）。要执行gee，您可以使用

gee

包。至于多重插补，在R中有很多方法，最流行的是

mice

软件包，但是，也可以用随机森林、K近邻等进行插补。总之，你所要求的听起来绝对可行，您应该首先尝试自己编写一些代码，如果遇到卡滞，则会出现堆栈溢出。@AdamB。感谢您的回复-我已经完成了多重插补并进行了GEE；然而，无论出于何种原因，Exp（B）及其置信区间的计算均未针对合并数据进行。我想知道是否有比我更熟练的人知道我可以对语法进行修改，以便对汇总数据进行计算（计算为e^[B和下限/上限值]），我试图通过在查询末尾添加（指数化）来修改语法，但它只对插补和原始数据进行修改；不是池数据在这种情况下，池数据到底是什么？如果您只需要从拟合模型中获取指数斜率、上限CI和下限CI，那么您应该能够从模型摘要中提取这些值，然后使用

exp（）

函数将它们指数化。