Description Usage Arguments Details Author(s) See Also Examples
Perform a priori goodness of fit using the two first components obtained with PCA.
1 
target 
a data frame or vector of the observed summary statistic. 
sumstat 
a matrix or data frame of the simulated summary statistics. 
index 
a vector of models names. It must be character and have the same
length as the number of row in 
cprob 

xlim, ylim 
optional, numeric vectors of length 2, giving the x and y coordinates ranges. 
... 
other parameters passed to 
The function performs PCA using the a priori simulated summary
statistics. It displays envelopes containing 1hprob
percent of
the simulations.
The projection of the observed summary statistics is displayed in
order to check if they are contained or not in the envelopes.
If the projection lies outside the envelope of a given model, it is
an indication of poor fit.
Louisiane Lemaire and Michael Blum
abc
, plot.gfit
, summary.gfit
,
gfit
1 2 3 4 5 6 7 8 9 10 11 12 13 14  ## human demographic history
require(abc.data)
data(human)
## five R objects are loaded. See ?human and vignette("abc") for details.
## Perform a priori goodness of fit for 3 different demographic models
## The envelopes containing 90% of the simulations are displayed.
## For the European data, a reasonable fit is only provided by the
## bottleneck model.
## The number of simulations is reduced to improve speed (do not do that
## with your own data)
index<c(1:5000,50001:55000,100001:105000)
gfitpca(target=stat.voight["italian",], sumstat=stat.3pops.sim[index,],
index=models[index], cprob=0.1)

Loading required package: abc.data
Loading required package: nnet
Loading required package: quantreg
Loading required package: SparseM
Attaching package: 'SparseM'
The following object is masked from 'package:base':
backsolve
Loading required package: MASS
Loading required package: locfit
locfit 1.59.1 20130322
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.