plot.egger {GLIDE} | R Documentation |
It is used to draw MR-Egger regression plot.
## S3 method for class 'egger' plot(x,qcutoff=0.2,xlab="Genetic association with the exposure", ylab="Genetic association with the outcome",...)
x |
x is a dataframe, it is the outcome of the glide function. |
qcutoff |
qcutoff is the q-value cutoff. |
xlab |
xlab is the default x label. |
ylab |
ylab is the default y label. |
... |
the rest of plot arguments. |
James Y. Dai and X. Wang
## Not run: ##first generate output using the glide function #load an example data data(simdata) simdat=simdata$simdat coeff=simdata$coeff formula=as.formula("outcome~age+sex+pc1+pc2+pc3") genotype_columns=which(grepl("^SNP",colnames(simdat))) #run glide out=glide(formula=formula,exposure_coeff=coeff,genotype_columns,data=simdat,np=100000, qcutoff=0.2,parallel=TRUE,corenumber=1,verbose=TRUE) ##draw the Egger plot plot.egger(out) ## End(Not run)