###################################################
### chunk number 1: loadLib
###################################################
#line 57 "vignettes/bioDist/inst/doc/bioDist.Rnw"
library("bioDist")
data(sample.ExpressionSet)
exData = t(exprs(sample.ExpressionSet))


###################################################
### chunk number 2: R
###################################################
#line 89 "vignettes/bioDist/inst/doc/bioDist.Rnw"
s1 = MIdist(exData)
s2 = as.matrix(s1)
dim(s2)
r1 = mutualInfo(exData)


###################################################
### chunk number 3: KLD
###################################################
#line 100 "vignettes/bioDist/inst/doc/bioDist.Rnw"
kl1 = KLdist.matrix(exData)
kl2 = KLD.matrix(exData, method="density", supp=range(exData))


###################################################
### chunk number 4: Tau
###################################################
#line 112 "vignettes/bioDist/inst/doc/bioDist.Rnw"
eS = sample.ExpressionSet[401:500,]
tauD = tau.dist(eS, sample=FALSE)
sp = spearman.dist(eS, sample=FALSE)


###################################################
### chunk number 5: closest
###################################################
#line 121 "vignettes/bioDist/inst/doc/bioDist.Rnw"
f1 = featureNames(eS)[1]
closest.top(f1, sp, 3)


###################################################
### chunk number 6: 
###################################################
#line 131 "vignettes/bioDist/inst/doc/bioDist.Rnw"
sessionInfo()