readRefTable: Read a reference table simulated from DIYABC

In abcrf: Approximate Bayesian Computation via Random Forests

Read a reference table simulated from DIYABC

Description

readRefTable reads a reference table simulated from DIYABC thanks to a .bin and a .txt file, respectively containing the reference table and its header.

Usage

readRefTable(filename = "reftable.bin", header = "header.txt", N = 0)

Arguments

filename	a .bin file from DIYABC containing the reference table.
header	a .txt file containing the header of the reference table in filename.
N	an integer indicating the number of observations to extract from the reference table. The default is 0 indicating that the whole reference table is used. Warning: if N is specified, nrecscen returns is obsolete.

Value

A list with the following components:

nrec	number of individuals of the reference table,
nscen	number of scenarios in the reference table,
nrecscen	number of individuals by scenario,
nparam	number of parameters by scenario,
scenarios	a vector of factor containing the scenario indices,
params	a matrix with the parameters,
stats	a matrix with the summary statistics.

References

Cornuet J.-M., Pudlo P., Veyssier J., Dehne-Garcia A., Gautier M., Leblois R., Marin J.-M. and Estoup A. (2014) DIYABC v2.0: a software to make Approximate Bayesian Computation inferences about population history using Single Nucleotide Polymorphism, DNA sequence and microsatellite data Bioinformatics doi: 10.1093/bioinformatics/btn514

abcrf documentation built on Aug. 9, 2022, 5:07 p.m.