library(hdrcde)
set.seed(3)
x = rnorm(100)
hdr(x, h = bw.nrd0(x))
$$hdr
$          [,1]      [,2]      [,3]      [,4]
$99% -1.7398509 1.7800070        NA        NA
$95% -1.4201507 1.5091742        NA        NA
$50% -0.8039541 0.6226473 0.6884934 0.7054537

$$mode
$[1] 0.1205937

#$falpha
#        1%         5%        50%
#0.05522019 0.13245133 0.34655653