graph - Group variables by clusters on heatmap in R -

-

i trying reproduce first figure of this paper on graph clustering: enter image description here

here sample of adjacency matrix: 

data=cbind(c(48,0,0,0,0,1,3,0,1,0),c(0,75,0,0,3,2,1,0,0,1),c(0,0,34,1,16,0,3,0,1,1),c(0,0,1,58,0,1,3,1,0,0),c(0,3,16,0,181,6,6,0,2,2),c(1,2,0,1,6,56,2,1,0,1),c(3,1,3,3,6,2,129,0,0,1),c(0,0,0,1,0,1,0,13,0,1),c(1,0,1,0,2,0,0,0,70,0),c(0,1,1,0,2,1,1,1,0,85)) colnames(data)=letters[1:nrow(data)] rownames(data)=colnames(data)

and these commands obtain following heatmap:

library(reshape) library(ggplot2) data.m=melt(data) data.m[,"rescale"]=round(rescale(data.m[,"value"]),3) p=ggplot(data.m,aes(x1, x2))+geom_tile(aes(fill=rescale),colour="white")  p=p+scale_fill_gradient(low="white",high="black") p+theme(text=element_text(size=10),axis.text.x=element_text(angle=90,vjust=0))

enter image description here

this similar plot on left of figure 1 above. differences (1) nodes not ordered randomly alphabetically, , (2) instead of having binary black/white pixels, using "shades of grey" palette able show strength of co-occurrence between nodes.

but point is hard distinguish cluster structure (and more true full set of 100 nodes). so, want order vertices clusters on heatmap. have membership vector community detection algorithm:

membership=c(1,2,4,2,5,3,1,2,2,3)

now, how can obtain heatmap similar plot on right of figure 1 above?

thanks lot in advance help

ps: have experimented r draw kmeans clustering heatmap , r: how display clustered matrix heatmap (similar color patterns grouped) not want.

turned out extremely easy. still posting solution others in case don't waste time on did.

the first part same before:

data.m=melt(data) data.m[,"rescale"]=round(rescale(data.m[,"value"]),3)

now, trick levels of factors of melted data.frame have ordered membership:

data.m[,"x1"]=factor(data.m[,"x1"],levels=levels(data.m[,"x1"])[order(membership)]) data.m[,"x2"]=factor(data.m[,"x2"],levels=levels(data.m[,"x2"])[order(membership)])

then, plot heat map (same before): 

p=ggplot(data.m,aes(x1, x2))+geom_tile(aes(fill=rescale),colour="white")  p=p+scale_fill_gradient(low="white",high="black") p+theme(text=element_text(size=10),axis.text.x=element_text(angle=90,vjust=0))

enter image description here

this time, cluster visible.


Comments

Post a Comment

Popular posts from this blog

node.js - Mongoose: Cast to ObjectId failed for value on newly created object after setting the value -

-
this question has answer here: what's mongoose error cast objectid failed value xxx @ path “_id”? 7 answers even though has been asked many times because entity has no valid objectid field. i came across situation objectid set value can casted objectid still error. the schema: var serviceschema = parammongoose.schema({ servicename:{ type: parammongoose.schema.types.objectid , ref: 'servicenames',required:true } }); pseudo code: servicefromdb=new service({servicename:'some name'}); servicefromdb.servicename='000000000000000000000001'; servicefromdb.save(function(paramerror,paramdata){ if(paramerror){ console.log('but but...',servicefromdb,paramerror) } }); the output of code is: but but... { servicename: 000000000000000000000001, _id: 55079a90286f49280364f78b } { [casterror: cast objectid faile...
Read more

gradle error "Cannot convert the provided notation to a File or URI" -

-
i have defined 2 system properties in gradle.properties: systemprop.builddir=build systemprop.cfgdir=build\\cfg and have following task defined in build.gradle: task clean(group:'clean',description:'clean project') << { ant.sequential { delete(dir: system.properties.builddir) mkdir(dir: system.properties.builddir) delete(dir: system.properties.cfgdir) mkdir(dir: system.properties.cfgdir) } } this generates following error: execution failed task ':clean'. > cannot convert provided notation file or uri: {dir=build}. following types/formats supported: - string or charsequence path, e.g 'src/main/java' or '/usr/include' - string or charsequence uri, e.g 'file:/usr/include' - file instance. - uri or url instance. but equivalent block of code not generate error , works expected: task clean(group:'clean',description:'clean project') << { ...
Read more

python - NameError: name 'subprocess' is not defined -

-
i'm trying make rng website: http://jamesdotcom.com/?p=417 . when run following code: captureimage = subprocess.popen( [ "fswebcam", "-r", "356x292", "-d", "/dev/video0", "static.jpg", "--skip", "10" ], stdout=devnull, stderr=devnull) captureimage.communicate() i following error message: nameerror: name 'subprocess' not defined what should do? thanks have imported first? import subprocess captureimage = subprocess.popen(["fswebcam", ...) subprocess included in python standard library since python 2.4. should not necessary install it.
Read more