Mercurial > repos > mingchen0919 > elastic_r_tidyverse
diff elastic_r_package.Rmd @ 9:06f38d3febdc draft default tip
planemo upload
| author | mingchen0919 |
|---|---|
| date | Wed, 28 Mar 2018 23:34:03 -0400 |
| parents | |
| children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/elastic_r_package.Rmd Wed Mar 28 23:34:03 2018 -0400 @@ -0,0 +1,110 @@ +--- +title: 'Tool Report' +output: html_document +--- + +<style> +pre code, pre, code { + white-space: pre !important; + overflow-x: scroll !important; + word-break: keep-all !important; + word-wrap: initial !important; +} +</style> + +```{r setup, include=FALSE, warning=FALSE, message=FALSE} +knitr::opts_chunk$set(error = TRUE) +``` + +## User input + +```{r, echo=FALSE} +knitr::kable(arguments) +``` + + +```{r, echo=FALSE} +# each tool execution runs one or multiple chained functions but generate only one single object. +# we save this object to an rdata file and output this file to galaxy history so that it can be used by other tools +# we can use this rdata output file's dataset id as the variable name of the saved object. +job_script_path = paste0(Sys.getenv('REPORT_FILES_PATH'), '/job-script.R') +tool_rdata_output = Sys.getenv('TOOL_RDATA_OUTPUT') +dataset_id = tail(strsplit(tool_rdata_output, '/')[[1]], 1) +dataset_num = gsub("(.+_)([0-9]+)\\.dat", "\\2", dataset_id) +rdata_id = paste0('rdata_', dataset_num) + +## build script +# the first line of the job script is 'rdata_NUM = ', where 'NUM' is the dataset number of the output rdata. +write(paste0(rdata_id, ' = '), file = job_script_path) +# loop through argument data frame to build up the job script. +for (i in 1: (nrow(arguments)-1)) { + row_type = arguments[i, 'row_type'] + switch (row_type, + # if it's a function row, the line has format 'function_name(' + func = write(paste0(arguments[i, 'function_name'], '('), + file = job_script_path, + append = TRUE ), + + + argument = { + # if it's an argument row and the next row is not an operator row, + # the line has format ' argument_name=argument_value,' + if (arguments[i+1, 'operator'] == "") { + write(paste0(' ', arguments[i, 'argument_name'], '=', arguments[i, 'argument_value'], ','), + file = job_script_path, + append = TRUE ) + } else { + # if it's an argument row and the next row IS an operator row, + # the line has format ' argument_name=argument_value'. note that there is not comma at the end. + write(paste0(' ', arguments[i, 'argument_name'], '=', arguments[i, 'argument_value']), + file = job_script_path, + append = TRUE ) + } + }, + + # if it is an operator row, the line has format ') operator' + operator = write(paste0(') ', arguments[i, 'operator']), + file = job_script_path, + append = TRUE ) + ) +} + +# the last line is missing a ')' +write(')', file = job_script_path, append = TRUE) +``` + + +```{bash, 'display script', results='asis', echo=FALSE} +echo '## Job script' +echo '' +echo '' +echo '```r' +cat ${REPORT_FILES_PATH}/job-script.R +echo '```' +``` + + +## Result + +```{r, 'run job script', echo=FALSE} +source(job_script_path) +# display result. +eval(parse(text = rdata_id)) +``` + + +```{r, 'display output directory contents', results='asis', echo=FALSE} +## after the job is done, we list all files from the output directory. +## full relative path to the output directory needs to be displayed. + +cat('##All output files') +cat('\n\n') +all_files = list.files(path = Sys.getenv('REPORT_FILES_PATH'), + full.names = TRUE, + recursive = TRUE) + +for (f in sub(Sys.getenv('REPORT_FILES_PATH'), '.', all_files) ) { + cat('* [', f, '](', f, ')\n') +} +cat('\n') +``` \ No newline at end of file
