Опосредованная модерация с фиктивными переменными lavaan

Question

Я пытаюсь настроить опосредованную модель умеренности, используя lavaan, но не знаю, как ее кодировать Я хочу воспроизвести модель процесса Хейса № 7 , но с двумя независимыми переменными и взаимодействием между ними. Независимые переменные являются фиктивными, поэтому я не думаю, что смогу просто умножить их вместе, чтобы получить термин взаимодействия. Но переменные модератор, посредник и результат являются непрерывными.

Вот некоторые примеры данных с моим кодом:

df <- data.frame( iv1 = c(rep(0, 50), rep(1, 50)),
                  iv2 = c(rep(c(0,1), 50)),
                  moderator = rnorm(mean = 5.19, n = 100, sd = 1.22),
                  mediator = rnorm(mean = 4.02, n = 100, sd = 1.90)
                  dv = rnorm(mean = 3.27, n = 100, sd = 1.35)
                  interaction1 = iv1*iv2
                  interaction2 = iv1*moderator
                  interaction3 = iv2*moderator
                  interaction4 = iv1*iv2*moderator)

library(lavaan)
mod <- '
        mediator ~ a1*moderator + a2*iv1 + a3*iv2 + a4*interaction1 +                   
a5*interaction2 + a6*interaction3 + a7*interaction4
        dv ~ c1*iv1 + c2*iv2 + c3*interaction1 + b*mediator

        direct:= c1 + c2 + c3
        indirect:= (a1 + a2 + a3) * (a4 + a5 + a6 + a7) * b
        total:= direct + indirect'
fit <- sem(mod, data = df, std.lv=T, se = "boot", bootstrap = 500)
summary(fit)

Спасибо!