如何在semPlots中为特定条目（如Q3、Q9）设置不同颜色

希望在semPlots生成的结构方程模型图中，为Q3和Q9等特定条目设置不同颜色以突出显示。尝试数周学习与反复试验仍未找到实现方法，以下是数据及代码：

library(lavaan)
library(semPlot)

baseData<- c(1,2,3,4,5) # categories
s <- 500 # sample size
set.seed(11)
catdata <- data.frame(
  Q1 = sample(baseData, size = s, replace = TRUE),
  Q2 = sample(baseData, size = s, replace = TRUE),
  Q3 = sample(baseData, size = s, replace = TRUE),
  Q4 = sample(baseData, size = s, replace = TRUE),
  Q5 = sample(baseData, size = s, replace = TRUE),
  Q6 = sample(baseData, size = s, replace = TRUE),
  Q7 = sample(baseData, size = s, replace = TRUE),
  Q8 = sample(baseData, size = s, replace = TRUE),
  Q9 = sample(baseData, size = s, replace = TRUE),
  Q10 = sample(baseData, size = s, replace = TRUE),
  Q11 = sample(baseData, size = s, replace = TRUE),
  Q12 = sample(baseData, size = s, replace = TRUE),
  Q13 = sample(baseData, size = s, replace = TRUE),
  Q14 = sample(baseData, size = s, replace = TRUE),
  Q15 = sample(baseData, size = s, replace = TRUE)
)

my_model <- '
Factor 1 =~ Q1 + Q2 + Q3 + Q4 + Q5
Factor 2 =~ Q6 + Q7 + Q8 + Q9 + Q10
Factor 3 =~ Q11 + Q12 + Q13 + Q14 + Q15
'

model_fit <- cfa(my_model, data = catdata)

plot <- semPaths(model_fit, whatLabels="std", layout="circle2", style = "lisrel", 
                        curve = 2, rotation = 3, shapeLat = "ellipse", shapeMan = "circle", sizeMan = 6,
                        nCharNodes = 0, intercepts = FALSE, residuals = FALSE,  exoCov = TRUE, edge.label.cex = 0.75, nDigits = 2, edge.label.color = "black")

可以通过semPaths()的colorMan参数为特定观测变量设置自定义颜色，步骤如下：

1. 确认观测变量顺序
   semPlots会按照模型中变量出现的顺序排列观测变量，你的模型中顺序为：Q1、Q2、Q3、Q4、Q5、Q6、Q7、Q8、Q9、Q10、Q11、Q12、Q13、Q14、Q15。若不确定顺序，可运行lavNames(model_fit, "ov")查看。

2. 构建颜色向量
   创建一个与观测变量数量匹配的颜色向量，默认颜色设为常规颜色（如黑色），再将目标变量对应的位置替换为突出颜色（如红色）。

3. 修改绘图代码
   在semPaths()中添加colorMan参数，传入构建好的颜色向量。

修改后的完整代码：

library(lavaan)
library(semPlot)

baseData<- c(1,2,3,4,5) # categories
s <- 500 # sample size
set.seed(11)
catdata <- data.frame(
  Q1 = sample(baseData, size = s, replace = TRUE),
  Q2 = sample(baseData, size = s, replace = TRUE),
  Q3 = sample(baseData, size = s, replace = TRUE),
  Q4 = sample(baseData, size = s, replace = TRUE),
  Q5 = sample(baseData, size = s, replace = TRUE),
  Q6 = sample(baseData, size = s, replace = TRUE),
  Q7 = sample(baseData, size = s, replace = TRUE),
  Q8 = sample(baseData, size = s, replace = TRUE),
  Q9 = sample(baseData, size = s, replace = TRUE),
  Q10 = sample(baseData, size = s, replace = TRUE),
  Q11 = sample(baseData, size = s, replace = TRUE),
  Q12 = sample(baseData, size = s, replace = TRUE),
  Q13 = sample(baseData, size = s, replace = TRUE),
  Q14 = sample(baseData, size = s, replace = TRUE),
  Q15 = sample(baseData, size = s, replace = TRUE)
)

my_model <- '
Factor 1 =~ Q1 + Q2 + Q3 + Q4 + Q5
Factor 2 =~ Q6 + Q7 + Q8 + Q9 + Q10
Factor 3 =~ Q11 + Q12 + Q13 + Q14 + Q15
'

model_fit <- cfa(my_model, data = catdata)

# 定义观测变量颜色：默认黑色，Q3和Q9设为红色
man_colors <- rep("black", 15)
man_colors[c(3, 9)] <- "red"

plot <- semPaths(model_fit, whatLabels="std", layout="circle2", style = "lisrel", 
                 curve = 2, rotation = 3, shapeLat = "ellipse", shapeMan = "circle", sizeMan = 6,
                 nCharNodes = 0, intercepts = FALSE, residuals = FALSE,  exoCov = TRUE, edge.label.cex = 0.75, nDigits = 2, edge.label.color = "black",
                 colorMan = man_colors) # 添加自定义颜色参数

额外说明 #

• 可替换"red"为任意支持的颜色值（如"blue"或十六进制颜色码"#FF6B6B"）。
• 若需要修改节点边框颜色，可使用border.colorMan参数，用法与colorMan一致。

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