Skip to content

Commit

Permalink
new EPIC manifest v4
Browse files Browse the repository at this point in the history
  • Loading branch information
perishky committed Sep 28, 2017
1 parent db55f65 commit e1b8526
Show file tree
Hide file tree
Showing 3 changed files with 14 additions and 10 deletions.
18 changes: 10 additions & 8 deletions data-raw/load-epic-manifest.r
Original file line number Diff line number Diff line change
@@ -1,7 +1,8 @@
load.epic.manifest <- function() {
## filename <- "ftp://webdata2:[email protected]/downloads/productfiles/methylationEPIC/infinium-methylationepic-v1-0-b1-manifest-file-csv.zip"
## filename <- "ftp://webdata2:[email protected]/downloads/productfiles/methylationEPIC/infinium-methylationepic-v1-0-b2-manifest-file-csv.zip"
filename <- "ftp://webdata2:[email protected]/downloads/productfiles/methylationEPIC/infinium-methylationepic-v1-0-b3-manifest-file-csv.zip"
## filename <- "ftp://webdata2:[email protected]/downloads/productfiles/methylationEPIC/infinium-methylationepic-v1-0-b3-manifest-file-csv.zip"
filename <- "ftp://webdata2:[email protected]/downloads/productfiles/methylationEPIC/infinium-methylationepic-v-1-0-b4-manifest-file-csv.zip"

on.exit(unlink(basename(filename)))
cat("Downloading", filename, "\n")
Expand All @@ -21,14 +22,15 @@ load.epic.manifest <- function() {
manifest$AddressB_ID <- sub("^0*", "", manifest$AddressB_ID)

## some addresses in the manifest don't appear in idat files
excluded.addresses <- readLines("epic-excludes.txt")
idx <- which(manifest$AddressA_ID %in% excluded.addresses
| manifest$AddressB_ID %in% excluded.addresses
| manifest$IlmnID %in% excluded.addresses
| manifest$IlmnID %in% excluded.addresses
| as.character(manifest$IlmnID) == "[Controls]")
## excluded.addresses <- readLines("epic-excludes.txt")
## idx <- which(manifest$AddressA_ID %in% excluded.addresses
## | manifest$AddressB_ID %in% excluded.addresses
## | manifest$IlmnID %in% excluded.addresses
## | manifest$IlmnID %in% excluded.addresses
## | as.character(manifest$IlmnID) == "[Controls]")
idx <- which(as.character(manifest$IlmnID) == "[Controls]")
if (length(idx) > 0)
manifest <- manifest[-idx,]
manifest <- manifest[-idx,]

manifest
}
Binary file modified inst/featuresets.rda
Binary file not shown.
6 changes: 4 additions & 2 deletions man/meffil.ewas.Rd

Some generated files are not rendered by default. Learn more about how customized files appear on GitHub.

0 comments on commit e1b8526

Please sign in to comment.