Creates a heatmap and returns a similarity matrix of the Jaccard Index (Row, Column or both dimensions) in order to compare 2 different biclustering results or compare the biclusters of a single result.
CompareResultJI(BCresult1, BCresult2 =NULL, type ="both", plot =TRUE)
Arguments
BCresult1: A S4 Biclust object. If only this input Biclust object is given, the biclusters of this single result will be compared.
BCresult2: A second S4 Biclust object to which BCresult1 should be compared. (default=NULL)
type: Of which dimension should the Jaccard Index be computed? Can be "row", "col" or "both" (default).
plot: Logical value if plot should be outputted (default=TRUE).
Returns
A list containing
SimMat: The JI Similarity Matrix between the compared biclusters.
MaxSim: A list containing the maximum values on each row (BCResult1) and each column (BCResult2).
Details
The Jaccard Index between two biclusters is calculated as following:
JI(BC1,BC2)=m12(m1+m2−m12)
in which
type="row" or type="col"
m1= Number of rows/columns of BC1
m2= Number of rows/columns of BC2
m12= Number of rows/columns of union of row/column membership of BC1 and BC2
type="both"
m1= Size of BC1 (rows times columns)
m2= Size of BC2 (rows times columns)
m12=m1+m2− size of overlapping BC of BC1 and BC2
Examples
## Not run:data <- matrix(sample(c(0,1),100*100,replace=TRUE,prob=c(0.9,0.1)),nrow=100,ncol=100)data[1:10,1:10]<-1# BC1data[11:20,11:20]<-1# BC2data[21:30,21:30]<-1# BC3data <- data[sample(1:nrow(data),nrow(data)),sample(1:ncol(data),ncol(data))]# Result 1result1 <- bibit(data,minr=5,minc=5)result1
# Result 2result2 <- bibit(data,minr=2,minc=2)result2
## Compare all BC's of Result 1 ##Sim1 <- CompareResultJI(BCresult1=result1,type="both")Sim1$SimMat
## Compare BC's of Result 1 and 2 ##Sim12 <- CompareResultJI(BCresult1=result1,BCresult2=result2,type="both",plot=FALSE)str(Sim12)## End(Not run)