Converting to GDS format
1
0
Entering edit mode
6.0 years ago
ebrown1955 ▴ 320

I'm trying to create nice looking LRR and BAF plots for some SNP array data that I have. I was hoping to achieve plots that also include a chromosomal ideogram. I was able to find "GWASTools" in R, however this requires my data to be in GDF format (genome data format). I have not been able to find how to create one of these files, or convert an illumina report or similar to it.

Can anyone share any expertise? Perhaps there are other tools in R that can work with a more palatable format that can also provide an ideogram plot.

Any help much appreciated!

genome array lrr baf snpintensity • 1.9k views
ADD COMMENT
0
Entering edit mode
6.0 years ago
ebrown1955 ▴ 320

Looks like I figured it out! From the GWASTools manual:

gdsfile <- tempfile()
path <- system.file("extdata", "illumina_raw_data", package="GWASdata")
data(illumina_snp_annot, illumina_scan_annot)
snpAnnot <- illumina_snp_annot[,c("snpID", "rsID", "chromosome","position", "alleleA", "alleleB")]
names(snpAnnot)[2] <- "snpName"
# subset of samples for testing
scanAnnot <- illumina_scan_annot[1:3, c("scanID", "genoRunID", "file")]
names(scanAnnot)[2] <- "scanName"
col.nums <- as.integer(c(1,2,12,13))
names(col.nums) <- c("snp", "sample", "a1", "a2")
diagfile <- tempfile()
ADD COMMENT

Login before adding your answer.

Traffic: 1989 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6