plotBicluster {fabia} | R Documentation |
plotBicluster
: R implementation of plotBicluster
.
plotBicluster(r,p,opp=FALSE,zlim=NULL,title=NULL,which=c(1, 2))
r |
the result of extract_bic. |
p |
the bicluster to plot. |
opp |
plot opposite bicluster, default = FALSE. |
zlim |
vector containing a low and high value to use for the color scale. |
title |
title of the plot. |
which |
which plots are shown: 1=data matrix with bicluster on upper left, 2=plot of the bicluster; default c(1, 2). |
One bicluster is visualized by two plots. The variable "which" indicates which plots should be shown.
Plot1 (which=1): The data matrix is sorted such that the bicluster appear at the upper left corner. The bicluster is marked by a rectangle.
Plot2 (which=2): Only the bicluster is plotted.
Implementation in R.
Sepp Hochreiter
fabia
,
fabias
,
fabiap
,
fabi
,
fabiasp
,
mfsc
,
nmfdiv
,
nmfeu
,
nmfsc
,
extractPlot
,
extractBic
,
plotBicluster
,
Factorization
,
projFuncPos
,
projFunc
,
estimateMode
,
makeFabiaData
,
makeFabiaDataBlocks
,
makeFabiaDataPos
,
makeFabiaDataBlocksPos
,
matrixImagePlot
,
fabiaDemo
,
fabiaVersion
#--------------- # TEST #--------------- dat <- makeFabiaDataBlocks(n = 100,l= 50,p = 3,f1 = 5,f2 = 5, of1 = 5,of2 = 10,sd_noise = 3.0,sd_z_noise = 0.2,mean_z = 2.0, sd_z = 1.0,sd_l_noise = 0.2,mean_l = 3.0,sd_l = 1.0) X <- dat[[1]] Y <- dat[[2]] resEx <- fabia(X,3,0.01,20) rEx <- extractBic(resEx) plotBicluster(rEx,p=1) ## Not run: #--------------- # DEMO1 #--------------- dat <- makeFabiaDataBlocks(n = 1000,l= 100,p = 10,f1 = 5,f2 = 5, of1 = 5,of2 = 10,sd_noise = 3.0,sd_z_noise = 0.2,mean_z = 2.0, sd_z = 1.0,sd_l_noise = 0.2,mean_l = 3.0,sd_l = 1.0) X <- dat[[1]] Y <- dat[[2]] resToy <- fabia(X,13,0.01,200) rToy <- extractBic(resToy) plotBicluster(rToy,p=1) #--------------- # DEMO2 #--------------- avail <- require(fabiaData) if (!avail) { message("") message("") message("#####################################################") message("Package 'fabiaData' is not available: please install.") message("#####################################################") } else { data(Breast_A) X <- as.matrix(XBreast) resBreast <- fabia(X,5,0.1,200) rBreast <- extractBic(resBreast) plotBicluster(rBreast,p=1) } ## End(Not run)