#'
#'@title Function to run TCSAM2013.
#'
#'@description This functions runs a TCSAM2013 model once.
#'
#'@details
#'This function creates a shell script ('./tmp.sh') in the
#'working directory and uses it to run a version of the TCSAM2013 model.\cr
#'Initial model parameters can be jittered based on the system clock time as a seed
#'to the random number generator. The seed and final objective function value are
#'saved for each model run in a csv file (the value of out.csv).
#'
#'@param os - 'win' or 'mac' or 'osx'
#'@param path - path for model output
#'@param model - TCSAM2013 model executable name
#'@param path2model - path to model executable
#'@param configFile - filename (including path) to model configuration file
#'@param pin - flag (T/F) to use/not use a default pin file, or an ascii filename to use a non-default one
#'@param hess - T/F to compute hessian (and .std file)
#'@param mcmc - T/F to run mcmc
#'@param minPhase - min phase to start estimation
#'@param maxPhase - max phase for estimation
#'@param jitter - T/F to jitter parameters
#'@param seed - seed for random number generator (or NULL)
#'@param plotResults - T/F to plot results using \code{plotTCSAM2013I.GG}
#'@param rmEXE - flag to remove model executable from run folder
#'@param cleanup - flag (T/F) to clean up unnecessary files
#'
#'@return - dataframe with 2 columns (name, value) with jitter seed (if jittered) and par file info, or NULL if par file does not exist
#'
#'@details If the path associated with \code{configFile} is a relative one, it should
#'be relative to the \code{path} variable.
#'
#'@export
#'
runTCSAM2013<-function(os='osx',
path='.',
model='tcsam2013alta',
path2model='',
configFile='',
pin=FALSE,
hess=FALSE,
mcmc=FALSE,
minPhase=NULL,
maxPhase=NULL,
jitter=FALSE,
seed=NULL,
plotResults=hess,
rmEXE=FALSE,
cleanup=TRUE){
#start timing
stm<-Sys.time();
#switch to run folder (create if necessary)
currdir<-getwd();
on.exit(setwd(currdir));
if (!file.exists(path)) dir.create(path,recursive=TRUE)
setwd(path);
cat("Running tcsam2013 model at '",path,"'.\n");
#set up copy commands
fn.par<-file.path(getwd(),"&&model.par");
fn.par<-gsub('&&model',tolower(model),fn.par)
run.cmds<-getRunCommands(os=os,
model=model,
path2model=path2model,
configFile=configFile,
pin=pin,
hess=hess,
mcmc=mcmc,
minPhase=minPhase,
maxPhase=maxPhase,
jitter=jitter,
seed=seed,
rmEXE=rmEXE,
cleanup=cleanup)
if (tolower(os)=='win'){
cat(run.cmds,file="tmp.bat")
Sys.chmod("tmp.bat",mode='7777')
system("tmp.bat",wait=TRUE);
} else {
cat(run.cmds,file="./tmp.sh")
Sys.chmod("./tmp.sh",mode='7777')
system("./tmp.sh",wait=TRUE);
}
#print timing-related info
etm<-Sys.time();
elt<-etm-stm;
cat("start time: ")
print(stm);
cat("end time: ")
print(etm);
cat("elapsed time: ")
print(elt);
#parse par file into dataframe
par<-paste0(model,'.par')
dfr<-getPar(par);
#get jitter info
if (jitter&(!is.null(dfr))) {
if (file.exists('jitterInfo.csv')){
tbl<-read.csv('jitterInfo.csv',header=TRUE);
dfr<-rbind(data.frame(name='seed',value=tbl$seed[1]),dfr);
}
}
if (plotResults){
lst<-getResLst(inp.dir=getwd());
plotTCSAM2013I.GG(lst,pdf=TRUE);
}
#return dataframe (and return to original folder as working directory)
return(dfr);
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.