man_plot_prob: Create Manhattan plot for posterior probability based models

In chenchunyu88/BATools: BATools: An R package that extends 'Bayesian Alphabet' for Whole Genome Prediction and Genome-wide association study

Description

Create Manhattan plot for posterior probability based models

Usage

man_plot_prob(ba, type = c("SNP", "Win"), col = c("black", "red"), ...)

Arguments

ba	an object of the class ba generated by baFit
type	can be c("SNP","Win") for single SNP or window based Manhattan plot
col	color for the Manhattan plot
...	other plot options

Value

Null. This function only prints to the screen.

Examples

## Not run: 
 rm(list=ls())
 library(BATools)
 data("Pig")
 #Standardize genotype matrix
 #or set your own starting values using 
 #init=list(df=5,scale=0.01,pi=1) 
 run_para=list(niter=2000,burnIn=1000,skip=10)
 print_mcmc=list(piter=500)
 update_para=list(df=FALSE,scale=TRUE,pi=F)
 op<-create.options(model="SSVS",method="MCMC",seed=1,priors=NULL,init=init,
                  update_para=update_para,run_para=run_para,save.at="SSVS",print_mcmc=print_mcmc)
 SSVS<-baFit(driploss~sex,data=PigPheno,geno=geno ,genoid = ~id,options = op,map=PigMap,GWA="Win")
 par(mfrow=c(1,2))
 man_plot_prob(SSVS)
 man_plot_prob(SSVS,type="Win")

## End(Not run)

chenchunyu88/BATools documentation built on May 19, 2019, 8:21 a.m.