#Using psych::describeBy function to compare means(etc) by group

#Initialize an empty data frame

Wk01<-data.frame(matrix(ncol = 0, nrow = 20))

Wk01$grp<-c("M","M","M","M","M","M","M","M","M","M","F","F","F","F","F","F","F","F","F","F")
Wk01$score<-c(10,3,3,7,3,3,8,5,1,8,10,13,6,7,8,15,11,9,8,13)

library(psych)

psych::describeBy(Wk01$score,Wk01$grp,type=3,mat=TRUE)

##     item group1 vars  n mean       sd median trimmed    mad min max range
## X11    1      F    1 10 10.0 2.943920    9.5   9.875 2.9652   6  15     9
## X12    2      M    1 10  5.1 2.960856    4.0   5.000 2.9652   1  10     9
##          skew  kurtosis        se
## X11 0.2821983 -1.458284 0.9309493
## X12 0.2593537 -1.579995 0.9363048