BIMSBbioinfo / pigx_rnaseq

Bulk RNA-seq Data Processing, Quality Control, and Downstream Analysis Pipeline
GNU General Public License v3.0
20 stars 11 forks source link

Code shown in Report does not define "params" #115

Closed smoe closed 2 years ago

smoe commented 2 years ago

To elucidate why the g:profiler results don't show in my hisat2 analysis' report, I thought I should execute the corresponding code prior to creating this issue. However, variables are missing which I presume to be defined somewhere in the code blocks above. The top code block however starts with

countDataFile <- params$countDataFile
colDataFile <- params$colDataFile
gtfFile <- params$gtfFile
caseSampleGroups <- params$caseSampleGroups
controlSampleGroups <- params$controlSampleGroups
covariates <- params$covariates
prefix <- params$prefix
workdir <- params$workdir

which likely refers to the settings.yaml, but the code how to read that is not shown.

borauyar commented 2 years ago

Yes, the Rmarkdown script is executed from another R script, where the parameters defined in the yaml header of the Rmd script are passed via rmarkdown::render. See the line here.

borauyar commented 2 years ago

116 by @smoe improves the report output for debugging purposes.