plotZ: plotZ

In jamesliley/subtest: Test whether two putative subgroups of a disease phenotype show evidence of differential genetic basis.

Description

Plot Z_a and Z_d scores efficiently, blocking out dense region near origin.

Usage

plotZ(z_ad, rlim = quantile(z_ad[, 2], 0.99), col = rep("black",
  dim(z_ad)[1]), mcol = NULL, abs = (min(z_ad) < 0), over = FALSE, ...)

Arguments

rlim	block out circle of radius rlim; only plot points with Z_d^2 + Z_a^2 > rlim^2.
col	colour of points or vector describing colour of each point of Z
mcol	colour of middle blocked-out region; defaults to first element of col
abs	set to TRUE to plot only upper right quadrant
over	set to TRUE to plot over existing plot; FALSE to draw new plot.
Z	an n x 2 array; Z[i,1], Z[i,2] are the Z_d and Z_a scores respectively for the ith SNP

Author(s)

James Liley

Examples

nn=100000
Z=abs(rbind(rmnorm(0.8*nn,varcov=diag(2)), rmnorm(0.15*nn,varcov=rbind(c(1,0),c(0,2^2))), rmnorm(0.05*nn,varcov=rbind(c(3^2,2),c(2,4^2)))));
plotZ(Z,rlim=2,col="red",cex=0.5)

jamesliley/subtest documentation built on May 18, 2019, 11:21 a.m.