Valid to normalise MethylumiSet using minfi?
1
0
Entering edit mode
Pem • 0
@pem-20182
Last seen 5.1 years ago

Hi,

I am working on a dataset previously read in using methylumIDAT. As this results in a MethylumiSet object, I am wondering whether it is valid to convert this to an RGChannelSet using methylumi:::methylumiToMinfi in order to normalise using the preprocessXXX functions available in the minfi library?

Thanks

normalization minfi methylumi • 716 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 12 hours ago
United States

You seem to be asking two questions simultaneously.

  1. Should I use an un-exported function to do something?
  2. Is it valid to do something that a package specifically has functions designed to accomplish?

The answer to #1 is usually no. Unexported functions are not exported for a reason. The conventional way to convert would be to follow the instructions in the vignette, in section 4, an example of which would be (using the data used in that vignette):

> library(methylumi)
> library(FDb.InfiniumMethylation.hg19)
> library(TCGAMethylation450k)
> idatPath <- system.file("extdata/idat", package = "TCGAMethylation450k")
> mset450k <- methylumIDAT(getBarcodes(path = idatPath), idatPath = idatPath)
0 HumanMethylation27 samples found
10 HumanMethylation450 samples found
Attempting to extract protocolData() from list...
Determining chip type from IDAT protocolData...
> rgset <- as(mset450k, "RGChannelSet")
Loading required package: rtracklayer
Fetching coordinates for hg19...
> rgset
class: RGChannelSet 
dim: 621928 10 
metadata(0):
assays(2): Green Red
rownames(621928): 14782418 12709357 ... 28673402 13742412
rowData names(0):
colnames(10): 6005486014_R01C01 6005486014_R02C02 ... 6005486023_R04C02
  6005486023_R05C01
colData names(1): barcode
Annotation
  array: IlluminaHumanMethylation450k

As for question #2, I think the answer has to be sure, why not?

ADD COMMENT

Login before adding your answer.

Traffic: 542 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