gistfile1.r

data(iris)

model.iris <- cfa()
F1: Sepal.Length, Sepal.Width, Petal.Length, Petal.Width

# 特殊記号はエスケープする必要がある
model.iris2 <- update(model.iris)
replace, lam\[Petal.Width:F1\], b2

model.iris3 <- update(model.iris)
replace, lam\[Petal.Width:F1\], b3

model.mg.iris <- multigroupModel(
    model.iris, model.iris2, model.iris3,
    groups=c("setosa", "versicolor", "virginica"))

sem.mg.iris <- sem(model.mg.iris, data=iris, group="Species",
    formula= ~ Sepal.Length + Sepal.Width + Petal.Length + Petal.Width)
summary(sem.mg.iris)