I have this following data which is a subset of large data, note that I have different number of id.
id <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3,