首页 理论教育 结局变量分析及SAS实现

结局变量分析及SAS实现

时间:2023-10-30 理论教育 版权反馈
【摘要】:10.4.1.3LTB法LTB法最早由Lanza,Tan,&Bray提出,后进一步发展适用于连续型、二分类、多分类、计数类型的结局变量。

结局变量分析及SAS实现

10.4.1.1 简单三步法

简单三步法虽然过程简单,但由于未处理分类误差,容易低估潜类别和结局变量之间的关系。但简单三步法的优势是在探讨潜类别对因变量影响时,可以控制其他协变量。

SAS程序:

SAS结果:

图10-17 多重线性回归

SAS结果解释:

图10-17表明,在控制其他协变量后,男大学生抑郁水平均值较低(β=-0.4212,p<0.001),与低风险组相比,高风险组(β=0.6223,p<0.001)的大学生抑郁水平均值较高。

10.4.1.2 BCH法

BCH法最早由Block,Croon&Hagenaars(2004)提出,只适用于分类结局变量,随后Vermunt(2010)和Vermunt&Magidson(2015)对该方法进行了改善,亦可应用于连续型结局变量。第一步,在不纳入结局变量的情况下进行LCA。第二步,依据后验概率进行加权,主要方法有非调整的莫代尔分配法(unadjusted modal assignment)、BCH-调整的莫代尔分配法(BCH-adjusted modal assignment)、非调整的比例分配法(unadjusted proportional assignment)、BCH-调整的比例分配法(BCH-adjusted proportional assignment)。第三步,根据权重变量计算结局变量在每个类别内的均值(连续型结局变量)或概率(分类结局变量)。BCH每次只纳入一个结局变量,如果对其他结果变量感兴趣,需要单独运行其他模型(Dziak et al,2017)。

SAS程序:

%INCLUDE

"C:\Users\Administrator\Desktop\LCA\SAS\LCA_Distal_BCH_v110-1-10jyv0q\LCA_Distal_BCH_v110.sas";

proc lca data=exe 10_4 outparam=conti_param outpost=conti_post;

id ID;

nclass 3;

items ACE1-ACE13;

categories 2 2 2 2 2 2 2 2 2 2 2 2 2;

seed 941622;

RHO prior=1;

nstarts 20;

maxiter 5000;

criterion 0.000001;

run;

%LCA_Distal_BCH(input_data=a,

param=conti_param,

post=conti_post,

id=ID,

distal=ZC,

metric=Continuous);

SAS程序解释:

首先,从https://www.methodology.psu.edu/downloads/distal/下载并调用LCA_Distal_BCH_v110.sas宏。通过maxiter指定最大迭代次数,distal指定结局变量名,metric指定结局变量的类型,可以是二分类变量(binary),连续型变量(continuous),计数变量(count),分类变量(categorical),在本例中为连续型变量。

SAS结果:

图10-18 BCH估计的抑郁水平的均值、标准误、95% CI

(www.xing528.com)

图10-19 Wald卡方检验

SAS结果解释:

图10-18展示了BCH估计的不同类别的大学生抑郁水平的均值、标准误、95% CI。图10-19通过Wald卡方检验比较了不同类别的大学生抑郁水均值的差异。具体来说,与高风险组相比,低风险组和暴力组中的大学生抑郁水平均值较低,差异具有统计学意义。

10.4.1.3 LTB法

LTB法最早由Lanza,Tan,&Bray(2013)提出,后进一步发展适用于连续型、二分类、多分类、计数类型的结局变量。当连续型结局变量的方差在不同类别内相等时,LTB法的估计是无偏的。但当存在异方差时,BCH法比LTB法更稳健、偏倚更小(Bakk&Vermunt,2016)。目前,LTB实际应用较少,主要用于和BCH法的结果进行比较(Dziak,et al.,2017)。从https://www.methodology.psu.edu/downloads/distal/下载并调用LCA_Distal_LTB_v110.sas宏。

SAS程序:

%INCLUDE

"C:\Users\Administrator\Desktop\LCA\SAS\LCA_Distal_LTB_v110-1-spex7u\LCA_Distal_LTB_v110.sas";

pro lcadata=exe 10_4 outparam=Conti_param;

id ID;

nclass 3;

items ACE1-ACE13;

categories 2 2 2 2 2 2 2 2 2 2 2 2 2;

covariates ZC;

seed 941622;

RHO prior=1;

nstarts 20;

maxiter 5000;

criterion 0.000001;

run;

%LCA_Distal_LTB(input_data=a,

param=Conti_param,

distal=ZC,

metric=2,

output_data set_name=Conti_res);

SAS程序解释:

首先将结局变量作为协变量纳入LCA,通过metric指定结局变量类型,metric=1表示二分类结局变量,metric=2连续型结局变量,metric=3表示计数结局变量,metric=4表示分类结局变量,在本例中结局变量为连续型变量。

SAS结果:

图10-20 β系数估计值(标准误)

图10-21 OR估计值(95% CI)

图10-22 β系数的p值

SAS结果解释:

图10-20、图10-21和图10-22分别展示了β系数估计值(标准误)、OR值(95% CI)及p值。与暴力组相比,高风险组中的大学生抑郁水平均值更高(β=0.5841,OR=1.7934,p<0.001)。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈