How does one subtract two msheatmaps to produce a result that highlights differences between the two heatmaps
1 view (last 30 days)
Show older comments
Is there a way to subtract two heatmaps to produce a third that shows only the differences between the two ?. The problem lies in the fact that mzXML2peaks function stores the ion intensities as a cell array and its impossible to subtract two different cell arrays in MATLAB. Even after resampling, the problem still persists. Tvec is just a brute force alignment solution that excludes the need for msalign functions. How do I subtract the two heat maps ?
sample1= mzxmlread('Sample1.mzXML') [ps,ts] = mzxml2peaks(sample1,'level',1); [MZs,Ys] = msppresample(ps,5000);
sample2= mzxmlread('Sample2.mzXML') [pg,tg] = mzxml2peaks(sample2,'level',1); [MZg,Yg] = msppresample(pg,5000);
plot( msheatmap(MZs,tvec,log(Ys2),'resolution',0.15)- msheatmap(MZg,tvec,log(Yg2),'resolution',0.15))
0 Comments
Answers (0)
See Also
Categories
Find more on Time Series in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!