Mapping between the probe set ID and the corresponding gene symbol
7 次查看(过去 30 天)
显示 更早的评论
Hi, iam using GSE21947 dataset to identify differential expressed gene . It uses [HG-U133A] Affymetrix Human Genome U133A Array. How to create (HuGeneFL_GeneSymbol_Map).
I would like to create a mapping between the probe set ID and the corresponding gene symbol provided as Map object in the MAT file eg:HuGeneFL_GeneSymbol_Map.
0 个评论
采纳的回答
Luuk van Oosten
2019-5-22
Dear Swarnambiga A,
Hope this is not too late (and what you are looking for, I'm more into proteins...), but I'll give it a shot:
If I understand correctly you are using Affymetrix array data, and you want to map the probe set ID to the corresponding gene.
There exists a bit of example code on how to this task in the documentation regarding exploring microarray gene expression data. Under the subsection 'Annotating Up-Regulated Genes Using Gene Onthology', there is explained what you want to do for some human data. They provide the map between the gene symbols associated with the GO IDs in a MAT-file, but also describe how you can get the latest data for that particlar human set from the GO database with annotations as follows:
GO = geneont('live',true);
HGann = goannotread('gene_association.goa_human',...
'Aspect','F','Fields',{'DB_Object_Symbol','GOid'});
HGmap = containers.Map();
for i = 1:numel(HGann)
key = HGann(i).DB_Object_Symbol;
if isKey(HGmap,key)
HGmap(key) = [HGmap(key) HGann(i).GOid];
else
HGmap(key) = HGann(i).GOid;
end
end
If you tweak that bit of code for your specific case, I think you will end up with what you want. If not, please provide some more info on what you tried (and what you want to get).
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Data Import and Management 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!