compareTwoDistributions: Function to compare two distributions and test a hypothesis,...

In clquezada/tRophicPosition: Bayesian Trophic Position Estimation with Stable Isotopes

Function to compare two distributions and test a hypothesis, in a Bayesian context

Description

Function to compare two distributions and test a hypothesis, in a Bayesian context

Usage

compareTwoDistributions(
  dist1 = NULL,
  dist2 = NULL,
  test = "<=",
  sample = NULL,
  round = 3,
  ...
)

Arguments

dist1	A collection of numerical values (posterior distribution).
dist2	A collection of numerical values (posterior distribution).
test	A logical operator which states what to test for. Might be "<", "<=", ">" or ">=".
sample	If sample is numeric, it will take 'sample' elements of each of the distributions.
round	integer to indicate number of decimals kept.
...	extra arguments are passed to compareTwoDistributions().

Value

probability given sum(dist1 "test" dist2) / length(dist1)

Examples

a <- rnorm(100, 2, 0.1)
b <- rnorm(100, 1.8, 0.1)
compareTwoDistributions(a, b, test = ">=")
compareTwoDistributions(a, b, test = "bhatt")

clquezada/tRophicPosition documentation built on Jan. 4, 2023, 12:31 p.m.