You know how skew and kurtosis are calculated (or look it up, e.g. in Wikipedia). So you can write your functions.

Note that different formulas for computing skewness and kurtosis can be used. See, for example, the options "type" of the function kurtosis() shown by using ?kurtosis() after loading the package e1071 by using library(e1071).

You can use package 'moments'. Check following link to get pdf of this package

http://cran.r-project.org/web/packages/moments/moments.pdf

Good luck!

You can also use the library "propagate". Install it in R, load it and use the functions kurtosis() and skewness() as follows:

> library (propagate)

> X skewness(X)

> kurtosis(X)

The psych package has a function describe(), which gives quite a few descriptives including skewness and kurtosis.

I believe it uses the methods similar to e1071, but you might want to confirm that.

Here is a good example;

#Simulation

n.sample

This link can help you

https://alstatr.blogspot.com/2013/06/measures-of-skewness-and-kurtosis.html