I have been running an unconstrained gllvm with a row effect, and of course, gllvm returns LV scores per row (= per sample). What i would like are LV scores on a per site (e.g. per Lake) basis. Is that possible? The ugliest way i could think of to achieve this is to take the mean of each LV per site.
doug
gllvm( y = y, family = binomial() , num.lv = lvchosen, studyDesign = Site, row.eff = ~ (1 | Lake), control.start = list(n.init = 10) )
hello!
I have been running an unconstrained gllvm with a row effect, and of course, gllvm returns LV scores per row (= per sample). What i would like are LV scores on a per site (e.g. per Lake) basis. Is that possible? The ugliest way i could think of to achieve this is to take the mean of each LV per site.
doug
gllvm( y = y, family = binomial() , num.lv = lvchosen, studyDesign = Site, row.eff = ~ (1 | Lake), control.start = list(n.init = 10) )