Learn R Programming

rubias (version 0.3.4)

per_locus_means_and_vars: 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)

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.

Arguments

par_list

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