est_dist_censored: Function to return either normal, lognormal, gamma, or...

In jentjr/manager: Tools for plotting and analyzing environmental data

Function to return either normal, lognormal, gamma, or non-parametric distribution estimate

Description

Function to return either normal, lognormal, gamma, or non-parametric distribution estimate

Usage

est_dist_censored(
  df,
  left_censored,
  method = "sf",
  alpha = 0.05,
  choices = c("norm", "lnorm"),
  combine_locations = FALSE,
  keep_data_object = FALSE
)

Arguments

df	groundwater data frame in tidy format
left_censored	logical vector.
method	default is "sf"
alpha	alpha
choices	vector of distributions to check. Default is c("norm, "lnorm")
combine_locations	TRUE/FALSE to estimate distribution by individual location, or grouped together. Default is FALSE.
keep_data_object	Default is FALSE

Examples

data("gw_data")
wells <- c("MW-1", "MW-2", "MW-3")
params <- c("Fluoride, total",
            "Arsenic, dissolved",
            "Nickel, dissolved")

gw_data <- gw_data %>%
to_censored() %>%
filter(location_id %in% wells, param_name %in% params) %>%
percent_lt() %>%
filter(percent_lt >=15, percent_lt <= 50)

gw_data %>%
group_by(location_id, param_name, default_unit) %>%
est_dist_censored(., alpha = 0.05, keep_data_object = FALSE)


gw_data %>%
group_by(param_name, default_unit) %>%
est_dist_censored(., combine_locations = TRUE, keep_data_object = FALSE)

jentjr/manager documentation built on Jan. 12, 2024, 9:45 p.m.