###################################################
### chunk number 1: 
###################################################
#line 40 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
#library(Biobase)
library(GlobalAncova)
library(globaltest)
sI <- sessionInfo()


###################################################
### chunk number 2: initialize
###################################################
#line 143 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
library(GlobalAncova)
library(globaltest)
library(golubEsets)
library(hu6800.db)
library(vsn)

data(Golub_Merge)
golubX <- justvsn(Golub_Merge)


###################################################
### chunk number 3: 
###################################################
#line 185 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
options(width=70)


###################################################
### chunk number 4: all
###################################################
#line 189 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gr <- as.numeric(golubX$ALL.AML=="ALL")
ga.all <- GlobalAncova(xx=exprs(golubX), group=gr, covars=NULL, perm=100)


###################################################
### chunk number 5: all2 eval=FALSE
###################################################
## #line 221 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
## GlobalAncova(xx=exprs(golubX), formula.full=~ALL.AML, formula.red=~1, model.dat=pData(golubX), perm=100)
## GlobalAncova(xx=exprs(golubX), formula.full=~ALL.AML, test.terms="ALL.AMLAML", model.dat=pData(golubX), perm=100)


###################################################
### chunk number 6: all.display
###################################################
#line 248 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
ga.all 


###################################################
### chunk number 7: gt.all
###################################################
#line 260 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.all <- gt(ALL.AML, golubX)


###################################################
### chunk number 8: 
###################################################
#line 263 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.all


###################################################
### chunk number 9: loadKEGG
###################################################
#line 276 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
kegg <- as.list(hu6800PATH2PROBE)


###################################################
### chunk number 10: 
###################################################
#line 280 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
cellcycle <- kegg[["04110"]]


###################################################
### chunk number 11: cc eval=FALSE
###################################################
## #line 291 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
## cellcycle <- kegg[["04110"]]


###################################################
### chunk number 12: cellcycle
###################################################
#line 297 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
ga.cc <- GlobalAncova(xx=exprs(golubX), group=gr, test.genes=cellcycle, method="both", perm=1000) 
ga.cc


###################################################
### chunk number 13: gt.cellcycle
###################################################
#line 304 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.cc <- gt(ALL.AML, golubX, subsets=cellcycle)
gt.cc


###################################################
### chunk number 14: 
###################################################
#line 308 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.cc


###################################################
### chunk number 15: adjust
###################################################
#line 321 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
ga.cc.BMPB <- GlobalAncova(xx=exprs(golubX), group=gr, covars=golubX$BM.PB, 
test.genes=cellcycle, method="both", perm=1000) 
ga.cc.BMPB


###################################################
### chunk number 16: gt.adjust
###################################################
#line 336 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.cc.BMPB <- gt(ALL.AML ~ BM.PB, golubX, subsets=cellcycle)
gt.cc.BMPB


###################################################
### chunk number 17: 
###################################################
#line 340 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.cc.BMPB


###################################################
### chunk number 18: 
###################################################
#line 351 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
data(vantVeer)
data(phenodata)	
data(pathways)	


###################################################
### chunk number 19: vV
###################################################
#line 376 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
data(vantVeer)
data(phenodata)	
data(pathways)
metastases <- phenodata$metastases
p53 <- pathways$p53_signalling	


###################################################
### chunk number 20: p53test
###################################################
#line 386 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
vV.1 <- GlobalAncova(xx=vantVeer, group=metastases, test.genes=p53, method="both", perm=1000)
vV.1


###################################################
### chunk number 21: grade
###################################################
#line 397 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
vV.3 <- GlobalAncova(xx=vantVeer, formula.full=~grade, formula.red=~1, model.dat=phenodata, test.genes=p53, method="both", perm=1000)
vV.3


###################################################
### chunk number 22: interact
###################################################
#line 413 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
signature.gene <- "AL137718"
model <- data.frame(phenodata, signature.gene=vantVeer[signature.gene, ])
vV.4 <- GlobalAncova(xx=vantVeer, formula.full=~signature.gene + ERstatus, formula.red=~ERstatus, model.dat=model, test.genes=p53, method="both", perm=1000)
vV.4


###################################################
### chunk number 23: coexpr
###################################################
#line 426 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
vV.5 <- GlobalAncova(xx=vantVeer, formula.full=~metastases + signature.gene + ERstatus, formula.red=~metastases + ERstatus, model.dat=model, test.genes=p53, method="both", perm=1000)
vV.5


###################################################
### chunk number 24: diffcoexpr
###################################################
#line 435 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
vV.6 <- GlobalAncova(xx=vantVeer, formula.full=~metastases * signature.gene + ERstatus, formula.red=~metastases + signature.gene + ERstatus, model.dat=model, test.genes=p53, method="both", perm=1000)
vV.6


###################################################
### chunk number 25: pathways
###################################################
#line 467 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
metastases <- phenodata$metastases
ga.pw <- GlobalAncova(xx=vantVeer, group=metastases, test.genes=pathways[1:4], method="both", perm=1000) 
ga.pw


###################################################
### chunk number 26: gt.pathways
###################################################
#line 477 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
gt.options(transpose = TRUE)
gt.pw <- gt(metastases, vantVeer, subsets=pathways[1:4])
gt.pw


###################################################
### chunk number 27: 
###################################################
#line 504 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
ga.pw.raw <- ga.pw[ ,"p.perm"] 
ga.pw.adj <- p.adjust(ga.pw.raw, "holm")
ga.result <- data.frame(rawp=ga.pw.raw, Holm=ga.pw.adj)
ga.result
gt.result <- p.adjust(gt.pw)
gt.result


###################################################
### chunk number 28: hypoth
###################################################
#line 534 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
if(require(Rgraphviz))
{
res <- GlobalAncova:::Hnull.family(pathways[1:4])
graph <- new("graphNEL", nodes=names(res), edgemode="directed")
graph <- addEdge(from=c(rep(names(res)[15],4),rep(names(res)[10],3),rep(names(res)[11],3),
  rep(names(res)[13],3),rep(names(res)[14],3),rep(names(res)[5],2),rep(names(res)[6],2),
  rep(names(res)[7],2),rep(names(res)[8],2),rep(names(res)[9],2),rep(names(res)[12],2)), 
  to=c(names(res)[10],names(res)[11],names(res)[13],names(res)[14],names(res)[5],names(res)[6],
  names(res)[8],names(res)[5],names(res)[7],names(res)[9],names(res)[6],names(res)[7],names(res)[12],
  names(res)[8],names(res)[9],names(res)[12],names(res)[1],names(res)[2],names(res)[1],names(res)[3],
  names(res)[1],names(res)[4],names(res)[2],names(res)[3],names(res)[2],names(res)[4],names(res)[3],
  names(res)[4]), graph, weights=rep(1, 28))

att <- list()
lab <- c("1","2","3","4","1-2","1-3","1-4","2-3","2-4","1-2-3","1-2-4","3-4","1-3-4","2-3-4","1-2-3-4")
names(lab) <- names(res)
att$label <- lab

plot(graph, nodeAttrs=att)
} else {
plot(1, 1, type="n", main="This plot requires Rgraphviz", 
     axes=FALSE, xlab="", ylab="")
}


###################################################
### chunk number 29: closed.test
###################################################
#line 562 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
ga.closed <- GlobalAncova.closed(xx=vantVeer, group=metastases, 
test.genes=pathways[1:4], previous.test=ga.pw, level=0.05, method="approx")


###################################################
### chunk number 30: names.closed.test
###################################################
#line 584 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
names(ga.closed)
rownames(ga.closed$test.results[[1]])
rownames(ga.closed$test.results[[1]]) <- NULL
ga.closed$test.results[1]
ga.closed$significant
ga.closed$not.significant


###################################################
### chunk number 31: gago
###################################################
#line 611 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
library(GO.db)
descendants <- get("GO:0007049", GOBPOFFSPRING)
gago <- GAGO(xx=exprs(golubX), formula.full=~ALL.AML, formula.red=~1, model.dat=pData(golubX),
      id=c("GO:0007049", descendants), annotation="hu6800", ontology="BP",
      multtest="focuslevel")


###################################################
### chunk number 32: 
###################################################
#line 618 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
head(gago)


###################################################
### chunk number 33: gakegg
###################################################
#line 640 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
GAKEGG(xx=exprs(golubX), formula.full=~ALL.AML, formula.red=~1, model.dat=pData(golubX),
      id=c("04110", "04210"), annotation="hu6800", multtest="BH")


###################################################
### chunk number 34: gabroad eval=FALSE
###################################################
## #line 659 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
## broad <- getBroadSets("your/path/to/msigdb_v.2.5.xml")
## GABroad(xx=exprs(golubX), formula.full=~ALL.AML, formula.red=~1, model.dat=pData(golubX),
##         id="chr5q33", collection=broad, annotation="hu6800")


###################################################
### chunk number 35: gabroad2 eval=FALSE
###################################################
## #line 667 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
## GABroad(xx=exprs(golubX), formula.full=~ALL.AML, formula.red=~1, model.dat=pData(golubX),
##         category=c("c1", "c3"), collection=broad, annotation="hu6800", multtest="holm")


###################################################
### chunk number 36: geneplot
###################################################
#line 719 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.genes(xx=vantVeer, group=metastases, test.genes=p53)
gt.vV <- gt(metastases, vantVeer, subsets=p53)
features(gt.vV, what="w")


###################################################
### chunk number 37: genes
###################################################
#line 727 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.genes(xx=vantVeer, group=metastases, test.genes=p53)


###################################################
### chunk number 38: gt_genes
###################################################
#line 741 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
features(gt.vV, what="w")


###################################################
### chunk number 39: geneplot
###################################################
#line 761 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.genes(xx=vantVeer, formula.full=~metastases, formula.red=~1, model.dat=phenodata, test.genes=p53, Colorgroup="ERstatus")


###################################################
### chunk number 40: genes2
###################################################
#line 767 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.genes(xx=vantVeer, formula.full=~metastases, formula.red=~1, model.dat=phenodata, test.genes=p53, Colorgroup="ERstatus")


###################################################
### chunk number 41: subjectsplot
###################################################
#line 803 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
#colnames(exprs(golubX)) <- pData(golubX)[ ,1]
Plot.subjects(xx=vantVeer, group=metastases, test.genes=p53, legendpos="bottomright")
subjects(gt.vV, what="w")


###################################################
### chunk number 42: subjects
###################################################
#line 811 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.subjects(xx=vantVeer, group=metastases, test.genes=p53, legendpos="bottomright")


###################################################
### chunk number 43: gt_subjects
###################################################
#line 824 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
subjects(gt.vV, what="w")


###################################################
### chunk number 44: subjectsplot2
###################################################
#line 844 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.subjects(xx=vantVeer, formula.full=~grade, formula.red=~1, model.dat=phenodata, test.genes=p53, Colorgroup="grade", legendpos="topleft")


###################################################
### chunk number 45: subjects2
###################################################
#line 850 "vignettes/GlobalAncova/inst/doc/GlobalAncova.rnw"
Plot.subjects(xx=vantVeer, formula.full=~grade, formula.red=~1, model.dat=phenodata, test.genes=p53, Colorgroup="grade", legendpos="topleft")