R/R/MAF.R

In qiuanzhu/xlink: Genetic Association Models for X-Chromosome SNPS on Continuous, Binary and Survival Outcomes

#' Select SNP by MAF.
#' 
#' \code{MAF} returns SNPs with higher MAF than default value. 
#' 
#' @param snp SNP name.
#' @param gender gender vairable.
#' @param male male information.
#' @param MAF_v minimum MAF value.
MAF <- function(snp, gender, male, MAF_v, data) {
    data <- stats::na.omit(data[, which(colnames(data) %in% c(snp, gender))])
    n <- length(data$snp)
    ind <- (data$gender == male)
    MAF_R <- (2 * sum(data$snp[ind] != 0) + 2 * sum(data$snp[!ind] == 2) + 1 * sum(data$snp[!ind] == 1))/(2 * n)
    if (MAF_R >= MAF_v) {
        return(c(snp, MAF_R))
    } else {
        return(NULL)
    }
}