Hi, > result <- lapply(unique(dataset$ID), function(uid) dataset[dataset$ID > == uid, ]) Ok, I have the element result as a list > names(result) <- unique(dataset$ID) Nothing happens. I don't have any submatrix... Matteo