Question: [MEDIPS] mergeFrames: `Error in rbind(output, cbind(chromosomes[i], new_start, new_stop)) `
0
gravatar for Pierre Lindenbaum
5.2 years ago by
Pierre Lindenbaum10 wrote:
Hello, I'm trying to run the MEDIPS package with the following script (loop over each chromosomes, compare each conditions ) everything works fine for the chr1 but at some point, the script raises an error for the chr2. (see below). Does anyone know what I've done wrong here ? Thank you Pierre (...) Preprocessing MEDIPS SET 6 in MSet2... Calculating calibration curve... Performing linear regression... Intercept: 1.23986548019807 Slope: 1.72956076992792 Calculating relative methylation score... Estimate methylation... Differential coverage analysis... Extracting count windows with at least 1 reads... Extracting non-zero coupling factor windows... Execute edgeR for count data of 1524540 windows... (Neglecting parameter 'type') Creating a DGEList object... Apply trimmed mean of M-values (TMM) for library sizes normalization... Estimating common dispersion... Estimating tagwise dispersion... Calculating differential coverage... Adjusting p.values for multiple testing... Please note, log2 ratios are reported as log2(MSet1/MSet2). Creating results table... Adding differential coverage results... Total number of windows: 5701362 Number of windows tested for differential methylation: 1524540 Remaining number of windows with adjusted p.value<=0.1: 0 Error in rbind(output, cbind(chromosomes[i], new_start, new_stop)) : number of columns of matrices must match (see arg 2) Calls: myMedipAnalysis -> MEDIPS.mergeFrames -> rbind Execution halted
ADD COMMENTlink modified 5.2 years ago by Matthias Lienhard140 • written 5.2 years ago by Pierre Lindenbaum10
Answer: [MEDIPS] mergeFrames: `Error in rbind(output, cbind(chromosomes[i], new_start, n
0
gravatar for Matthias Lienhard
5.2 years ago by
Max Planck Institute for molecular Genetics, Berlin, Germany
Matthias Lienhard140 wrote:
Hi Pierre, as you are looping over the chromosomes and your code works for chr1 and fails for chr2, my first guess would be that there are no regions for chr2. Can you check for that? Best, Matthias Am 15.01.2014 05:13, schrieb Pierre Lindenbaum: > Hello, I'm trying to run the MEDIPS package with the following script > (loop over each chromosomes, compare each conditions ) > > > > > everything works fine for the chr1 but at some point, the script > raises an error for the chr2. (see below). > > Does anyone know what I've done wrong here ? > > Thank you > > Pierre > > > (...) > Preprocessing MEDIPS SET 6 in MSet2... > Calculating calibration curve... > Performing linear regression... > Intercept: 1.23986548019807 > Slope: 1.72956076992792 > Calculating relative methylation score... > Estimate methylation... > Differential coverage analysis... > Extracting count windows with at least 1 reads... > Extracting non-zero coupling factor windows... > Execute edgeR for count data of 1524540 windows... > (Neglecting parameter 'type') > Creating a DGEList object... > Apply trimmed mean of M-values (TMM) for library sizes normalization... > Estimating common dispersion... > Estimating tagwise dispersion... > Calculating differential coverage... > Adjusting p.values for multiple testing... > Please note, log2 ratios are reported as log2(MSet1/MSet2). > Creating results table... > Adding differential coverage results... > Total number of windows: 5701362 > Number of windows tested for differential methylation: 1524540 > Remaining number of windows with adjusted p.value<=0.1: 0 > Error in rbind(output, cbind(chromosomes[i], new_start, new_stop)) : > number of columns of matrices must match (see arg 2) > Calls: myMedipAnalysis -> MEDIPS.mergeFrames -> rbind > Execution halted > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD COMMENTlink written 5.2 years ago by Matthias Lienhard140
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: 228 users visited in the last hour