R

aggregate(x = xdata$species, by = list(xdata$id), FUN = paste, collapse = ",")
aggregate(data = xdata, species ~ id, FUN = paste, collapse = ",")