Hi guys,
I can't seem to find an answer to this issue. It has been flagged before but without a successful solution.
cwp <- CentWaveParam(peakwidth=c(5,20), ppm=15, snthresh=5, prefilter = c(0,0)) xod <- findChromPeaks(raw, param = cwp)
xod <- adjustRtime(xod, param = ObiwarpParam())
Sample number 6 used as center sample. Aligning SJ45.CDF against SJ50.CDF ... Error: BiocParallel errors element index: 1, 2, 3, 4, 5, 6, ... first error: Dimensions of profile matrices of files SJ50.CDF and SJ45.CDF do not match!
I also tried running adjustRtime on the raw data files, but without success.
Strangely enough, when I run the old XCMS pipeline with rector instead of adjustRtime it all works, but I'd really like to try the new functions!
Any insights?
R version 3.4.0 (2017-04-21) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows >= 8 x64 (build 9200)
Hi Johannes,
Thanks for the reply.
Here's the full sessionInfo.
Sadly I won't be able to share the data - not my data and 12GB.
I did try various binSizes, but that didn't work :-( Could it be somewhere in the initial upload steps? Maybe it has to be binned before being RT corrected?
Good news is from the version you are using (1.52.0) to the current version (2.99.10) there have been many changes, improvements and fixes. Unfortunately we couldn't merge all fixes into the Bioconductor 3.5 release version because of some nasty git/svn complications.
Now, Bioconductor release 3.6 is on its way (release next week) with the updated and (presumably fixed) xcms. In the meantime you could try to install the current Bioconductor developmental version (3.6); should be save now, since it will be released soon:
Thank you! Will give this a shot and report back.
Ooo lala it worked :-)
Thanks for your help! Oh and if you have any recommendations for settings adjusted for Q-TOF for these new functions, let me know. I'm using as much as I can from my settings for the older xcmsSets, but the new workflow is quite a bit different and there's not much info out there as to adjustments. For example, how do I extract the results? Diffreport obviously won't work with the new expression sets.
The settings and parameters are essentially the same as in the old interface - just a little revamped to fit todays Bioconductor (and reuse functionality and classes from other packages). Eventually have a look at the new xcms vignette in the updated package - I try to explain some additional stuff and diagnostic plots there too.