/working/auteur.extended/R/cleanup.files.R

http://github.com/eastman/auteur · R · 24 lines · 18 code · 4 blank · 2 comment · 3 complexity · 5c9ff91d7e133923fe9c737f22dc11af MD5 · raw file 

    

  1. #utility for converting text files to .rda to compress output from rjmcmc
    2. 

  2. #author: JM EASTMAN 2011
    3. 

  3. 

  4. cleanup.files <-
    5. 

  5. function(parmBase, model, fileBase){
    6. 

  6. 	if(model=="BM") {
    7. 

  7. 		parms=c("rates","rate.shifts")
    8. 

  8. 	} else if(model=="OU") {
    9. 

  9. 		parms=c("optima","optima.shifts")
    10. 

  10. 	}
    11. 

  11. 	
    12. 

  12. 	posteriorsamples=lapply(parms, function(x) {
    13. 

  13. 							y=read.table(paste(parmBase,paste(x,"txt",sep="."),sep="/"))
    14. 

  14. 							names(y)=y[1,]
    15. 

  15. 							y=y[-1,]
    16. 

  16. 							row.names(y)=1:nrow(y)
    17. 

  17. 							return(y)
    18. 

  18. 							})
    19. 

  19. 	
    20. 

  20. 	names(posteriorsamples)=parms
    21. 

  21. 	save(posteriorsamples, file=paste(parmBase,paste(fileBase,"posteriorsamples","rda",sep="."),sep="/"))
    22. 

  22. 	unlink(paste(parmBase,paste(parms,"txt",sep="."),sep="/"))
    23. 

  23. }
    24.