使用semTools中parcelAllocation函数创建随机项目包裹的问题

问题背景 #

要构建包含MAAS、PFS、SFS三个潜变量的结构方程模型，需为各潜变量创建随机项目包裹，数据集为Mydata：

• MAAS含15个项目（Mydata$MAAS_1至Mydata$MAAS_15），需生成3个各含5个项目的包裹；
• PFS含7个项目（Mydata$PFS_1至Mydata$PFS_7），需生成3个各含2-3个项目的包裹；
• SFS含7个项目（Mydata$SFS_1至Mydata$SFS_7），需生成3个各含2-3个项目的包裹。

使用semTools包的parcelAllocation函数实现时，遇到两个问题：

1. 无法查看项目到包裹的具体分配结果；
2. 不知道如何生成项目平均分类型的包裹（默认是求和包裹）。

解决方案 #

1. 查看包裹分配结果 #

parcelAllocation返回的对象中，parcelList元素存储了所有随机分配的包裹列表（每个子元素对应一次分配结果），可通过以下方式查看：

查看单次分配结果（以第1次为例）

# 查看MAAS第1次包裹分配详情
print(maas.parcels$parcelList[[1]])

# 查看PFS第1次包裹分配详情
print(pfs.parcels$parcelList[[1]])

# 查看SFS第1次包裹分配详情
print(sfs.parcels$parcelList[[1]])

遍历查看所有分配结果

如果需要对比多次分配的差异，可循环遍历parcelList：

# 遍历MAAS的100次分配结果
for (i in seq_along(maas.parcels$parcelList)) {
  cat(sprintf("第 %d 次分配结果:\n", i))
  print(maas.parcels$parcelList[[i]])
  cat("\n")
}

2. 生成项目平均分包裹 #

parcelAllocation默认生成求和包裹，要得到平均分包裹，需基于分配结果手动计算均值并添加到原数据集。以下是具体实现：

单潜变量处理示例（MAAS）

假设选择第1次分配结果生成平均分包裹：

# 获取MAAS第1次分配的项目列表
maas_parcel_assign <- maas.parcels$parcelList[[1]]

# 计算每个包裹的平均分并添加到Mydata（忽略缺失值）
Mydata$maas.par1 <- rowMeans(Mydata[, maas_parcel_assign$maas.par1], na.rm = TRUE)
Mydata$maas.par2 <- rowMeans(Mydata[, maas_parcel_assign$maas.par2], na.rm = TRUE)
Mydata$maas.par3 <- rowMeans(Mydata[, maas_parcel_assign$maas.par3], na.rm = TRUE)

批量处理所有潜变量

可封装函数简化操作，同时处理缺失值：

# 定义生成平均分包裹的函数
add_mean_parcels <- function(data, parcel_assign, parcel_names) {
  for (parcel in parcel_names) {
    # 计算均值，忽略缺失值
    data[[parcel]] <- rowMeans(data[, parcel_assign[[parcel]]], na.rm = TRUE)
  }
  return(data)
}

# 处理MAAS
Mydata <- add_mean_parcels(Mydata, maas.parcels$parcelList[[1]], maas.parcel.names)

# 处理PFS
Mydata <- add_mean_parcels(Mydata, pfs.parcels$parcelList[[1]], pfs.parcel.names)

# 处理SFS
Mydata <- add_mean_parcels(Mydata, sfs.parcels$parcelList[[1]], sfs.parcel.names)

3. 补充：固定包裹大小（可选） #

如果需要严格控制每个包裹的项目数量（比如MAAS必须5个/包裹，PFS为3、2、2），可在调用parcelAllocation时通过parcel.size参数指定：

# MAAS指定每个包裹5个项目
maas.parcels <- parcelAllocation(maas.mod.parcels, data = Mydata, 
                  nAlloc = 100, parcel.names = maas.parcel.names, 
                  item.syntax = item.syntaxMaas, std.lv = TRUE,
                  parcel.size = c(5, 5, 5)) # 明确各包裹项目数

# PFS指定包裹大小为3、2、2
pfs.parcels <- parcelAllocation(pfs.mod.parcels, data = Mydata, 
                                nAlloc = 100, parcel.names = pfs.parcel.names, 
                                item.syntax = item.syntaxPfs, 
                                std.lv = TRUE, parallel = "snow",
                                parcel.size = c(3, 2, 2))

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