diff --git a/OysterHatcheryMicrobiome_all16SFigures.Rmd b/OysterHatcheryMicrobiome_all16SFigures.Rmd
index 0e70c7c..e50bdf3 100644
--- a/OysterHatcheryMicrobiome_all16SFigures.Rmd
+++ b/OysterHatcheryMicrobiome_all16SFigures.Rmd
@@ -586,7 +586,7 @@ allwaterorder<-as.data.frame(c(meta_table_wo, abund_table_wo))
 allwaterorder$Sample<-allwaterorder$taxa
 
 # group data and calculate means of bacillales, oceanospirillales, vibrionales
-grouped_orderdata <- plyr::group_by(allwaterorder, Trial, Day, Treatment)
+grouped_orderdata <- group_by(allwaterorder, Trial, Day, Treatment)
 data_means <- dplyr::summarise(grouped_orderdata,
                          meanoc = mean(Bacteria.Proteobacteria.Gammaproteobacteria.Oceanospirillales, na.rm=TRUE), 
                          stdevoc = sd(Bacteria.Proteobacteria.Gammaproteobacteria.Oceanospirillales, na.rm=TRUE),