###################################################
### chunk number 1: set width
###################################################
#line 29 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
options(width=60)
options(continue=" ")
try(X11.options(type="xlib"), silent=TRUE)


###################################################
### chunk number 2: packages
###################################################
#line 69 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
library(biclust)
library(eisa)
library(ALL)
data(ALL)


###################################################
### chunk number 3: filter
###################################################
#line 81 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
library(GO.db)
library(hgu95av2.db)
gotable <- toTable(GOTERM)
myterms <- unique(gotable$go_id[gotable$Term %in% c("immune system process")])
myprobes <- unique(unlist(mget(myterms, hgu95av2GO2ALLPROBES)))
ALL.filtered <- ALL[myprobes,]


###################################################
### chunk number 4: filtered
###################################################
#line 91 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
nrow(ALL.filtered)


###################################################
### chunk number 5: seed
###################################################
#line 96 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
set.seed(0xf00)


###################################################
### chunk number 6: bcplaid
###################################################
#line 102 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
Bc <- biclust(exprs(ALL.filtered), BCPlaid(),
              fit.model = ~m + a + b, verbose = FALSE)


###################################################
### chunk number 7: bcplaid result
###################################################
#line 108 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
class(Bc)
Bc


###################################################
### chunk number 8: convert bc to isa
###################################################
#line 122 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
Bc <- annotate(Bc, ALL.filtered)
modules <- as(Bc, "ISAModules")
modules


###################################################
### chunk number 9: biclust enrichment
###################################################
#line 135 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
library(KEGG.db)
KEGG <- ISAKEGG(modules)
sigCategories(KEGG)[[2]]
unlist(mget(sigCategories(KEGG)[[2]], KEGGPATHID2NAME))


###################################################
### chunk number 10: heatmap eval=FALSE
###################################################
## #line 147 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
## col <- ifelse(grepl("^B", ALL.filtered$BT), "white", "black")
## modcol <- col[ getSamples(modules, 2)[[1]] ]
## ISA2heatmap(modules, 2, ALL.filtered, 
##             ColSideColors=modcol)


###################################################
### chunk number 11: heatmap-real
###################################################
#line 156 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
#line 147 "vignettes/eisa/inst/doc/EISA_biclust.Rnw#from line#156#"
col <- ifelse(grepl("^B", ALL.filtered$BT), "white", "black")
modcol <- col[ getSamples(modules, 2)[[1]] ]
ISA2heatmap(modules, 2, ALL.filtered, 
            ColSideColors=modcol)
#line 157 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"


###################################################
### chunk number 12: profilePlot eval=FALSE
###################################################
## #line 174 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
## profilePlot(modules, 2, ALL, plot="both")


###################################################
### chunk number 13: profilePlot-real
###################################################
#line 179 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
#line 174 "vignettes/eisa/inst/doc/EISA_biclust.Rnw#from line#179#"
profilePlot(modules, 2, ALL, plot="both")
#line 180 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"


###################################################
### chunk number 14: psoptions
###################################################
#line 193 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
ps.options(fonts=c("serif", "mono"))


###################################################
### chunk number 15: GOtreeplot eval=FALSE
###################################################
## #line 196 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
## library(GO.db)
## GO <- ISAGO(modules)
## gog <- gograph(summary(GO$CC)[[2]])
## summary(gog)
## gographPlot(gog)


###################################################
### chunk number 16: GOtreeplot-real
###################################################
#line 205 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
#line 196 "vignettes/eisa/inst/doc/EISA_biclust.Rnw#from line#205#"
library(GO.db)
GO <- ISAGO(modules)
gog <- gograph(summary(GO$CC)[[2]])
summary(gog)
gographPlot(gog)
#line 206 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"


###################################################
### chunk number 17: html
###################################################
#line 218 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
CHR <- ISACHR(modules)
htmldir <- tempdir()
ISAHTML(eset=ALL.filtered, modules=modules, target.dir=htmldir, 
        GO=GO, KEGG=KEGG, CHR=CHR, condPlot=FALSE)


###################################################
### chunk number 18: html-show eval=FALSE
###################################################
## #line 224 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
## if (interactive()) {
##   browseURL(URLencode(paste("file://", htmldir, "/index.html", sep="")))
## }


###################################################
### chunk number 19: mnplot eval=FALSE
###################################################
## #line 236 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
## group <- ifelse(grepl("^B", ALL.filtered$BT), "B-cell", "T-cell")
## ISAmnplot(modules, 2, ALL.filtered, norm="raw", group=group)


###################################################
### chunk number 20: mnplot-real
###################################################
#line 242 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
#line 236 "vignettes/eisa/inst/doc/EISA_biclust.Rnw#from line#242#"
group <- ifelse(grepl("^B", ALL.filtered$BT), "B-cell", "T-cell")
ISAmnplot(modules, 2, ALL.filtered, norm="raw", group=group)
#line 243 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"


###################################################
### chunk number 21: load isa data
###################################################
#line 260 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
data(ALLModules)
ALLModules


###################################################
### chunk number 22: isa to biclust
###################################################
#line 267 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
BcMods <- as(ALLModules, "Biclust")
BcMods


###################################################
### chunk number 23: coherence
###################################################
#line 277 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
data <- exprs(ALL[featureNames(ALLModules),])
constantVariance(data, BcMods, 1)
additiveVariance(data, BcMods, 1)
multiplicativeVariance(data, BcMods, 1)
signVariance(data, BcMods, 1)


###################################################
### chunk number 24: coherence vs robustness
###################################################
#line 287 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
cV <- sapply(1:BcMods@Number, function(x) constantVariance(data, BcMods, x))
aV <- sapply(1:BcMods@Number, function(x) additiveVariance(data, BcMods, x))
mV <- sapply(1:BcMods@Number, function(x) multiplicativeVariance(data, BcMods, x))
sV <- sapply(1:BcMods@Number, function(x) signVariance(data, BcMods, x))
rob <- ISARobustness(ALL, ALLModules)


###################################################
### chunk number 25: pairs eval=FALSE
###################################################
## #line 297 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
## panel.low <- function(x, y) {
##   usr <- par("usr")
##   m <- c((usr[2]+usr[1])/2, (usr[4]+usr[3])/2)
##   text(m[1], m[2], adj=c(1/2,1/2), cex=1.5,
##        paste(sep="\n", "Correlation:", round(cor(x,y),2)))
## }
## pairs( cbind(cV, aV, mV, sV, rob), lower.panel=panel.low )


###################################################
### chunk number 26: pairs-real
###################################################
#line 308 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
#line 297 "vignettes/eisa/inst/doc/EISA_biclust.Rnw#from line#308#"
panel.low <- function(x, y) {
  usr <- par("usr")
  m <- c((usr[2]+usr[1])/2, (usr[4]+usr[3])/2)
  text(m[1], m[2], adj=c(1/2,1/2), cex=1.5,
       paste(sep="\n", "Correlation:", round(cor(x,y),2)))
}
pairs( cbind(cV, aV, mV, sV, rob), lower.panel=panel.low )
#line 309 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"


###################################################
### chunk number 27: sessioninfo
###################################################
#line 328 "vignettes/eisa/inst/doc/EISA_biclust.Rnw"
toLatex(sessionInfo())