xQTLbiolinks (version 1.6.2)

xQTLvisual_PZPlot: Compare P-values reported to P-values calculated from Z statistics derived from the reported beta and standard error.

Description

Compare P-values reported to P-values calculated from Z statistics derived from the reported beta and standard error.

Usage

xQTLvisual_PZPlot(
  summaryDT,
  binCutLogP = 4,
  binNumber = 2000,
  distribution_func = "pnorm"
)

Value

a list containing a data.frame of estimated pvalues and A ggplot2 object

Arguments

summaryDT

A data.frame with three cols: pval, beta, se.

binCutLogP

To speed up the rendering process of the plot for tens of millions of GWAS variants, variants with a p-value below a specified threshold (binCutLogP) are randomly sampled for display.

binNumber

The number of points randomly selected for plotting.

distribution_func

"pnorm"(default) or "pchisq"

Examples

# \donttest{
url1 <- "http://bioinfo.szbl.ac.cn/xQTL_biolinks/xqtl_data/gwasDFsub_MMP7.txt"
sumDT <- data.table::fread(url1, sep="\t")
xQTLvisual_PZPlot(sumDT[,.(pValue, beta, se)], distribution_func="pchisq")
# }