DMRcate design
0
0
Entering edit mode
7 months ago
Bioinfonext ▴ 460

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)

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.

Many thanks

R Epigenetics EPIC DMRcate • 258 views
ADD COMMENT

Login before adding your answer.

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