compareestimators: Compare relatedness estimators.

In related: related: an R package for analyzing pairwise relatedness data based on codominant molecular markers

Description

Generates simulated individuals of known relatedness, based on observed allele frequencies, then estimates relatedness using 4 methods (lynchli, lynchrd, quellergt, wang) and creates a boxplot of results so that the performance of the different estimators can easily be compared.

Usage

compareestimators(filename, ninds)

Arguments

filename	An object generated by the readgenotypedata function, containing genotype and allele frequency data.
ninds	A whole number indicating the number of pairs of individuals of known relatedness you want to generate for each relatedness value (e.g., entering 100 would result in 100 parent-offspring pairs, 100 full-sib pairs, 100 half-sib pairs, and 100 unrelated pairs.)

Value

A figure containing box plots comparing the values obtained from each estimator for each known relatedness value.

Author(s)

Tim Frasier <[e-mail]>

See Also

readgenotypedata, familysim, cleanuprvals

Examples

	## Not run: 
		#---Load data set---#
		data(GenotypeData)
		input <- readgenotypedata(GenotypeData)

		#---Compare relatedness estimators---#
		compareestimators(input, 100)
	
## End(Not run)	

related documentation built on May 2, 2019, 6:49 p.m.