diff --git a/R/applyQC.R b/R/applyQC.R index 268ca70..1c1e08c 100755 --- a/R/applyQC.R +++ b/R/applyQC.R @@ -266,7 +266,7 @@ cleanData <- function(indir, name, qcdir=indir, if (is.null(hweTh)) { stop("filterHWE is TRUE but hweTh not specified") } else { - hwe <- c("--hwe", hweTh) + hwe <- c("--hwe midp", hweTh) } } if (filterMAF) { diff --git a/R/markerQC.R b/R/markerQC.R index ca18f28..a960a7b 100755 --- a/R/markerQC.R +++ b/R/markerQC.R @@ -436,7 +436,7 @@ check_hwe <- function(indir, name, qcdir=indir, hweTh=1e-5, interactive=FALSE, } suffix <- "" sys::exec_wait(path2plink, - args=c("--bfile", prefix, "--hardy", "--out", + args=c("--bfile", prefix, "--hardy midp", "--out", paste(out, suffix, sep="")), std_out=showPlinkOutput, std_err=showPlinkOutput) } else {