I am analyzing EPIC array data using the Champ program in Bioconductor. I have successfully generated DMRs between my 4 groups and an having trouble using the DMR.GUI function to view the DRMs. It looks like the program is unsuccessfully generating the contrast matrix using all 4 groups instead of the 2 groups I specified. My code is below. Any help would be appreciated!
DMR.GUI(myDMR, beta = myNorm, pheno = sampleRightOrder$Sample_Group, runDMP = TRUE, compare.group=c("1ECLP", "3ENC"), arraytype = "EPIC") !!! important !!! Since we just upgrated champ.DMP() function, which is now can support multiple phenotypes. Here in DMR.GUI() function, if you want to use "runDMP" parameter, and your pheno contains more than two groups of phenotypes, you MUST specify compare.group parameter as compare.group=c("A","B") to get DMP value between group A and group B.
[ Section 1: Calculate DMP Start ]
You pheno is factor type. Your pheno information contains following groups. >> <1ECLP>:7 samples. <3ENC>:8 samples. <4NEC>:8 samples. <2NECLP>:8 samples. Your pheno contains more than 2 phenotypes, please use compare.group to specify only two of them here. Your compare.group is in accord with your pheno, which is good, now we are about to extract information for your compare.group. Calculating DMP [===========================]
[<<<<< ChAMP.DMP START >>>>>]
!!! Important !!! New Modification has been made on champ.DMP():
(1): In this version champ.DMP() if your pheno parameter contains more than two groups of phenotypes, champ.DMP() would do pairewise differential methylated analysis between each pair of them. But you can also specify compare.group to only do comparasion between any two of them. (2): champ.DMP() now support numeric as pheno, and will do linear regression on them. So covariates like age could be inputted in this function. You need to make sure your inputted "pheno" parameter is "numeric" type.
[ Section 1: Check Input Pheno Start ]
You pheno is factor type. Your pheno information contains following groups. >> <1ECLP>:7 samples. <3ENC>:8 samples. [The power of statistics analysis on groups contain very few samples may not strong.] pheno contains only 2 phenotypes Your compare.group parameter is in accord with your pheno. Two pheno types has been added into Compare List. 1ECLPto3ENC compare group : 1ECLP, 3ENC
[ Section 1: Check Input Pheno Done ]
[ Section 2: Find Differential Methylated CpGs Start ]
Start to Compare : 1ECLP, 3ENC Contrast Matrix Contrasts Levels p2NECLP-p1ECLP p1ECLP -1 p2NECLP 1 p3ENC 0 p4NEC 0 Coefficients not estimable: p2NECLP p4NEC Error in contrasts.fit(fit, contrast.matrix) : trying to take contrast of non-estimable coefficient In addition: Warning message: Partial NA coefficients for 717600 probe(s)