nsec {bayesnec}R Documentation

Extracts the predicted NSEC value as desired from an object of class bayesnecfit or bayesmanecfit.

Description

Extracts the predicted NSEC value as desired from an object of class bayesnecfit or bayesmanecfit.

Usage

nsec(
  object,
  sig_val = 0.01,
  precision = 1000,
  posterior = FALSE,
  x_range = NA,
  hormesis_def = "control",
  xform = NA,
  prob_vals = c(0.5, 0.025, 0.975)
)

Arguments

object

An object of class bayesnecfit or bayesmanecfit returned by bnec.

sig_val

Probability value to use as the lower quantile to test significance of the predicted posterior values. against the lowest observed concentration (assumed to be the control), to estimate NEC as an interpolated NOEC value from smooth ECx curves.

precision

The number of unique x values over which to find NSEC - large values will make the NSEC estimate more precise.

posterior

A logical value indicating if the full posterior sample of calculated NSEC values should be returned instead of just the median and 95 credible intervals.

x_range

A range of x values over which to consider extracting NSEC.

hormesis_def

A character vector, taking values of "max" or "control". See Details.

xform

A function to apply to the returned estimated concentration values.

prob_vals

A vector indicating the probability values over which to return the estimated NSEC value. Defaults to 0.5 (median) and 0.025 and 0.975 (95 percent credible intervals).

Details

For hormesis_def, if "max", then NSEC values are calculated as a decline from the maximum estimates (i.e. the peak at nec); if "control", then ECx values are calculated relative to the control, which is assumed to be the lowest observed concentration.

Value

A vector containing the estimated NSEC value, including upper and lower 95% credible interval bounds.

See Also

bnec

Examples


library(bayesnec)

data(manec_example)
nsec(manec_example)



[Package bayesnec version 2.0.2 Index]