Why the order of sample names in the beta matrix affect the analysis results of Champ.DMP? What should I do?
0
0
Entering edit mode
Jinyue • 0
@1e365b75
Last seen 5 weeks ago
United States

Hello. I am using Champ.DMP in Champ package to find differential methylation cpg sites in my recent work. I first sorted the pd files by the phenotype column (that is, Sample_Group 1, 2, 3), but then I found that the order of the columns in the beta matrix will affect the cpg sites finally found. When I matched the column names of the beta matrix identical with the row names of the PD files, I even got the 'Error: ChAMP.DMP Have not detected even one significant CpGs.'. Which should be the correct operation? What should I do? My code is below. Any help would be appreciated so much!


# only sort pd file by pheno, that is, Sample_ Group

newpd <- newpd[order(newpd$Sample_Group),] myDMP <- champ.DMP(beta = correctedbeta, pheno = newpd$Sample_Group, arraytype="EPIC", adjPVal = 0.01)

[<<<<< 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 character type.
Your pheno information contains following groups. >>
<1>:55 samples.
<2>:52 samples.
<3>:56 samples.
[The power of statistics analysis on groups contain very few samples may not strong.]
pheno contains 3 phenotypes
compare.group parameter is NULL, EACH PAIR of phenotypes will be added into Compare List.
1_to_2 compare group : 1, 2
1_to_3 compare group : 1, 3
2_to_3 compare group : 2, 3

[ Section 1:  Check Input Pheno Done ]

[ Section 2:  Find Differential Methylated CpGs Start ]

-----------------------------
Start to Compare : 1, 2
Contrast Matrix
Contrasts
Levels p2-p1
p1    -1
p2     1
You have found 130 significant MVPs with a BH adjusted P-value below 0.01.
Calculate DMP for 1 and 2 done.
-----------------------------
Start to Compare : 1, 3
Contrast Matrix
Contrasts
Levels p3-p1
p1    -1
p3     1
You have found 802 significant MVPs with a BH adjusted P-value below 0.01.
Calculate DMP for 1 and 3 done.
-----------------------------
Start to Compare : 2, 3
Contrast Matrix
Contrasts
Levels p3-p2
p2    -1
p3     1
You have found 1947 significant MVPs with a BH adjusted P-value below 0.01.
Calculate DMP for 2 and 3 done.

[ Section 2:  Find Numeric Vector Related CpGs Done ]

[ Section 3:  Match Annotation Start ]

[ Section 3:  Match Annotation Done ]

[<<<<<< ChAMP.DMP END >>>>>>]
[===========================]
[You may want to process DMP.GUI() or champ.GSEA() next.]

# sort pd file by pheno, and then make column of beta matrix identical with the row of pd file

newpd <- newpd[order(newpd$Sample_Group),] correctedbeta <- correctedbeta[,newpd$Sample_Name]
myDMP <- champ.DMP(beta = correctedbeta, pheno = newpd$Sample_Group, arraytype="EPIC", adjPVal = 0.01) [<<<<< 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 character type. Your pheno information contains following groups. >> <1>:55 samples. <2>:52 samples. <3>:56 samples. [The power of statistics analysis on groups contain very few samples may not strong.] pheno contains 3 phenotypes compare.group parameter is NULL, EACH PAIR of phenotypes will be added into Compare List. 1_to_2 compare group : 1, 2 1_to_3 compare group : 1, 3 2_to_3 compare group : 2, 3 [ Section 1: Check Input Pheno Done ] [ Section 2: Find Differential Methylated CpGs Start ] ----------------------------- Start to Compare : 1, 2 Contrast Matrix Contrasts Levels p2-p1 p1 -1 p2 1 You have found 0 significant MVPs with a BH adjusted P-value below 0.01. Calculate DMP for 1 and 2 done. ----------------------------- Start to Compare : 1, 3 Contrast Matrix Contrasts Levels p3-p1 p1 -1 p3 1 You have found 0 significant MVPs with a BH adjusted P-value below 0.01. Calculate DMP for 1 and 3 done. ----------------------------- Start to Compare : 2, 3 Contrast Matrix Contrasts Levels p3-p2 p2 -1 p3 1 You have found 0 significant MVPs with a BH adjusted P-value below 0.01. Calculate DMP for 2 and 3 done. [ Section 2: Find Numeric Vector Related CpGs Done ] Error in champ.DMP(beta = a, pheno = newpd$Sample_Group, arraytype = "EPIC",  :
ChAMP.DMP Have not detected even one significant CpGs. You may try other threshold.

ChAMP limma Bioconductor ChAMPdata • 175 views
0
Entering edit mode

I can't understand your last error, it is said Error in champ.DMP(beta = a, pheno = newpd$Sample_Group, arraytype = "EPIC",, why do you have beta=a in input matrix ? You may have some mistakes in your code ADD REPLY 0 Entering edit mode Thanks for your comment! Sorry, I made an mistake while pasting the error information here. The correct information should be "Error in champ.DMP(beta = correctedbeta, pheno = newpd$Sample_Group, arraytype = "EPIC", : ChAMP.DMP Have not detected even one significant CpGs. You may try other threshold.". I still met the problem that the result was inconsistent with the previous one when I changed the column order of the beta matrix.