axials_combine {kindisperse}

Combine axial distributions to produce a mixed distribution

Description

combines axial distributions to produce a mixed distribution. This is useful in settings where you have two separate distributions (e.g. FS & HS) with their own axial values, but you want to average them appropriately so that they can be compared to e.g. a mixed distribution of full & half cousins which cannot be distinguished via kinship determination methods and thus are best treated as an even mixture of the two categories. Different to adding dispersal events.

Usage

axials_combine(axvals)

Arguments

axvals

numeric. vector of axial distribution values from different kernels that are to be combined

Value

numeric Returns the axial value that results from combining the input axial values under an additive variance framework.

Examples

fullax <- axials(c(2, 4, 5), composite = 2)
halfax <- axials(c(6, 5, 7), composite = 2)
sibax <- axials_combine(c(fullax, halfax))