# code to produce the heatmap figure in the supplementary paper
Y<-as.matrix(read.table("CaltechY.dat"))
diag(Y)<-1
heatmap(Y,scale="none",col=c(0,1))