Summary Analysis

rm(list=ls())
setwd("~/Dropbox/NSCT/112年計畫113年八月/phygarch11rate/rcode/simulationunif/")

treetype.array <- c("balanced","left","coalescent","bdtree") #c("bdtree")
taxasize.array <- c(16, 32, 64, 128) #c(22)
for(taxasizeIndex in 1:length(taxasize.array)){
  #taxasizeIndex<-4
  taxasize<- taxasize.array[taxasizeIndex]
  for(treetypeIndex in 1:length(treetype.array)){
    #treetypeIndex<-1
    treetype<-treetype.array[treetypeIndex]
    jobtoload<-paste(treetype,"taxasize", taxasize, ".rda",sep="")
    print(jobtoload)
    if (exists(jobtoload)){
    try(load(jobtoload))
     pngtosave<-paste(treetype,"taxasize", taxasize, "trace.png",sep="")
     png(filename = pngtosave)
     par(mfrow=c(1,3))
     plot(Theta.vec[(n_iter/2):n_iter,"omega"],main="omega")
     plot(Theta.vec[(n_iter/2):n_iter,"alpha"],main="alpha",type="l")
     plot(Theta.vec[(n_iter/2):n_iter,"beta"],type="l",main="beta")
     dev.off()
    #
    pngtosave<-paste(treetype,"taxasize", taxasize, "hist.png",sep="")
    png(filename = pngtosave)
    par(mfrow=c(1,3))
    hist(Theta.vec[(n_iter/2):n_iter,"omega"],main="omega")
    hist(Theta.vec[(n_iter/2):n_iter,"alpha"],main="alpha")
    hist(Theta.vec[(n_iter/2):n_iter,"beta"],main="beta")
    dev.off()
    }

    treetype.array <- c("balanced","left","coalescent","bdtree") #c("bdtree")
  }
}

[1] "balancedtaxasize16.rda"
[1] "lefttaxasize16.rda"
[1] "coalescenttaxasize16.rda"
[1] "bdtreetaxasize16.rda"
[1] "balancedtaxasize32.rda"
[1] "lefttaxasize32.rda"
[1] "coalescenttaxasize32.rda"
[1] "bdtreetaxasize32.rda"
[1] "balancedtaxasize64.rda"
[1] "lefttaxasize64.rda"
[1] "coalescenttaxasize64.rda"
[1] "bdtreetaxasize64.rda"
[1] "balancedtaxasize128.rda"
[1] "lefttaxasize128.rda"
[1] "coalescenttaxasize128.rda"
[1] "bdtreetaxasize128.rda"