R package for computing Polygenic Risk Scores (PRS) from dosages (DS). Supports VCF and BCF input formats.
Clone repository using git clone
and run the following command:
R CMD INSTALL Rprs
library(Rprs)
p <- prs_pt("weights.txt", "dosages.vcf.gz", "samples.txt", pvalues = c(1.0, 5e-8))
write.table(p, col.names = TRUE, row.names = FALSE, sep = "\t", quote = F)
weight_col
function argument. Default is WEIGHT_
.library(Rprs)
p <- prs_gw("weights.txt", "dosages.vcf.gz", "samples.txt", weight_col = "LDpred_")
write.table(p, col.names = TRUE, row.names = FALSE, sep = "\t", quote = F)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.