DMRcate design
1
0
Entering edit mode
Jitendra ▴ 10
@nabiyogesh-11718
Last seen 4 months ago
United Kingdom

Hi all,

I am trying to use DMRcate first time for epic array data.

I have done the limma association analysis using below model design but not sure how to design model metrix for DMRcate for same type of BMI and CpG sites association analysis; Could you help me if below DMRcate design correct or wrong?

Limma linear regression;

#model matrix

    var<-model.matrix(~BMI + as.factor(Gender) + Age +CD8T +CD4T +NK + Bcell +Mono ,data=targets2)

    fit<-lmFit(mval,var)


    fit2<-eBayes(fit,trend=TRUE, robust=TRUE)

    probe<-topTable(fit2,adjust="BH",coef=2,num=Inf)

DMRcate analysis;

library(DMRcate)
myMs <- logit2(myBetas)

myMs.noSNPs <- rmSNPandCH(myMs, dist=2, mafcut=0.05)

design <- model.matrix(~BMI + as.factor(Gender) + Age +CD8T +CD4T +NK + Bcell +Mono ,data=targets2)

myannotation <- cpg.annotate(myMs.noSNPs, analysis.type="differential",design=design, coef=2)

Many thanks,

Epigenetics EpigeneticsWorkflow DMRcate • 1.3k views
ADD COMMENT
2
Entering edit mode
@james-w-macdonald-5106
Last seen 12 hours ago
United States

The help page for cpg.annotate says this:

design: Study design matrix. Identical context to differential
          analysis pipeline in 'limma'. Must have an intercept if
          'contrasts=FALSE'.  Applies only when 'analysis.type %in%
          c("differential", "ANOVA", "diffVar")'.
ADD COMMENT
0
Entering edit mode

Thanks James, but I did not understand it clearly. Could I used same Limma design matrix for DMRcate?

ADD REPLY
1
Entering edit mode

That's what identical context means.

ADD REPLY
0
Entering edit mode

Thanks, do I also need to mention 'contrasts=FALSE' in cpg.annotate.And also what to mentions in analysis.type="differential",/ "ANOVA", /"diffVar"" as we are looking for BMI associated CpG regions methylation patterns. Could anyone please help with this.

ADD REPLY

Login before adding your answer.

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