| GWASTrack-class | R Documentation |
GWASTrack creates an IGV manhattan track GWAS data
GWASTrack( trackName, table, chrom.col, pos.col, pval.col, color = "darkBlue", trackHeight = 50, visibilityWindow = 1e+05 )
trackName |
A character string, used as track label by igv, we recommend unique names per track. |
table |
data.frame of 6 or more columns |
chrom.col |
numeric, the column number of the chromosome column |
pos.col |
numeric, the column number of the position column |
pval.col |
numeric, the column number of the GWAS pvalue colum |
color |
A css color name (e.g., "red" or "#FF0000" |
trackHeight |
track height, typically in range 20 (for annotations) and up to 1000 (for large sample vcf files) |
visibilityWindow |
Maximum window size in base pairs for which indexed annotations or variants are displayed. Defaults: 1 MB for variants, whole chromosome for other track types. |
A GWASTrack object
file <- system.file(package="igvR", "extdata", "gwas-5k.tsv")
tbl.gwas <- read.table(file, sep="\t", header=TRUE, quote="")
dim(tbl.gwas)
track <- GWASTrack("gwas 5k", tbl.gwas, chrom.col=12, pos.col=13, pval.col=28)
if(interactive()){
igv <- igvR()
setGenome(igv, "hg38")
setBrowserWindowTitle(igv, "GWAS demo")
displayTrack(igv, track)
Sys.sleep(1) # pause before zooming in
showGenomicRegion(igv, "chr6:32,240,829-32,929,353")
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.