[BioC] [MEDIPS] mergeFrames: `Error in rbind(output, cbind(chromosomes[i], new_start, new_stop)) `

[Matthias Lienhard]

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 )
>
>
> https://gist.github.com/lindenb/6297130d57146ec2c2c0
>
> 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