R/combine_data_spc_function.R

In tydarnell/GeneticMediation: An R Package for Causal Mediation Analysis of ROSMAP Data

#' Combine Data SNP-Peak-Clinical
#'
#' @description SNP-Peak-Clinical data: combine projid, SNPs, PCs, membership,
#'   age, gender, peaks, and outcome data for a chromosome
#' @param chr chromosome name
#' @keywords dataprep spc
#' @export

combine_data_spc <- function(chr) {
    load(paste("data/mediator/mediation_", chr, "_new.RData", sep = ""))
    load(paste("data/peak.data/peak.data_", chr, ".RData", sep = ""))
    SNP.cov.data = merge(SNP.pc.data, phenotype.data[, c("projid", "gender")], by = "projid", all = F)
    disease.data = phenotype.data[, c("projid", "plaq_n_sqr")]
    
    med.data = Reduce(function(x, y) merge(x, y, all = F, by = "projid"), list(SNP.cov.data, peak.data, 
        disease.data))
    med.data = unique(med.data)
    med.data$membership = factor(med.data$membership)
    med.data$gender = factor(med.data$gender)
    med.data
}