Question: Error in Champ for Bioconductor analyzing Methlayion array data
0
gravatar for erinkello
10 days ago by
erinkello0
erinkello0 wrote:

Hello,

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)

ADD COMMENTlink modified 10 days ago by Yuan Tian40 • written 10 days ago by erinkello0
Answer: Error in Champ for Bioconductor analyzing Methlayion array data
0
gravatar for Yuan Tian
10 days ago by
Yuan Tian40
Shanghai Institute for Biology Science, Shanghai, China
Yuan Tian40 wrote:

Hi:

The error should be caused by NA value in your pheno parameter. Please try to remove them.

Best Yuan Tian

ADD COMMENTlink written 10 days ago by Yuan Tian40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 259 users visited in the last hour