ksi =~ NA*x1 + x2 + x3
eta =~ y1 + y2 + y3
eta ~ label(c('gamma', 'gamma'))*ksi
ksi ~~ c(1, 1)*ksi
