Clear Filters
Clear Filters

How to import data from MAIRplot to workspace using code?

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

More Answers (0)

Categories

Find more on Bioinformatics Toolbox in Help Center and File Exchange

Products


Release

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!