Compute probability of observations across many potential vectors This will try (1/stepsize)^nflips possible vectors, computing the probability of the observation for each

try_many_vectors(
  nheads,
  nflips,
  number_samples = 1000,
  stopval = 1e-05,
  log = FALSE,
  possibilities = get_possibilities(nheads, nflips)
)

Arguments

nheads

Number of heads

nflips

Total number of flips (heads and tails)

number_samples

How many vectors to sample

stopval

How large a difference in probability is considered close enough between the flat model and others

log

If TRUE return log transformed probabilities.

possibilities

All possible sequences of flips that lead to the observed number of heads

Value

The likelihood of the data (or log likelihood if log=TRUE)