Skip to contents

Function used to specify univariate normal density estimator for get_hdr_1d() and layer functions (e.g. geom_hdr_rug()).

Usage

method_norm_1d()

Details

For more details on the use and implementation of the method_*_1d() functions, see vignette("method", "ggdensity").

Examples

# Normal estimators are useful when an assumption of normality is appropriate
df <- data.frame(x = rnorm(1e3))

ggplot(df, aes(x)) +
  geom_hdr_rug(method = method_norm_1d()) +
  geom_density()


# Can also be used with `get_hdr_1d()` for numerical summary of HDRs
res <- get_hdr_1d(df$x, method = method_norm_1d())
str(res)
#> List of 3
#>  $ df_est:'data.frame':	512 obs. of  4 variables:
#>   ..$ x               : num [1:512] -3.54 -3.53 -3.51 -3.5 -3.49 ...
#>   ..$ fhat            : num [1:512] 0.00102 0.00106 0.0011 0.00115 0.0012 ...
#>   ..$ fhat_discretized: num [1:512] 1.28e-05 1.33e-05 1.39e-05 1.45e-05 1.51e-05 ...
#>   ..$ hdr             : num [1:512] 1 1 1 1 1 1 1 1 1 1 ...
#>  $ breaks: Named num [1:5] 0.0178 0.0595 0.1721 0.3089 Inf
#>   ..- attr(*, "names")= chr [1:5] "99%" "95%" "80%" "50%" ...
#>  $ data  :'data.frame':	1000 obs. of  2 variables:
#>   ..$ x             : num [1:1000] -0.8861 -1.9223 1.6197 0.5193 -0.0558 ...
#>   ..$ hdr_membership: num [1:1000] 0.8 0.95 0.95 0.5 0.5 0.5 0.5 0.95 0.99 0.5 ...