df1 <- data.frame(
  CompanyID=c("Drinkers","Drinkers","Drinkers","Drinkers","Drinkers","Drinkers","Drinkers","Drinkers"
              ,"Drinkers","Drinkers", "Liquders","Liquders","Liquders","PelletCoffeeCo","PelletCoffeeCo"),
  Email= c("john@coffee.com", "john@coffee.com","john@coffee.com","john@coffee.com", "john@coffee.com",
           "john@coffee.com", "john@coffee.com", "john@coffee.com", "john@coffee.com", "john@coffee.com",
           "george@liquid.com","george@liquid.com","george@liquid.com","stacy@pelletcoffee.com",
           "stacy@pelletcoffee.com"),
  Day= c("1","2","3","4","5","6","7","8","9","10","1","2","3","1","2"),
  var1= c(4,5,5,5,2,3,2,7,6,5,7,6,6,2,3),
  var2= c(4,5,5,4,5,2,2,1,6,7,8,1,2,7,5),
  var3= c(4,2,7,5,3,4,5,1,2,3,4,5,7,4,3),
  demographics= c("C","C","C","C","C","C","C","C","C","C","C","C","C","C","C"))

df2 <- data.frame(
  Day = c("1","1","1","1","1","1","1","1","1","1","1","1","1","1","1")
  var1= c(7,4,2,4,6,8,3,4,5,3,6,6,1,2,3),
  var2= c(4,2,7,5,3,4,5,1,2,3,4,5,7,4,3),
  var3= c(1,2,4,6,8,3,1,2,3,4,5,7,4,3,7))

df<- data.frame(
  CompanyID=c("Drinkers","Drinkers","Drinkers","Drinkers","Drinkers","Drinkers","Drinkers","Drinkers"
              ,"Drinkers","Drinkers", "Liquders","Liquders","Liquders","PelletCoffeeCo","PelletCoffeeCo", "NA",
              "NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA"),
  Email= c("john@coffee.com", "john@coffee.com","john@coffee.com","john@coffee.com", "john@coffee.com",
           "john@coffee.com", "john@coffee.com", "john@coffee.com", "john@coffee.com", "john@coffee.com",
           "george@liquid.com","george@liquid.com","george@liquid.com","stacy@pelletcoffee.com",
           "stacy@pelletcoffee.com","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA"),
  Day= c("1","2","3","4","5","6","7","8","9","10","1","2","3","1","2","1","1","1",
         "1","1","1","1","1","1","1","1","1","1","1","1"),
  var1= c(4,5,5,5,2,3,2,7,6,5,7,6,6,2,3,7,4,2,4,6,8,3,4,5,3,6,6,1,2,3),
  var2= c(4,5,5,4,5,2,2,1,6,7,8,1,2,7,5,4,2,7,5,3,4,5,1,2,3,4,5,7,4,3),
  var3= c(4,2,7,5,3,4,5,1,2,3,4,5,7,4,3,1,2,4,6,8,3,1,2,3,4,5,7,4,3,7),
  demographics= c("C","C","C","C","C","C","C","C","C","C","C","C","C","C","C","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA","NA"))