Add text to explain waht is going on.

Load some data

dat = read.csv("./irished2.csv")

head(dat)
##   X sex DVRT edlevel lvcert fathocc schltype  logDVRT
## 1 1   1  113       3      0      28        1 4.727388
## 2 2   1  101       1      0      28        9 4.615121
## 3 3   1  110       9      1      69        1 4.700480
## 4 4   1  121       5      0      57        1 4.795791
## 5 5   1   82       4      0      18        2 4.406719
## 6 6   1   85       4      0      28        2 4.442651

Calculate mean

mean(dat$DVRT)
## [1] 100.152

Make a histogram

hist(dat$DVRT, col = 'darkorange')