如何解决适合 Lavaan CFA 的模型这段代码有什么问题?
我使用 R 中的 lavaan 包实现了验证性因素分析。我编写的代码产生了以下错误:
lavaan 0.6-7 在 2044 次迭代后没有正常结束
有人可以帮我吗?
##PACKAGES
library(lavaan)
library(semTools)
library(semPlot)
library(readxl)
library(dplyr)
##DATASET
dataset_ita <- read_excel("dataset_ita.xlsx",col_types = "numeric")
dataset_eng <- read_excel("dataset_eng.xlsx",col_types = "numeric")
dataset <- rbind(dataset_ita,dataset_eng)
##CONFIRMATORY FACTOR ANALYSIS
model <- '
WinePurchaseBehaviour =~ X5+X6+X7+X8+X9+X10+X11
WineConsumptionBehaviour =~ X12+X13+X14+X15+X16+X17+X18+X19+X20+X21+X22+X23+X24+X25
WineClubInterest =~ X26+X27+X28+X29+X30+X31+X64+X65+X66+X67
WineInvolvement =~ WinePurchaseBehaviour+WineConsumptionBehaviour+WineClubInterest
GeneralInvolvement =~ X32+X33+X34+X35+X36+X37+X38+X39+X40
FeatureInvolvement =~ X44+X46+X47+X48+X49
RitualInvolvement =~ X41+X42+X43+X45+X50
AppExperience =~ X51+X52+X53+X54+X55+X56+X57+X58
SensoryExperience =~ X59+X60+X61+X62+X63
ProductInvolvement =~ GeneralInvolvement+FeatureInvolvement+AppExperience+SensoryExperience
Purchase =~ X68+X69+X70+X71+X72+X73
Purchase =~ WineInvolvement+ProductInvolvement'
analysis <- cfa(model,data = dataset,se = "robust.sem")
summary(analysis,fit.measures=TRUE)
The datasets can be found here
解决方法
我快速查看了数据,似乎 RitualInvolvement
因子不起作用 - 它甚至不能单独起作用。 FeatureInolvement
因子使用 X47
作为系数设置为 1 效果更好,因此在列表中首先使用它。下面的模型,没有估计 RitualInvolvement
因子估计并为我收敛。
model <- '
WinePurchaseBehaviour =~ X5+X6+X7+X8+X9+X10+X11
WineConsumptionBehaviour =~ X12+X13+X14+X15+X16+X17+X18+X19+X20+X21+X22+X23+X24+X25
WineClubInterest =~ X26+X27+X28+X29+X30+X31+X64+X65+X66+X67
WineInvolvement =~ WinePurchaseBehaviour+WineConsumptionBehaviour+WineClubInterest
GeneralInvolvement =~ X32+X33+X34+X35+X36+X37+X38+X39+X40
FeatureInvolvement =~ X47+X46+X44+X48+X49
# RitualInvolvement =~ X41+X42+X43+X45+X50
AppExperience =~ X51+X52+X53+X54+X55+X56+X57+X58
SensoryExperience =~ X59+X60+X61+X62+X63
ProductInvolvement =~ GeneralInvolvement+FeatureInvolvement+AppExperience+SensoryExperience
Purchase =~ X68+X69+X70+X71+X72+X73
Purchase =~ WineInvolvement+ProductInvolvement'
analysis <- cfa(model,data = dataset,se = "robust.sem")
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。