How to import data from MAIRplot to workspace using code?
3 views (last 30 days)
Show older comments
Mohammad Naser
on 7 Aug 2023
Commented: Mohammad Naser
on 25 Aug 2023
I'm using the MAIRplot in Matlab to differentiate between up-regulated and down-regulated genes. I can export these populations manually using the Export button. I want to run it for multiple data files and store the list in the workspace using code. Any suggestion?
Accepted Answer
Kausthub
on 23 Aug 2023
Hi Mohammad Naser,
I understand that you are looking for a way to store the up-regulated and down-regulated genes in the workspace using code rather than using the “Export” button since you are dealing with multiple data files.
A possible solution to the issue you are facing would be to utilize the outputs of the mairplot() function.
[Intensity, Ratio, H] = mairplot(DataX, DataY)
mairplot() has Ratio as one of the return values from which we can calculate the up-regulated and down regulated genes. Up-regulated genes have ratio greater than threshold and down-regulated genes has ratio lesser than negative of the threshold.
% Calculate threshold from Fold Change
threshold = log2(FoldChange); % threshold = 1; when foldChange = 2
% Find the respective indices
upIdx = ratio > threshold; % upIdx = ratio > 1; when foldChange = 2
downIdx = ratio < -threshold; % downIdx = ratio < -1; when foldChange = 2
diffIdx = upIdx | downIdx;
upIndices = find(upIdx);
downIndices = find(downIdx);
diffIndices = find(diffIdx);
% Find the respective labels
upLabels = maStruct.Names(upIdx);
downLabels = maStruct.Names(downIdx);
diffLabels = maStruct.Names(diffIdx);
Please refer to the following MATLAB documentation to know more about “mairplot” function
I hope this helps!
0 Comments
More Answers (0)
See Also
Categories
Find more on Bioinformatics Toolbox 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!