another practice

library(tseries)
library(ggplot2)
library(zoo)

## 
## Attaching package: 'zoo'

## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric

library(moments)

n.sample <- rnorm(n = 10000, mean = 55, sd = 4.5)

kurtosis(n.sample)

## [1] 3.030709

skewness(n.sample)

## [1] -0.01807928

data <- data.frame(n.sample)
ggplot(data, aes(x=n.sample), bandwidth=2)+ geom_histogram(colour='red')+ xlab(expression(bold('Simulated Samples')))+ ylab(expression(bold('Density')))

## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.