Mercurial > repos > ethevenot > multivariate
comparison test/test-multi @ 4:5526f8258e8a draft default tip
planemo upload for repository https://github.com/workflow4metabolomics/multivariate.git commit 0f382a5296aae9bfc77df06b0a5ad493eb3c01f3
| author | ethevenot |
|---|---|
| date | Wed, 28 Feb 2018 09:59:25 -0500 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| 3:e91de3b04320 | 4:5526f8258e8a |
|---|---|
| 1 #!/bin/bash | |
| 2 | |
| 3 # Constants {{{1 | |
| 4 ################################################################ | |
| 5 | |
| 6 PROG_PATH=$(dirname $0) | |
| 7 | |
| 8 # MAIN {{{1 | |
| 9 ################################################################ | |
| 10 | |
| 11 $PROG_PATH/../multivariate_wrapper.R dataMatrix_in $PROG_PATH/../test-data/input-dataMatrix.tsv sampleMetadata_in $PROG_PATH/../test-data/input-sampleMetadata.tsv variableMetadata_in $PROG_PATH/../test-data/input-variableMetadata.tsv respC age predI 1 orthoI 1 testL FALSE sampleMetadata_out $PROG_PATH/outputSampleMetadata.tsv variableMetadata_out $PROG_PATH/outputVariableMetadata.tsv | |
| 12 | |
| 13 # Computed numbers are not always the same. We cannot use diff directly on files. | |
| 14 | |
| 15 # 0) Check they have the same number of lines | |
| 16 nlines=$(wc -l <"$PROG_PATH/outputSampleMetadata.tsv") | |
| 17 ref_nlines=$(wc -l <"$PROG_PATH/../test-data/output-sampleMetadata.tsv") | |
| 18 if [ $nlines != $ref_nlines ] ; then | |
| 19 echo "Incorrect output sample metadata." >&2 | |
| 20 exit 1 | |
| 21 fi | |
| 22 | |
| 23 # 1) We check the header line | |
| 24 head -n 1 "$PROG_PATH/outputSampleMetadata.tsv" > "$PROG_PATH/header.tsv" | |
| 25 head -n 1 "$PROG_PATH/../test-data/output-sampleMetadata.tsv" > "$PROG_PATH/ref_header.tsv" | |
| 26 if ! diff "$PROG_PATH/header.tsv" "$PROG_PATH/ref_header.tsv" ; then | |
| 27 echo "Incorrect output sample metadata." >&2 | |
| 28 exit 1 | |
| 29 fi | |
| 30 | |
| 31 # 2) We remove computed values and compare what's left | |
| 32 awk 'BEGIN{FS=OFS="\t"} {$7=$8=$9="";sub("\t\t","\t")}1' "$PROG_PATH/outputSampleMetadata.tsv" > "$PROG_PATH/nocomput.tsv" | |
| 33 awk 'BEGIN{FS=OFS="\t"} {$7=$8=$9="";sub("\t\t","\t")}1' "$PROG_PATH/../test-data/output-sampleMetadata.tsv" > "$PROG_PATH/ref_nocomput.tsv" | |
| 34 if ! diff "$PROG_PATH/nocomput.tsv" "$PROG_PATH/ref_nocomput.tsv" ; then | |
| 35 echo "Incorrect output sample metadata." >&2 | |
| 36 exit 1 | |
| 37 fi |
