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