readRefTable {abcrf}

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