Code in R is provided for multilevel slopes-as-outcomes analysis. The current implementation minimizes either the full negative log-likelihood or the restricted negative log-likelihood using the scoring method. We use a simple repeated measures example to illustrate the analysis.