How to import data from MAIRplot to workspace using code?
显示 更早的评论
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?
采纳的回答
Kausthub
2023-8-23
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 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
