I suggest to add a shortcut multiple group SEM constructor that looks like this:
model = MultipleGroupsSem(
specification = partable
data = dat,
by = column
)
which creates an SEMEnsemble with groups corresponding to the number of unique entries in column column in dataset dat.