pipeline/R/3_subset_zmat.R

library(dplyr)
library(purrr)

args <- commandArgs(trailingOnly=TRUE)
normbeta <- readRDS(args[1])
keep_snps <- readRDS(args[2])
out_data <- args[3]


X <- normbeta %>%
      filter(snp %in% keep_snps$snp)

saveRDS(X, file=out_data)
jean997/mrScan documentation built on Dec. 20, 2024, 3:39 a.m.