find.outliers function

This identifies all of the SNPs outside of the smoothed quantiles in the dataset.

find.outliers(df, boot.out, ci.df = NULL, file.name = NULL)

Arguments

• df: Provide the dataframe with Ht and Fst values.
• boot.out: Bootstrap output. You must provide this.
• ci.df: List of confidence intervals. You may provide this in addition to bootstrap output to save a small amount of time.
• file.name: You may provide a file name to output the outliers to a csv file. Otherwise, the function will only return the outliers.

Returns

• out: A list of the outlier loci

Examples

## Not run:

    gfile<-system.file("extdata", "example.genepop.txt",package = 'fsthet')
    gpop<-my.read.gpop(gfile)
    fsts<-calc.actual.fst(gpop)
    boot.out<-as.data.frame(t(replicate(10, fst.boot(gpop))))
    outliers<-find.outliers(fsts,boot.out)
## End(Not run)