【问题标题】:Model fit for Lavaan CFA. What is the problem with this code?模型适合 Lavaan CFA。这段代码有什么问题?
【发布时间】:2020-12-23 11:11:28
【问题描述】:

我已经使用 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

【问题讨论】:

  • 欢迎来到 SO。如果我们有重现问题的代码和数据,我们会更容易为您提供帮助。有一个有用的post 可以制作一个可能有用的可重现的 R 示例。
  • 谢谢!我通过添加可重现的 R 示例和 github 存储库的链接来编辑问题,以便下载验证性因子分析所需的数据集。

标签: r r-lavaan


【解决方案1】:

我快速查看了数据,似乎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")

【讨论】:

  • 非常感谢!同时我发现,通过从RitualInvolvement中去掉变量X43,模型收敛了……这听起来很奇怪!!我一直在努力!!
猜你喜欢
  • 2015-04-25
  • 2018-09-09
相关资源
最近更新 更多