## Get levels of species list
taxlist::levels(Easplist)
## Add aggregate as new taxonomic level
levels(Easplist) <- c(
"form", "variety", "subspecies", "species",
"complex", "aggregate", "genus", "family"
)
summary(Easplist)
Run the code above in your browser using DataLab