lambda_mid_csv<-read.csv("lambda_mid_rep.csv") library(rptR) library(lme4) rptGaussian(formula = lambda_mid~(1|Bird.ID),data = lambda_mid_csv,grname = "Bird.ID")