原文链接：http://tecdat.cn/?p=12761

原文出处：拓端数据部落公众号

我们使用整容手术数据说明两种中心化类型。将此文件加载到SPSS中。假设我们要中心化的变量BDI。

数据中心化

首先，我们需要找出BDI的平均得分。我们可以使用一些简单的描述性统计信息。选择进入对话框。选择BDI并将其拖到标有Variable（s）的框中，然后单击并仅选择均值。

结果输出告诉我们平均值为23.05：

我们使用此值将变量中心化。通过选择访问计算命令。在出现的对话框中，在标有“目标变量”的框中输入名称BDI_Centred，然后单击并为变量指定一个更具描述性的名称。选择变量BDI并将其拖动到标记为Numeric Expression的区域，然后单击，然后键入平均值（23.05）。完成的对话框如图所示。

单击，将创建一个名为BDI_Centred的新变量，该变量以BDI的平均值为中心。这个新变量的均值应约为0：运行一些描述性统计数据。

可以在语法窗口中通过输入以下内容执行相同的操作：

COMPUTE BDI_Centred = BDI−23.05.
EXECUTE.

组均值中心化

组均值中心化要复杂得多。第一步是创建一个包含组均值的文件。让我们再试一次以获取BDI分数。我们希望将此变量在Clinic的2级变量中中心化。我们首先需要知道每个组中的平均BDI，并以SPSS形式保存该信息。我们要使用aggregate命令。在此对话框（图3）中，我们要选择Clinic并将其拖动到标有Break Variable（s）的区域。这说明着将使用变量Clinic来分割数据文件（换句话说，当计算平均值时，它将对每个诊所分别进行处理）。然后，我们需要选择BDI并将其拖动到标记为变量汇总的区域。一旦选择了此变量，默认值就是SPSS将创建一个名为BDI_mean的新变量，这是BDI的平均值（显然是由Clinic分割）。我们需要将此信息保存在一个文件中，以便以后使用。

默认情况下，SPSS会将名称为aggr.sav的文件保存在默认目录中。如果您想将其保存在其他位置或使用其他名称，则单击以打开一个普通的文件系统对话框，可以在其中命名文件并导航至要保存在其中的目录。单击创建此新文件。

如果打开生成的数据文件，则会看到它仅包含两列，其中一列带有一个数字，用于指定数据来自的诊所（共有10个诊所），第二个包含每个诊所内的平均BDI得分。

当SPSS创建汇总数据文件时，它将按从最低到最高的顺序对诊所进行排序（无论它们在数据集中的顺序如何）。因此，为了使我们的工作数据文件与该聚合文件匹配，我们需要确保从诊所1到诊所10也订购了来自各个诊所的所有数据。这可以通过使用sort cases命令轻松完成。

要访问sort cases命令，请选择select。出现的对话框如图所示。选择您要对文件进行排序的变量（在本例中为Clinic），并将其拖动到标有“排序依据”的区域（或单击）。可以选择按升序排列文件（诊所1到诊所10），或降序排列（前往诊所1的诊所10）。单击对文件排序。

下一步是在汇总文件中使用这些临床方法，以将BDI变量放在我们的主文件中。为此，我们需要使用match files命令，可以通过选择进行访问。这将打开一个对话框，其中列出了所有打开的数据文件（在我的情况下，除了我正在工作的文件之外，其他所有文件都没有打开，因此该空间为空白）或询问您选择SPSS数据文件。单击并导航到您决定存储聚合值文件的位置（在我的情况下为aggr.sav）。选择此文件，然后单击以返回到对话框。然后单击进入下一个对话框。

在下一个对话框中，我们需要匹配两个文件，这只是告诉SPSS两个文件已连接。为此，请单击。然后，我们还需要专门连接Clinic变量上的文件。为此，select告诉SPSS无效的数据集（即，汇总分数文件）应视为与关键变量上的工作数据文件匹配的值表。我们需要选择此关键变量是什么。我们要匹配Clinic变量上的文件，因此在“排除的变量”列表中选择此变量，并将其拖到标有“关键变量”的空间（或单击）。

数据编辑器现在应包含一个新变量BDI_mean，其中包含我们文件aggr.sav中的值。基本上，SPSS已匹配诊所变量的文件，因此BDI_mean中的值对应于各个诊所的平均值。因此，当临床变量为1时，BDI_mean已设置为25.19，但是当临床变量为2时，BDI_mean已设置为31.32。我们可以再次在compute命令中使用这些值来使BDI中心化。通过选择访问计算命令。在出现的对话框（图7）中，在标有“目标变量”的框中输入名称BDI_Group_Centred，然后单击并为变量指定一个更具描述性的名称。选择变量BDI并将其拖到标有“数字表达式”的区域，然后单击，然后键入“ BDI_mean”或选择此变量并将其拖到标有“目标变量”的框中。单击，将创建一个新变量，其中包含以组为中心的均值。

另外，可以使用以下语法来完成所有操作：

AGGREGATE/OUTFILE='C:\Users\Dr. Andy Field\Documents\Academic\Data\aggr.sav'/BREAK=Clinic/BDI_mean=MEAN(BDI).
SORT CASES BY Clinic(A).
MATCH FILES /FILE=*/TABLE='C:\Users\Dr. Andy Field\Documents\Academic\Data\aggr.sav'/BY Clinic.
EXECUTE.
COMPUTE BDI_Group_Centred=BDI − BDI_mean.
EXECUTE.

要访问“重组数据向导”，请选择。向导中的步骤如图所示。在第一个对话框中，您需要说是否要将变量转换为案例，还是将案例转换为变量。我们在不同的列（变量）中具有不同的时间级别，并且希望它们在不同的行（案例）中，因此我们需要选择。单击以移至下一个对话框。该对话框询问您是要从旧数据文件的不同列中在新数据文件中仅创建一个新变量，还是要创建多个新变量。

在我们的案例中，我们将创建一个代表生活满意度的变量。默认，SPSS在新数据文件中创建一个名为id的变量，该变量告诉您数据来自哪个样本（即原始数据文件的哪一行）。它通过使用原始数据文件中的案例编号来实现。然后从数据文件中选择一个变量以充当新数据文件中的标签。

其余对话框非常简单。接下来的两个处理索引变量。SPSS创建一个新变量，该变量将告诉你数据源自哪一列。在我们有四个时间点的情况下，这将意味着变量只是一个从1到4的数字序列。

多层（等级）线性模型

将BDI，年龄和性别包括在内作为固定效果预测指标。

选择，然后通过从变量列表中选择Clinic并将其拖动到标有Subjects的框中来指定变量（或单击）。

单击以移至主对话框。首先，我们必须指定结果变量，即手术后的生活质量（QoL），因此选择Post_QoL并将其拖动到标有因变量的空间（或单击）。

我们需要将预测变量作为固定效应添加到我们的模型中，因此单击，按住Ctrl并在标记为Factors和Covariates的列表中选择Base_QoL，Surgery，Age，性别，Reason和BDI。

现在，我们需要请求随机截距和随机斜率以达到手术效果。

单击并选择。单击以返回到主对话框。在主对话框中，单击并请求参数估计和协方差参数的检验。单击以返回到主对话框。要运行分析。输出如下：

就此新模型的整体拟合而言，我们可以使用对数似然统计：

卡方统计的临界值为7.81（p <.05，df = 3）；因此，这一变化意义重大。包括这三个预测变量可以改善模型的拟合度。年龄，F（1，150.83）= 37.32，p <.001，BDI，F（1，260.83）= 16.74，p <.001，显着预测了手术后的生活质量，但性别没有，F（1，264.48 ）= 0.90，p = 0.34。包括这些因素的主要区别在于，Reason的主要影响变得不显着，并且Reason×Surgery交互作用变得更加重要（其b从4.22，p = .013变为5.02，p = .001）。

我们可以通过拆分并运行更简单的分析来分解此交互，如本文所述（没有交互和Reason的主要影响，但包括Base_QoL，Surgery，BDI，Age和Gender）。如果进行这些分析，将获得输出中所示的参数表。对于那些只为改变外观而进行手术的患者，手术显着预测了手术后的生活质量，b = –3.16，t（5.25）= –2.63，p = .04。与不包括年龄，性别和BDI的情况不同，这种影响现在很明显。负系数表明，与对照组相比，这些人的手术后生活质量较低。但是，对于那些通过手术解决身体问题的人，手术并不能显着预测生活质量，b = 0.67，t（10.59）= 0.58，p = 0.57。从本质上讲，年龄，性别和BDI的纳入对后一组几乎没有什么影响。但是，该斜率是正的，表明接受手术治疗的人的生活质量得分比等候名单上的得分高（尽管不是很明显！）。因此，相互作用的影响反映了在进行身体问题手术的患者（轻微的正斜率）和仅出于虚荣心进行手术的患者（负的斜率）中手术斜率作为生活质量预测指标的差异。

改变外貌的手术

身体问题的手术

最受欢迎的见解

1.基于R语言的lmer混合线性回归模型

2.R语言用Rshiny探索lme4广义线性混合模型（GLMM）和线性混合模型（LMM）

3.R语言线性混合效应模型实战案例

4.R语言线性混合效应模型实战案例2

5.R语言线性混合效应模型实战案例

6.线性混合效应模型Linear Mixed-Effects Models的部分折叠Gibbs采样

7.R语言LME4混合效应模型研究教师的受欢迎程度

8.R语言中基于混合数据抽样(MIDAS)回归的HAR-RV模型预测GDP增长

9.使用SAS，Stata，HLM，R，SPSS和Mplus的分层线性模型HLM

拓端tecdat|SPSS中的多层（等级）线性模型Multilevel linear models研究整容手术数据相关推荐

拓端tecdat荣获掘金社区入驻新人奖
2021年7月,由掘金发起了"入驻成长礼"颁奖活动.本次活动邀请到知名开发者.服务机构代表等业界人士. 据了解,掘金社区"新入驻创作者礼"主要对已经积累了一定历 ...
拓端tecdat荣获2022年度51CTO博主之星
相信技术,传递价值,这是51CTO每一个技术创作者的动力与信念,2022 年度,拓端tecdat 作为新锐的数据分析咨询公司,在51CTO平台上,不断的输出优质的技术文章,分享前沿创新技术,输出最佳生 ...
拓端tecdat|bilibili视频流量数据潜望镜
最近我们被客户要求撰写关于bilibili视频流量的研究报告,包括一些图形和统计输出. 最新研究表明,中国有超过7亿人在观看在线视频内容.Bilibili,被称为哔哩哔哩或简称为B站,是中国大陆第二个 ...
拓端tecdat|R语言向量误差修正模型 (VECMs)分析长期利率和通胀率影响关系
最近我们被客户要求撰写关于向量误差修正模型的研究报告,包括一些图形和统计输出. 向量自回归模型估计的先决条件之一是被分析的时间序列是平稳的.但是,经济理论认为,经济变量之间在水平上存在着均衡关系,可以 ...
拓端tecdat|R语言线性回归和时间序列分析北京房价影响因素可视化案例
最近我们被客户要求撰写关于北京房价影响因素的研究报告,包括一些图形和统计输出. 目的房价有关的数据可能反映了中国近年来的变化: 人们得到更多的资源(薪水),期望有更好的房子人口众多独生子女政策: ...
拓端tecdat|R语言用LOESS(局部加权回归)季节趋势分解（STL）进行时间序列异常检测
最近我们被客户要求撰写关于LOESS(局部加权回归)的研究报告,包括一些图形和统计输出. 这篇文章描述了一种对涉及季节性和趋势成分的时间序列的中点进行建模的方法.我们将对一种叫做STL的算法进行研究, ...
拓端tecdat|R语言逻辑回归（Logistic回归）模型分类预测病人冠心病风险
最近我们被客户要求撰写关于冠心病风险的研究报告,包括一些图形和统计输出. 相关视频:R语言逻辑回归(Logistic回归)模型分类预测病人冠心病风险逻辑回归Logistic模型原理和R语言分类预测冠 ...
用SPSS估计HLM多层（层次）线性模型模型
原文链接:http://tecdat.cn/?p=3230 作为第一步,从一个不包含协变量的空模型开始(点击文末"阅读原文"获取完整代码数据). 每所学校的截距,β 0J,然后设置 ...
如何给图片赋值_医学数据的变量类型及在SPSS中的赋值方法（医学统计前的重要步骤）——【杏花开医学统计】...
杏花开生物医药统计一号在手,统计无忧! 关注医学数据的变量类型及在SPSS中的赋值方法关键词:医学数据,变量类型,赋值,spss 下方为视频版和音频版,含软件操作步骤导读不同的医学数 ...
社交平台中的会员等级制度建立需要注意哪些问题
只是讲一点自己的理解,不一定对,希望大神们来喷,社交软件中的用户等级需要区别于其他类型下的用户等级,社交软件内的用户等级在我的理解当中,很大的作用是为了活跃用户和筛选优质和资深用户.同时,能够在平台内 ...

拓端tecdat|SPSS中的多层（等级）线性模型Multilevel linear models研究整容手术数据