code for GRUP data project, some plot with data subsetting