我正在使用lavaan并且只观察变量(没有潜在变量)。 我想在模型中包含一个交互术语,但不知道如何做到这一点。
这就是我所拥有的
model4 <-'
interac =~ var1 * var2
Ent ~ age
presu ~ age + interac
protein ~ age + fat
fat ~ age
tempo ~ age +interac+protein
score ~sex+education+presu+tempo
'
fit <- sem(model4, data=mydata)
summary(fit4, fit.measures=TRUE)
(所有变量在开始之前都已缩放,因为我遇到了一些问题,其中一些变量比其他变量大100倍)。
我想知道这是否正确?我不会在回归中产生互动的主要影响吗?不应该包括这些吗? 当我在回归中直接添加交互项(var1 * var2)时,我得到1作为估计,所以这一定是错误的......
答案 0 :(得分:0)
不,这是不正确的。对于清单变量交互,您有两种选择:
1 - 在 lavaan 之外创建交互项,例如:
mydata$interac <- mydata$var1 * mydata$var2
或
2 - 使用 :
运算符:
model4 <-'
Ent ~ age
presu ~ age + var1:var2 #interaction and age as predictors
protein ~ age + fat
fat ~ age
tempo ~ age + var1:var2 + protein #interaction, age and protein as predictors
score ~sex+education+presu+tempo
'
fit <- sem(model4, data=mydata)
summary(fit4, fit.measures=TRUE)