### Basic descriptives, except for mode (see below)

library(psych)

describe(dat1$V4)

##    vars   n   mean    sd median trimmed   mad min max range skew kurtosis  se
## X1    1 100 100.02 14.97   99.5   99.71 14.08  66 154    88 0.37     0.64 1.5

library(DescTools)

Mode(dat1$V4) #Note multiple modes with a frequency of 3 below

##  [1]  91  93  95  97  98  99 101 102 104 105 108
## attr(,"freq")
## [1] 3