This function performs the Bayesian Sign test

bayesianSign.test(
  x,
  y = NULL,
  s = 1,
  z_0 = 0,
  rope.min = -0.01,
  rope.max = 0.01,
  weights = c(s/2, rep(1, length(x))),
  n.samples = 1e+05
)

Arguments

x

First vector of observations

y

Second vector of observations

s

Prior pseudo-observation probabilitie

z_0

Prior pseudo-observation

rope.min

Inferior limit of the rope considered

rope.max

Superior limit of the rope considered

weights

A priori weights

n.samples

Number of samples of the distributio

Value

List with probabilities for each region and a sample of posterior distribution.

Examples

bs <- bayesianSign.test(results$random.forest, results$KNN) bs.stronger.prior <- bayesianSign.test(results$random.forest, results$KNN, s=3, z_0 = 0.5)