适合 Lavaan CFA 的模型。这段代码有什么问题?
Model fit for Lavaan CFA. What is the problem with this code?
我已经使用 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)
我快速查看了数据,似乎 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")
我已经使用 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)
我快速查看了数据,似乎 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")