轻松使用聚合

aggregate(cover~species+plotID, data=df_original, FUN=sum) 

使用 data.table

as.data.table(df_original)[, sum(cover), by = .(plotID, species)]