R/topologyprofile.R

"topologyprofile"<-function(basename){
   seqlength<-scan(paste(basename,".phy",sep=""),nlines=1,quiet= TRUE)[2]
   tps <- rJava::.jnew("TopologyProfile", paste(basename,".post",sep=""),paste(basename,".topoprob",sep=""),paste(basename,".tree",sep=""),as.integer(seqlength))
   #system(paste("java TopologyProfile ",basename,".post ",basename,".topoprob ",basename,".tree ",seqlength,sep=""))
}
vnminin/rbrothers documentation built on May 3, 2019, 6:38 p.m.