Question: Unpaired Data with DMRMark
0
gravatar for lizlitkowski
7 months ago by
lizlitkowski0 wrote:

I am trying to find differentially methylated regions using the R package DMRMark. My data is unpaired data between controls and cases and I went back to the example in the DMRMark Reference Manual to make sure I had the steps. In the code below, the first 5 lines are directly from the manual.

I was trying to use the output from the 'reformData' function as the input to the 'DMRMark' function but I tried it 3 ways and I get errors. I expect the first option to work but it gives the error that I have 'NA's in my data. The reformData function creates the 'NA's so that error was a particular surprise.

Any thoughts on what I am doing wrong here?

# The case with One more sample from Tumour group
# The second Tumour sample is the extra one
mv2 <- matrix(1:25,5)
mv2[,2] <- mv2[,2] + 5
patient <- factor(c(1,3,1:3))
type = c(rep("Normal",2),rep("Tumour",3))
pd <- model.matrix(~patient + type + 0)
exmpl = reformData(mv2, pd)
L = c(300,2578,608,362,48)
set.seed(13)
starting = 1 
par <- DMRMark(exmpl, L, starting)
par <- DMRMark(exmpl, L, starting, pd = pd)
par <- DMRMark(mv2, L, starting, pd = pd)
methylation • 229 views
ADD COMMENTlink written 7 months ago by lizlitkowski0
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 2.3.0
Traffic: 1495 users visited in the last hour