powered by
Quality control using SNPRelate functions.
snprelate_qc( gdata, samples_nas = 0.03, ibs = 0.99, keep_ids = NULL, snps_nas = 0.01, maf = 0.05, tagsnp = 0.8, n_cores = 1 )
List of gdata, Genotype data object, and df_qc, QC info data frame
Genotype data object
NA threshold for samples, default 3 pct
Samples identity by state threshold, default 99 pct
Samples ids to keep even if IBS is higher than threshold. Used for monozygotic twins.
NA threshold for SNPs, default 1 pct
Minor allele frequency threshold, default 5 pct
TagSNP r2 correlation threshold, default 0.8
Number of cores
library(snplinkage) gds_path <- save_hgdp_as_gds() gdata <- load_gds_as_genotype_data(gds_path) qc <- snprelate_qc(gdata, tagsnp = .99)
Run the code above in your browser using DataLab