Compare fitted f2, f3 and f4 f-statistics of an admixture graph with estimated ones
compare.fitted.fstats(fstats, fitted.graph, n.worst.stats = 5)
A matrix with 3 columns for each test (row names of the matrix corresponding to the test):
The estimated f-statistics (mean across block-Jackknife samples)
The fitted f-statistics (obtained from the fitted grah parameters
A Z-score measuring the deviation of the fitted values from the estimated values in units of standard errors (i.e., Z=(fitted.value-target.value)/se(target.value))
Object of class fstats containing estimates of fstats (as obtained with compute.fstats)
Object of class fitted graph (as obtained with fit.graph function).
The number of worst statistics to be displayed in the terminal
Compare fitted and estimated f-statistics may allow identifying problematic edges on the graph.
See compute.fstats
and fit.graph