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,

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")'.