I am using lavaan and have only observed variables (no latent variables). I would like to include an interaction term in the model, but not sure how to do this.
This is what I have
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)
(all variables have been scaled before starting, because I had some issues with some variables being 100 times larger than others).
I am wondering whether this is correct? I don't have the main effects of the interaction in the regression? Shouldn't these be included? When I add the interaction term directly in the regression (var1*var2), I get 1 as estimates, so that must be wrong...