Question: Creating a GenomicRatioSet with non minfi normalized data
0
gravatar for hrishi27n
2.9 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)

Beta <-as.matrix(read.csv("testGenomic.csv", header=T,sep=",", row.names=1))

Mval <- as.matrix(read.csv("testGenomicMVal.csv", header=T,sep=",", row.names=1))

pdata <-read.table("MaleTargets.csv",header=TRUE,sep=",")

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

 

ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by hrishi27n20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 186 users visited in the last hour