# data
data(data.JSM)
cat(data.JSM$genepop, file="JSM_genepop.txt", sep="\n")
# remove markers designated by their names
clip.genepop("JSM_genepop.txt", "JSM_genepop_clipped_name.txt",
remove.list=c("Sni21","Sni26"))
# remove markers with high major allele frequencies (in this example, > 0.5)
clip.genepop("JSM_genepop.txt", "JSM_genepop_clipped_maf.txt",
major.af=0.5)
# remove markers both by their names and by major allele frequencies
clip.genepop("JSM_genepop.txt", "JSM_genepop_clipped_both.txt",
remove.list=c("Sni21","Sni26"), major.af=0.5)
Run the code above in your browser using DataLab