How does one subtract two msheatmaps to produce a result that highlights differences between the two heatmaps

1 view (last 30 days)
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))

Answers (0)

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!