compare_posterior( x, comparison = "nec", ecx_val = 10, type = "absolute", hormesis_def = "control", sig_val = 0.01, precision, x_range = NA, make_newdata = TRUE, ... )
The posterior predictions to compare, takes values of "nec", "nsec", "ecx" or "fitted".
The desired percentage effect value. This must be a value between 1 and 99 (for type = "relative" and "absolute"), defaults to 10.
character vector, taking values of
"relative", "absolute" (the default) or "direct". See Details.
character vector, taking values
of "max" or "control". See Details.
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.
The number of unique x values over which to find ECx -- large values will make the ECx estimate more precise.
A range of x values over which to consider extracting ECx.
Only used if
comparison = "fitted". Should the
user allow the package to create
newdata for predictions?
If so, arguments
x_range will be used. Defaults
to TRUE. See details.
Further arguments that control posterior predictions via
type "relative" is calculated as the percentage decrease
from the maximum predicted value of the response (top) to the minimum
predicted value of the response. Type "absolute" (the default) is
calculated as the percentage decrease from the maximum value of the
response (top) to 0 (or bot for a 4 parameter model fit). Type "direct"
provides a direct estimate of the x value for a given y.
Note that for the current version, ECx for an "nechorme" (NEC Hormesis)
model is estimated at a percent decline from the control.
hormesis_def, if "max", then ECx or NSEC values -- i.e.,
depending on argument
comparison -- are calculated
as a decline from the maximum estimates (i.e. the peak at NEC);
if "control", then ECx or NSEC values are calculated relative to the
control, which is assumed to be the lowest observed concentration.
make_newdata is only used if
comparison = "fitted". It is relevant to those who want the package
to create a data.frame from which to make predictions. This is done via
bnec_newdata and uses arguments
make_newdata = FALSE and no additional
newdata argument is provided (via
...), then the predictions
are made for the raw data. Else, to generate predictions for a specific
user-specific data.frame, set
make_newdata = FALSE and provide
an additional data.frame via the
newdata argument. For guidance
on how to structure
newdata, see for example