How do I find the pooled standard error of the mean in R?

Question

I have the following result from doing some descriptive statistics on my data:

> psych::describeBy(mydata_young$length, mydata_young$Strain)

 Descriptive