Question: Creating a GenomicRatioSet with non minfi normalized data
0
2.8 years ago by
hrishi27n20
hrishi27n20 wrote:

I have methylation data that is pre-processed and normalized with an unpublished normalization method, I am hoping to run DMR and DMP analysis on this dataset using minfi. Currently, I have Beta values, M-values, P-values and pData(sample_info, age,sex etc), I am trying to create the GenomicRatioSet with the data that I have, the idea is that once I have the GenomicRatioSet class I could use minfi to run the further analysis. Does the GenomicRatioSet function also add the genomic co-ordinates? so far I have written the following: Is this the right way to do this? The execution of the following was stopped due to an error. Would really appreciate your help and guidance on this.

Code:

library(minfi)

ratio <-GenomicRatioSet(Beta = Beta, M=Mval, pData=pdata)

print(ratio)

Error

Error in rownames<-(*tmp*, value = c("X10261170275986964480", "X10333227870024892416",  :

invalid rownames length

Calls: GenomicRatioSet ... SummarizedExperiment -> .local -> rownames<- -> rownames<-

Execution halted