PRS_PGx_CT: Construct PGx PRS unadjusted or using clumping and...

In PRSPGx: Construct PGx PRS

Construct PGx PRS unadjusted or using clumping and thresholding

Description

Shrink prognostic and predictive effect sizes simutaneously by 2-df (main and interaction) p-value cutoff (PRS-PGx-CT turns out to be PRS-PGx-Unadj when setting p-value cutoff = 1)

Usage

PRS_PGx_CT(
  PGx_GWAS,
  G_reference,
  pcutoff = 1e-04,
  clumping = TRUE,
  p1 = 1e-04,
  d1 = 250000,
  r1 = 0.8
)

Arguments

PGx_GWAS	a numeric matrix containing PGx GWAS summary statistics, including SNP ID, MAF, position, β, α, 2-df p-value, and N
G_reference	a numeric matrix containing the individual-level genotype information from the reference panel (e.g., 1KG)
pcutoff	a numeric value indicating the p-value cutoff
clumping	a logical flag indicating should clumping be performed
p1	a numeric value indicating p-value threshold to decide flag SNPs in clumping
d1	a numeric value indicating window size in clumping
r1	a numeric value indicating correlation in clumping

Details

PRS-PGx-CT only needs PGx summary statistics

Value

A numeric list, the first sublist contains estimated prognostic effect sizes, the second sublist contains estimated predictive effect sizes, the third sublist contains 2-df p-values

Author(s)

Song Zhai

References

Zhai, S., Zhang, H., Mehrotra, D.V. & Shen, J. Paradigm Shift from Disease PRS to PGx PRS for Drug Response Prediction using PRS-PGx Methods (submitted).

Examples

data(PRSPGx.example); attach(PRSPGx.example)
coef_est <- PRS_PGx_CT(PGx_GWAS, G_reference, pcutoff = 0.01, clumping = TRUE)
summary(coef_est$coef.G)
summary(coef_est$coef.TG)

PRSPGx documentation built on July 20, 2022, 5:07 p.m.