I finished the aov and post-hoc test. And I need the data of means and Standard deviation for reporting. While I can locate the means by model.tables(nes2000.aov,type='means"). I don't know how to calculate standard deviation in R.
I tried the below command I found online But it pops on error "Error in aggregate.data.frame(lhs, mf[-1L], FUN = FUN, ...) . I don't know what went wrong.
argument "FUN" is missing, with no default