手把手教你用lavaan做路径分析，零基础也能快速上手

关注

原创最新推荐文章于 2026-09-23 11:12:05 发布 1k 阅读 AI 写同款· GEO 优化 ›

基于R语言的结构方程：lavaan简明教程 [中文翻译版] lavaan简明教程 [中文翻译版] 译者注：此文档原作者为比利时Ghent大学的Yves Rosseel博士，lavaan亦为其开发，完全开源、免费。我在学习的时候顺手翻译了一下，向Yves的开源精神致敬。此翻译因偷懒部分删减，但也有增加，有错误请留言「转载请注明出处」目录 lavaan简明教程 [中文翻译版] 目录摘要在开始之前安装lavaan包模型语法例1：验证性... 阅读详情

第一章：路径分析与结构方程模型入门

路径分析与结构方程模型（Structural Equation Modeling, SEM）是社会科学、心理学及行为科学研究中广泛使用的多变量统计方法。它结合了因子分析和路径分析的优势，能够同时处理测量误差与潜在变量之间的复杂关系，适用于验证理论模型的拟合程度。

核心概念解析

潜在变量：无法直接观测但可通过多个指标变量推断的抽象构念，如“满意度”或“认知能力”
观测变量：实际测量得到的数据，作为潜在变量的代理指标
路径系数：表示变量间直接影响强度的估计值，通常通过标准化回归权重解释
模型识别：确保参数有唯一解的数学条件，常见规则是自由度非负且参数数量不超过样本矩

基础建模步骤

明确研究假设并绘制理论模型图
收集数据并检验变量的信效度
使用软件指定模型结构并进行参数估计
评估模型拟合指数，如 CFI > 0.95、RMSEA < 0.06
根据修正指数调整模型（若理论允许）

R语言实现示例


# 加载lavaan包用于SEM分析
library(lavaan)

# 定义潜变量与路径关系
model <- '
  # 测量模型
  visual  =~ x1 + x2 + x3
  textual =~ x4 + x5 + x6
  speed   =~ x7 + x8 + x9

  # 结构模型
  textual ~ visual
  speed   ~ visual
'

# 拟合模型
fit <- sem(model, data = HolzingerSwineford1939)

# 输出拟合摘要
summary(fit, standardized = TRUE, fit.measures = TRUE)

graph LR A[观测变量x1] --> B(潜在变量: 视觉能力) A2[观测变量x2] --> B A3[观测变量x3] --> B B --> C[文本能力] B --> D[速度能力] C --> E[结果输出] D --> E

拟合指数	推荐阈值	含义说明
CFI	> 0.95	比较拟合指数，越接近1越好
RMSEA	< 0.06	近似误差均方根，反映模型简约性
SRMR	< 0.08	标准化残差均值，衡量整体偏差

第二章：lavaan基础与数据准备

2.1 理解路径分析与结构方程模型的基本概念

路径分析是结构方程模型（SEM）的基础，用于检验变量间的直接与间接因果关系。它通过图形化路径图表示假设关系，并利用协方差矩阵进行统计验证。

模型构成要素

结构方程模型包含测量模型与结构模型两部分。前者连接潜变量与观测变量，后者描述潜变量之间的因果路径。

潜变量（Latent Variables）：无法直接测量的构念，如“用户满意度”
观测变量（Observed Variables）：可通过数据采集的实际指标
路径系数：表示变量间影响强度的估计参数

典型模型表达式

model <- '
  # 测量模型
  Satisfaction =~ sat1 + sat2 + sat3
  Quality =~ qual1 + qual2 + qual3

  # 结构模型
  Satisfaction ~ Quality + Service
'

该代码定义了一个包含测量与结构关系的SEM模型，Satisfaction 作为潜变量由三个观测项构成，并受 Quality 和 Service 影响。路径方向用 ~ 表示，左侧为结果变量，右侧为预测变量。

2.2 R语言中lavaan包的安装与初步使用

安装lavaan包

在R环境中，可通过CRAN直接安装lavaan包。执行以下命令完成安装：

install.packages("lavaan")

该命令会自动下载并配置lavaan及其依赖包（如MASS、stats等），确保结构方程建模所需的核心功能完整可用。

加载与简单模型定义

安装完成后，需加载包以启用相关函数：

library(lavaan)

随后可定义一个简单的验证性因子分析模型。例如：

model <- '
  # 潜变量定义
  visual =~ x1 + x2 + x3
  textual =~ x4 + x5 + x6
'

上述代码中，=~表示潜变量对观测变量的加载关系，注释通过#添加，提升模型脚本可读性。此结构为后续数据拟合提供语法基础。

2.3 数据读入、清洗与描述性统计分析

数据读入与初步观察

使用Pandas进行数据读入是分析的第一步。常见格式如CSV、Excel可直接加载：

import pandas as pd
df = pd.read_csv('data.csv')
print(df.head())

该代码读取CSV文件并显示前5行，便于快速了解数据结构。

数据清洗关键步骤

清洗包括处理缺失值、去重和类型转换。例如：

df.drop_duplicates(inplace=True)
df['age'].fillna(df['age'].median(), inplace=True)
df['date'] = pd.to_datetime(df['date'])

上述操作确保数据一致性与完整性，为后续分析奠定基础。

描述性统计分析

通过内置方法快速获取统计量：

统计指标	含义
count	非空值数量
mean	平均值
std	标准差

调用 df.describe() 可输出数值型字段的汇总统计，辅助识别异常分布。

2.4 变量类型识别与正态性检验

变量类型的初步识别

在数据分析前期，准确识别变量类型是建模基础。常见变量分为数值型（连续、离散）、分类型和有序型。可通过 pandas 快速查看：

import pandas as pd
data = pd.read_csv('dataset.csv')
print(data.dtypes)

该代码输出各列数据类型，float64 和 int64 通常为数值型，object 多为分类变量，需进一步判断。

正态性检验方法

对连续变量需检验是否服从正态分布，常用 Shapiro-Wilk 检验：

from scipy.stats import shapiro
stat, p = shapiro(data['age'])
print(f'Statistic: {stat:.3f}, p-value: {p:.3f}')

当 p 值小于 0.05 时，拒绝原假设，数据非正态。此时应考虑变换（如 log）或使用非参数方法。

2.5 构建理论模型前的变量关系探索

在构建理论模型之前，深入分析变量间的潜在关系是确保模型科学性的关键步骤。通过探索性数据分析（EDA），可以识别出协变量之间的相关性、共线性及非线性交互效应。

变量相关性热力图分析

使用皮尔逊相关系数评估连续变量间的线性关联，并通过热力图可视化结果：


import seaborn as sns
import matplotlib.pyplot as plt

corr_matrix = df.corr(numeric_only=True)
sns.heatmap(corr_matrix, annot=True, cmap='coolwarm', center=0)
plt.title("Variable Correlation Heatmap")
plt.show()

该代码段生成变量间相关性的热力图，其中 `annot=True` 显示具体数值，`cmap='coolwarm'` 提供直观的颜色映射，便于快速识别高相关性变量对，为后续特征选择提供依据。

变量筛选建议

剔除与其他变量相关性高于0.9的特征，以降低多重共线性风险；
保留与目标变量相关性显著的输入变量；
考虑引入交互项或多项式特征以捕捉非线性模式。

第三章：构建与设定路径模型

3.1 使用lavaan语法定义潜变量与观测变量

在结构方程模型中，lavaan包通过简洁的语法实现潜变量与观测变量之间的关系定义。用户使用“=~”操作符表示潜变量由哪些观测变量构成。

基本语法结构


model <- '
  # 潜变量定义
  visual  =~ x1 + x2 + x3
  textual =~ x4 + x5 + x6
  speed   =~ x7 + x8 + x9
'

上述代码中，visual 是一个潜变量，由观测变量 x1 到 x3 构成。“=~”左侧为潜变量，右侧为对应的观测指标。默认情况下，lavaan将第一个载荷固定为1以识别模型。

参数说明

“=~”：表示测量关系，即潜变量生成观测变量
“+”：连接多个观测变量
注释：使用“#”添加说明，提升模型可读性

3.2 设定直接效应、间接效应与协方差路径

在结构方程模型中，明确变量间的因果关系是建模的核心。设定路径时需区分直接效应与间接效应，并合理引入协方差路径以反映潜在相关性。

直接效应与间接效应的定义

直接效应指一个变量对另一变量的直接影响，而间接效应则通过中介变量传递。例如，X 影响 M，M 再影响 Y，则 X→Y 为直接效应，X→M→Y 构成间接效应。

路径设定示例


# 使用lavaan语法设定路径
model <- '
  # 直接效应
  Y ~ c*X
  # 间接路径
  M ~ a*X
  Y ~ b*M
  # 间接效应
  ab := a*b
  # 协方差路径
  X ~~ M
'

上述代码中，c 表示 X 对 Y 的直接效应，a*b 计算间接效应，X ~~ M 允许自变量与中介变量存在协方差，提升模型拟合合理性。

3.3 模型识别与参数约束原则

在构建可解释性强、泛化能力高的机器学习模型时，模型识别与参数约束是关键环节。合理的参数约束不仅能提升训练稳定性，还能有效防止过拟合。

参数正则化策略

常用的约束手段包括L1和L2正则化，通过在损失函数中引入惩罚项限制参数规模：

loss = mse_loss + lambda_l1 * torch.norm(param, 1) + lambda_l2 * torch.norm(param, 2)

其中，L1范数促进稀疏性，适用于特征选择；L2范数抑制极端权重，增强鲁棒性。

约束类型对比

约束类型	适用场景	优势
L1正则化	高维稀疏数据	自动特征筛选
L2正则化	参数敏感模型	提升数值稳定性

硬约束与软约束

硬约束：直接限定参数取值范围，如梯度裁剪
软约束：通过损失函数间接引导参数分布

实际应用中常结合两者以兼顾收敛性与模型性能。

第四章：模型拟合评估与结果解读

4.1 执行模型拟合并查看输出结果

在完成数据预处理与特征工程后，进入模型训练的核心阶段。首先调用 `fit()` 方法对训练数据执行拟合操作：


model.fit(X_train, y_train, epochs=50, batch_size=32, validation_split=0.2)

该代码启动神经网络的训练流程，其中 `epochs=50` 表示完整遍历训练集50次，`batch_size=32` 控制每批处理32个样本，`validation_split=0.2` 自动保留20%数据用于验证性能，避免过拟合。

输出日志解析

训练过程中将逐轮输出损失值与评估指标，典型日志如下：

Epoch 1/50 - loss: 0.6821 - val_loss: 0.5234
Epoch 2/50 - loss: 0.4910 - val_loss: 0.4567
...

观察损失是否稳定下降，是判断模型收敛的关键依据。

4.2 关键拟合指标解读：CFI, TLI, RMSEA, SRMR

在结构方程模型中，模型拟合优度通过多个统计指标综合评估。常用的指标包括CFI、TLI、RMSEA和SRMR，它们从不同角度反映模型与数据的匹配程度。

拟合指标及其解释

CFI（Comparative Fit Index）：比较独立模型与目标模型的卡方值，通常大于0.95表示良好拟合；
TLI（Tucker-Lewis Index）：考虑自由度调整，对复杂模型更敏感，建议值 ≥ 0.90；
RMSEA（Root Mean Square Error of Approximation）：反映近似误差，小于0.06为理想；
SRMR（Standardized Root Mean Residual）：残差均值标准化，越小越好，一般 < 0.08 可接受。

典型输出示例


CFI = 0.978
TLI = 0.965
RMSEA = 0.042 [90% CI: 0.030, 0.053]
SRMR = 0.038

该结果表明模型具备良好的整体拟合性，所有指标均处于可接受范围内。

4.3 路径系数显著性检验与标准化输出

在结构方程模型中，路径系数的显著性检验是验证变量间关系强度的关键步骤。通常采用Z检验或t检验评估估计参数是否显著不为零。

显著性检验流程

计算路径系数的标准误
构造Z统计量：系数除以其标准误
对比显著性水平（如p < 0.05）判断显著性

标准化结果输出示例


from semopy import Model
model = Model()
model.load_syntax(syntax)
model.fit(data)
print(model.summary())

上述代码执行模型拟合并输出包含路径系数及其显著性指标的汇总表。summary()方法自动提供非标准化与标准化系数、标准误、Z值和p值。

结果解读参考表

路径	系数	p值	显著性
X → Y	0.47	0.003	**
M → Y	0.32	0.061	ns

4.4 修正模型：基于MI指数的路径增删策略

在动态网络建模中，路径结构的合理性直接影响模型表达能力。为优化拓扑结构，引入互信息（MI）指数作为路径重要性度量指标，通过量化节点间信息传递强度，识别冗余或缺失连接。

MI指数驱动的路径裁剪

当两节点间MI值低于阈值 θ_min，视为信息耦合弱，触发删除机制：


if mutual_info[u][v] < theta_min:
    graph.remove_edge(u, v)

该策略有效抑制过拟合，提升模型泛化性能。

低密度路径补充机制

针对局部信息流通阻塞区域，检测MI梯度突变点，动态插入新边：

计算邻域MI方差 σ²_MI
若 σ²_MI > τ，激活路径增补模块
引入带权重的反向传播校正ΔW

参数	含义	默认值
θ_min	最小MI阈值	0.15
τ	方差触发阈值	0.22

第五章：总结与进阶学习建议

持续构建项目以巩固技能

实际项目是检验技术掌握程度的最佳方式。建议从微服务架构入手，例如使用 Go 语言构建一个具备 JWT 鉴权、REST API 和 PostgreSQL 数据库连接的用户管理系统。

// 示例：Go 中的简单 HTTP 路由处理
func main() {
    r := mux.NewRouter()
    r.HandleFunc("/api/users", getUsers).Methods("GET")
    r.HandleFunc("/api/users", createUser).Methods("POST")
    http.ListenAndServe(":8080", r)
}
// 添加中间件如日志记录和身份验证可提升系统健壮性