This is a complex follow up to this question Split a dataframe into a list of nested data frames and matrices I now split the original data frame with 2 grouping variables.