Cross validation Balanced Tree

The following code directly generate the plot for cross validation using balanced tree. All relevant files from simulations can be accessed at http://www.tonyjhwueng.info/ououcir/crossvalidationV3/balancedtree/.

If not load well, please try again.

rm(list=ls())
par(mfrow=c(2,2))
par(mar=c(2,2,1.5,0.5))
size<-array(c(64,128,256,512),c(4))
for(sizeIndex in 1: length(size)){
#  sizeIndex<-1
  print(sizeIndex)
  loaddata<-paste("http://www.tonyjhwueng.info/ououcir/crossvalidationV3/balancedtree/",size[sizeIndex],"cvfullmodel.RData",sep = "")
  print(loaddata)
  load(url(loaddata))
  sdata<-as.matrix(s$conf.matrix$tol0.01)
  barplot(t(sdata),col=c("pink","skyblue","orange","green"),main=paste("taxa size ", size[sizeIndex]),ylim=c(0,100) ,xlab="", ylab="",beside=TRUE)
  #print(s$conf.matrix$tol0.01)
  }

## [1] 1
## [1] "http://www.tonyjhwueng.info/ououcir/crossvalidationV3/balancedtree/64cvfullmodel.RData"

## [1] 2
## [1] "http://www.tonyjhwueng.info/ououcir/crossvalidationV3/balancedtree/128cvfullmodel.RData"

## [1] 3
## [1] "http://www.tonyjhwueng.info/ououcir/crossvalidationV3/balancedtree/256cvfullmodel.RData"

## [1] 4
## [1] "http://www.tonyjhwueng.info/ououcir/crossvalidationV3/balancedtree/512cvfullmodel.RData"