Code is here:
brown<-melt(brown, id=c("date","variety", "GDD")) brown$variety <- as.factor(brown$variety) #set variety as