data(mexchn)
fo <- list(self = xsayself ~ 1,
vign = xsay3 ~ 1,
tau = ~1 )
out1 <- chopit( fo, mexchn, options= anchors.options( vign.var = "hetero", normalize="hilo" ))
summary(out1)
Run the code above in your browser using DataLab