Here my data
mydat=structure(list(id = c(15010124001, 15010153006, 15010169005, 15010228019, 15010229028, 6010001012, 6010012023, 6010014015, 6010015008, 60