Description Usage Arguments Author(s) References Examples
Compares the QST of a single phenotypic trait to the mean FST of series of marker loci. It calculates the distribution of QST  FST under a model assuming neutrality of both the phenotypic trait and the genetic markers from which FST is estimated. Returns the simulated estimates of Qst  Fst under neutrality following the procedure described in Gilbert and Whitlock (2014) and Whitlock & Guillaume (2009). Also returns the simulated estimates of Fst and Qst used to compute the null distribution.
1 2  QstFstComp(fst.dat, qst.dat, numpops, nsim = 1000, AFLP = FALSE,
breeding.design, dam.offspring.relatedness = 0.25, output = "concise")

fst.dat 
A data frame with the first column indicating population of origin and
the following columns representing genotypes at loci; see the
README https://github.com/kjgilbert/QstFstComp/blob/master/README.md for further description.
If using AFLPs, this is a data frame of qhat values, with pops in columns, loci in rows
and the corresponding qhat variances in the following columns, and 
qst.dat 
the input table of the breeding data

numpops 
number of populations in the sample 
nsim 
number of simulation replicates to perform to create the null distributions and bootstraps @param AFLP whether or not to use AFLP data @param breeding.design the breeding design used when collecting the trait data There are two options for breeding design:
@param dam.offspring.relatedness relatedness between offspring in the dam model, default is 1/4, i.e. halfsib @param output whether to output full, concise, or without writing out vector of resampled QF values, see details below @return Returns either a concise list of a subset of results or a full list with all possible results. Both output options write the vector of QstFst values to a text file unless "_nowrite" is appended to the option, e.g. "concise_nowrite" returns only the concise output without writing output to a text file. Concise list returns (default)
Full list returns
Be mindful of the fact that with a small number of loci, bootstrapped confidence intervals of Fst can be less accurate. 
Kimberly J Gilbert & Michael C Whitlock
Gilbert KJ and MC Whitlock (2015) Qst Fst comparisons with unbalanced halfsib designs. Molecular Ecology Resources, 15(2), 262267.
Whitlock MC and F Guillaume (2009) Testing for spatially divergent selection: Comparing Qst to Fst. Genetics, 183:10551063.
1 2 3 4 5 6 7 8  ## using balanced halfsib sire trait data and biallelic marker data
data(hssire) # trait data
data(biallelic) # marker data
QstFstComp(biallelic, hssire, numpops=15, nsim=100, breeding.design="half.sib.sire", output="full")
data(hsdam)
data(aflp)
QstFstComp(aflp, hsdam, numpops=15, nsim=100, AFLP=TRUE, breeding.design="half.sib.dam", output="concise")

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.