pivec for flexdog EM algorithm.The key idea here is choosing the pi's so that the two modes have equal probability.
initialize_pivec(
ploidy,
mode,
model = c("hw", "bb", "norm", "ash", "f1", "s1", "f1pp", "s1pp", "f1ppdr", "s1ppdr",
"flex", "uniform", "custom")
)The ploidy of the species. Assumed to be the same for each individual.
The mode if model = "ash". If not provided,
flexdog will estimate the mode. This is the starting point of
the allele frequency if model = "hw". This should be
NULL for all other options of model.
What form should the prior (genotype distribution) take? See Details for possible values.
A vector of numerics. The initial value of pivec
used in flexdog_full.
flexdog for where this is used.