varying = list(as.character(unique(dat[,1]))))[-1] This line should be : varying = list(as.character(unique(dat[,1][-1]))
varying = list(as.character(unique(dat[,1]))))[-1]
This line should be :
varying = list(as.character(unique(dat[,1][-1]))