R语言update()函数更新measEq.syntax生成的模型公式失败

错误原因 #

1. 初始调用measEq.syntax时未开启均值结构（meanstructure=TRUE）：模型默认不包含均值项，后续添加~1形式的均值约束会导致函数解析多组语法时出现参数长度判断异常，触发'length = 3' in coercion to 'logical(1)'错误。
2. 误用change.syntax参数：该参数用于替换原有模型语法，而非添加新约束，会破坏measEq.syntax生成的多组模型结构，进一步引发解析错误。

正确解决方案 #

方案一：开启均值结构后用add.syntax追加约束 #

先在生成初始模型时开启均值结构，再通过add.syntax参数添加均值约束，避免破坏原有模型结构：

HS.model <- ' visual  =~ x1 + x2 + x3
              textual =~ x4 + x5 + x6
              speed   =~ x7 + x8 + x9
              general =~ visual + textual + speed'

library(semTools)
# 生成带均值结构的多组模型语法
syntax.config <- measEq.syntax(configural.model = HS.model,
                               data = HolzingerSwineford1939,
                               ID.fac = "ul",
                               group = "school",
                               meanstructure = TRUE) # 必须开启均值结构

# 定义均值约束语句
fixMeans.c <- '
textual ~ c(0, 0)*1
speed ~ c(0, 0)*1
general ~ c(0, 0)*1
'
# 用add.syntax追加约束，而非替换
syntax.config <- update(syntax.config, add.syntax = fixMeans.c)

方案二：直接在初始模型中包含均值约束 #

如果不需要分步构建模型，可以直接在初始模型语法中写入均值约束，再生成多组模型：

HS.model_with_means <- ' visual  =~ x1 + x2 + x3
                        textual =~ x4 + x5 + x6
                        speed   =~ x7 + x8 + x9
                        general =~ visual + textual + speed
                        # 嵌入均值约束
                        textual ~ c(0,0)*1
                        speed ~ c(0,0)*1
                        general ~ c(0,0)*1'

library(semTools)
syntax.config <- measEq.syntax(configural.model = HS.model_with_means,
                               data = HolzingerSwineford1939,
                               ID.fac = "ul",
                               group = "school",
                               meanstructure = TRUE)

内容的提问来源于stack exchange，提问作者Sinval