x<- c(93, 60, 55, 91, 107, 86, 113, 90, 97, 92, 74, 144, 93, 116, 90, 87, 125, 89, 84, 104, 92, 109, 77, 77, 81, 124, 112, 125, 105, 128, 103, 96, 98, 91)

length(x)

mean(x)

median(x)

sd(x)

library("DescTools")

MeanAD(x, center = median)

bw <-(34**-0.2)

bw

kde = kdensity(x, bw, kernel = "uniform", normalized = FALSE)

kde(93)