Basic Examples (2)
Basic Examples
(2)
Gene expression (nTPM) in different human tissues:
In[1]:=
Dataset@ResourceData
Out[1]=
Expression of 50 strongest expressing genes in human tissues ( color indicates data not available):
In[2]:=
RNAExpressionData=ResourceData;tissues=RNAExpressionData[[1]];ensemblIDGeneAssoc=ResourceData,"Ensembl ID gene name association";ensemblIDs=Keys@RNAExpressionData[[2;;-1]];genes=ensemblIDGeneAssoc/@ensemblIDs;v=Values@RNAExpressionData;totExp=Table[Total@v[[i]]/.{"NA"->0},{i,2,Length@genes+1}];sv=Sort[totExp];numGenes=50;pos=Flatten@Map[Position[totExp,#]&,sv[[Length@genes-numGenes;;Length@genes]]];highestExpressionGenes=ensemblIDs[[pos]];highestExpressionGeneNames=Map[ensemblIDGeneAssoc,highestExpressionGenes];ftx=MapThread[{#1,#2}&,{Range@Length@highestExpressionGeneNames,highestExpressionGeneNames}];fty=MapThread[{#1,Rotate[Capitalize@#2,Pi/2]}&,{Range@Length@tissues,tissues}];array=v[[pos+1]];mA={Min@array/."NA"->0,Max@array/."NA"->0};ArrayPlot[array,ColorFunction->"BlueGreenYellow",FrameTicks->{ftx,fty},ImageSize->1000,Frame->True,FrameStyle->Directive[Black,12],PlotLegends->BarLegend[{"BlueGreenYellow",mA}],PlotLabel->Style["Highest expressing genes in human tissues (nTPM)",Black,Bold,20]]
Out[18]=
Visualizations (3)
Visualizations
(3)
Analysis (2)
Analysis
(2)