generate_initialPop: Calculation of population structure at t1

Description Usage Arguments Details Author(s)

Description

Calculates possible population structures at time t1 so that the observed sample could have been drawn at random from these structures.

Usage

1
generate_initialPop<-function(iniSample, n_target_sim, alphaDirichlet)

Arguments

iniSample

vector containing the absolute frequencies of the observed cultural variant types at t1

n_target_sim

number of population structures to be generated

alphaDirichlet

vector of size (k+1), where k stands for the number of variant types present at t1, containing the parameters of the Dirichlet distribution

Details

This function will generate a vector of size (k+1) containing relative frequencies of the variant types present at t1 in the population. The variant type (k+1) contains variants not observed in the sample. The estimation is based on the Dirichlet distribution and is independent of total population size. The resulting population structure has the property that the observed sample at t1 could have been drawn at random from it. Importantly, iniSample has to contain absolute frequencies.

Author(s)

Clementine Straub


annenna4/culturechange documentation built on May 30, 2019, 6:54 p.m.