Sabre was removed from CRAN because certain bugs haven't been fixed by the authors and it doesn't seem to be maintained any longer. Perhaps you could switch to another package, fsay, lme4? What kind of models do you want to fit?
The basis for mixed model analysis in lme4 is a model formula of the type Y~ x+(1|z) with the function glmer you can specify the family (gaussian, binomial, poisson) so lme4 can fit bivariate models.