Compute bootstrapped parameters and their related indices such as Confidence Intervals (CI) and p-values.

## Usage

```
bootstrap_parameters(
model,
iterations = 1000,
centrality = "median",
ci = 0.95,
ci_method = "quantile",
test = "p-value",
...
)
```

## Arguments

- model
Statistical model.

- iterations
The number of draws to simulate/bootstrap.

- centrality
The point-estimates (centrality indices) to compute. Character (vector) or list with one or more of these options:

`"median"`

,`"mean"`

,`"MAP"`

or`"all"`

.- ci
Value or vector of probability of the CI (between 0 and 1) to be estimated. Default to

`.95`

(`95%`

).- ci_method
The type of index used for Credible Interval. Can be

`"ETI"`

(default, see`eti()`

),`"HDI"`

(see`hdi()`

),`"BCI"`

(see`bci()`

),`"SPI"`

(see`spi()`

), or`"SI"`

(see`si()`

).- test
The indices to compute. Character (vector) with one or more of these options:

`"p-value"`

(or`"p"`

),`"p_direction"`

(or`"pd"`

),`"rope"`

,`"p_map"`

,`"equivalence_test"`

(or`"equitest"`

),`"bayesfactor"`

(or`"bf"`

) or`"all"`

to compute all tests. For each "test", the corresponding**bayestestR**function is called (e.g.`bayestestR::rope()`

or`bayestestR::p_direction()`

) and its results included in the summary output.- ...
Arguments passed to or from other methods.

## Details

This function first calls `bootstrap_model()`

to generate
bootstrapped coefficients. The resulting replicated for each coefficient
are treated as "distribution", and is passed to `bayestestR::describe_posterior()`

to calculate the related indices defined in the `"test"`

argument.

Note that that p-values returned here are estimated under the assumption of
*translation equivariance*: that shape of the sampling distribution is
unaffected by the null being true or not. If this assumption does not hold,
p-values can be biased, and it is suggested to use proper permutation tests
to obtain non-parametric p-values.

## Using with **emmeans**

The output can be passed directly to the various functions from the
**emmeans** package, to obtain bootstrapped estimates, contrasts, simple
slopes, etc. and their confidence intervals. These can then be passed to
`model_parameter()`

to obtain standard errors, p-values, etc. (see
example).

Note that that p-values returned here are estimated under the assumption of
*translation equivariance*: that shape of the sampling distribution is
unaffected by the null being true or not. If this assumption does not hold,
p-values can be biased, and it is suggested to use proper permutation tests
to obtain non-parametric p-values.