per_locus_means_and_vars: Compute the mean and variance of the single-locus genotype...

View source: R/per_locus_means_and_vars.R

per_locus_means_and_varsR Documentation

Compute the mean and variance of the single-locus genotype likelihoods for each collection

Description

This assumes that you have compiled params for a reference data set and then it just calls rcpp_per_locus and then summarizes the results.

Usage

per_locus_means_and_vars(par_list)

Arguments

par_list

genetic data converted to the param_list format by tcf2param_list. This should be include genotypes only for the reference individuals.

Value

Returns a list with two components, mean and var, each one a matrix that has C (number of collections) rows and L (number of loci) columns, giving the mean (or variance) of the genotype likelihoods in the individuals in that collection at that locus.


benmoran11/rubias documentation built on Feb. 1, 2024, 10:52 p.m.