I am running the following multilevel regression using lme function in package nlme:
addRandomSlope<-lme(IC ~ AU+CU+IU, data =combined_database, random=~AU