cpg.annotate model requirements
1
0
Entering edit mode
kadm • 0
@kadm-21455
Last seen 2.0 years ago

I am using ChAMP to process EPIC methylation array data. While delving into the source code for champ.DMR for another reason, I noticed that it defined the linear model used in the DMRcate method as ~ pheno. Further delving revealed that champ.DMR calls DMRcate::cpg.annotate, which itself requires a model with an intercept (e.g. it is has this in the function stopifnot(colnames(design)[1] == "(Intercept)")).

Is there any way to circumvent this and input a model into DMRcate (and thus, into champ.DMR) that has no intercept, as in 0 + pheno? This is necessary as other methods (like champ.DMP) already use a no intercept model.

DMRcate ChAMP • 335 views
ADD COMMENT
0
Entering edit mode

You can specify contrasts=TRUE and pass a contrasts matrix (as you would in limma) to cpg.annotate; this will allow you to fit non-intercept models.

Tim

ADD REPLY
0
Entering edit mode
Yuan Tian ▴ 240
@yuan-tian-13904
Last seen 21 hours ago
United Kingdom

Hi, Kadm:

Sorry that ChAMP did not designed parameter to set intercept in DMRcate, thus I suggest you directly run DMRcate package with specific settings. The result should not be too hard to organize.

Best Yuan Tian

ADD COMMENT

Login before adding your answer.

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