Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Topic: Looking for help: Dimensions of profile matrices do not match ! #400

Open
wants to merge 537 commits into
base: XChromatograms
Choose a base branch
from

Conversation

genefisher
Copy link

Hey ,

I met some problem in using XC MS, I can not fix by myself now. I am looking for help in here.

I divided my samples(650 samples) into two sub groups named mock and treatment and then put them into one folder to run XCMS.

When I run mock(325 samples) independently, the XC MS runs well and it works fine for me.

When I run treatment(325 samples) independently, the the XC MS also runs well.

However, when I run mock and treatment samples(650 samples) together, the XC MS shows error as follows:

Control0h_L1-50_BD7_01_38715 Create profile matrix with method 'bin' and step 0.1 ... OK
Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-51_BE7_01_38716 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-52_BF7_01_38717 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-53_BG7_01_38718 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-54_BH7_01_38719 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-55_BA8_01_38720 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-56_BB8_01_38721 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-57_BC8_01_38722 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-58_BD8_01_38723 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-59_BE8_01_38724 Create profile matrix with method 'bin' and step 0.1 ... OK
Control0h_L1-6_BG1_01_38666 Create profile matrix with method 'bin' and step 0.1 ... OK
Found gaps: cut scantime-vector at 1399.01 seconds
Control0h_L1-60_BF8_01_38725 Create profile matrix with method 'bin' and step 0.1 ... OK

Error in .local(object, ...) :

Dimensions of profile matrices do not match !

xset3 <- group(xset2, method="density", mzwid=0.01, sleep=0.001,

  •            minfrac=0.1, minsamp=1, bw=5)
    

Error in group(xset2, method = "density", mzwid = 0.01, sleep = 0.001, :
object 'xset2' not found

xset4 <- fillPeaks(xset3)
Error in fillPeaks(xset3) : object 'xset3' not found
an <- xsAnnotate(xset4)
Error in xsAnnotate(xset4) : object 'xset4' not found
anF <- groupFWHM(an, perfwhm = 0.6)
Error in groupFWHM(an, perfwhm = 0.6) : object 'an' not found
anI <- findIsotopes(anF, mzabs=0.01)
Error in findIsotopes(anF, mzabs = 0.01) : object 'anF' not found
anIC <- groupCorr(anI, cor_eic_th=0.75)
Error in groupCorr(anI, cor_eic_th = 0.75) : object 'anI' not found
anFA <- findAdducts(anIC, polarity="positive")
Error in findAdducts(anIC, polarity = "positive") :
object 'anIC' not found
write.csv(getPeaklist(anFA), file="data.csv")
Error in getPeaklist(anFA) : object 'anFA' not found

Many thanks for you in advance!

@sneumann
Copy link
Owner

Hi, I think something went wrong here, can you close
and open as issue rather than pull request ?
Also, please include the output of sessionInfo()
Yours, Steffen

@codecov-io
Copy link

codecov-io commented Sep 12, 2019

Codecov Report

Merging #400 into XChromatograms will not change coverage.
The diff coverage is n/a.

Impacted file tree graph

@@               Coverage Diff               @@
##           XChromatograms     #400   +/-   ##
===============================================
  Coverage           49.17%   49.17%           
===============================================
  Files                  79       79           
  Lines               19731    19731           
===============================================
  Hits                 9702     9702           
  Misses              10029    10029           

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update 5cd366a...5cd366a. Read the comment docs.

@github-pages github-pages bot temporarily deployed to github-pages April 30, 2020 11:02 Inactive
jorainer and others added 27 commits May 11, 2020 10:37
- Export `Param` class.
- Don't drop initial feature definitions in the `featureDefinitions<-` method
  allowing to add additional columns to this data frame.
- The `filterFile` method gains parameter `keepFeatures` to force keeping
  correspondence result even due filtering by samples (files).
…d code not new expected behavior and type should be 'double' rather than 'numeric'
- `featureChromatograms` gets parameters `n` and `value` to extract EICs from
  only the top n samples with highest intensities.
- Add related unit tests and documentation.
- Add `normalize,Chromatogram` and `normalize,Chromatograms` to normalize
  intensities in a chromatogram object.
- Add related documentation and unit tests.
PR to resolve issue #464 by fixing group matrix type
- Fix issue in conversion from `XCMSnExp` to `xcmsSet` ensuring that a phenodata
  `data.frame` will **always** be set. Thanks to Andris Jankevics for reporting.
Several updates, fixes etc
- Improve performance of splitting of objects for paralle processing.
- Should improve the speed of `findChromPeaks`, `refineChromPeaks`,
  `chromPeakSpectra` and `reconstructChromPeakSpectra`.
- Further improve splitting an `XCMSnExp` object by file by first splitting the
  chrom peak data.
- Use `extractROWS` instead of `[` for `DataFrame`.
- Small update to refineChromPeaks to use bplapply instead of bpmapply.
sneumann and others added 30 commits April 18, 2022 11:29
[DOC] Replace links to MassSpecWavelet::sav.gol() with links to peakDetectionCWT()
- Fix problem with `fillChromPeaks` after `MatchedFilterParam` peak detection on
  very sparse data (i.e. data with many empty spectra).
Fix problem in fillChromPeaks with sparse data
- Reimplement the `reconstructChromPeakSpectra` function for DIA data analysis
  to fix an issue with chrom peaks in overlapping isolation windows for SWATH
  data and to generally improve performance.
Improve performance for SWATH/DIA data analysis
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.