comparison generate_test_data @ 1:ace84aac38ce draft

planemo upload for repository https://github.com/ASaiM/galaxytools/tree/master/tools/humann2/ commit 6753a2c0fec9222f5977cab43a855c63fa9f7bee

0:8ed686bd4229

    --output-format 'tsv' \
    --output-max-decimals '10' \
    --output-basename 'humann2' \
    --pathways 'unipathway'

# humann2_regroup_table
humann2_regroup_table \
    --input 'test-data/demo_genefamilies.tsv' \
    -o 'test-data/regrouped_gene_families_to_infogo1000.tsv' \
    --function 'sum' \

    --function 'max' \
    --sort-by 'name'

# humann2_rename_table
humann2_rename_table \
    --input 'test-data/demo_genefamilies.tsv' \
    -o 'test-data/renamed_genefamilies.tsv' \
    --names 'metacyc-pwy'

# humann2_split_table
mkdir split_table_tmp_dir
humann2_split_table \
    --input 'test-data/joined_pathway_coverage_abundance.tsv' \
    -o 'split_table_tmp_dir'
cp 'split_table_tmp_dir/demo_Abundance.tsv' 'test-data/'
cp 'split_table_tmp_dir/demo_Coverage.tsv' 'test-data/'

1:ace84aac38ce

    --output-format 'tsv' \
    --output-max-decimals '10' \
    --output-basename 'humann2' \
    --pathways 'unipathway'

# humann2_associate
humann2_associate \
    --input 'test-data/hmp_pathabund.txt' \
    --focal-metadatum 'STSite' \
    --last-metadatum 'STSite' \
    --focal-type 'categorical' \
    --fdr '0.2' \
    --output 'test-data/associate_output.txt'

# humann2_barplot
humann2_barplot \
    --input 'test-data/hmp_pathabund.txt' \
    --focal-feature ANAGLYCOLYSIS-PWY \
    --sort none \
    --last-metadatum STSite \
    --focal-metadatum STSite \
    --scaling none \
    --width 5 \
    --dimensions 8 4 \
    --output test-data/barplot1.png

humann2_barplot \
    --input 'test-data/hmp_pathabund.txt' \
    --focal-feature ANAGLYCOLYSIS-PWY \
    --sort sum \
    --last-metadatum STSite \
    --focal-metadatum STSite \
    --scaling pseudolog \
    --width 5 \
    --dimensions 8 4 \
    --output test-data/barplot2.pdf

humann2_barplot \
    --input 'test-data/hmp_pathabund.txt' \
    --focal-feature ANAGLYCOLYSIS-PWY \
    --sort metadata \
    --last-metadatum STSite \
    --focal-metadatum STSite \
    --exclude-unclassified \
    --scaling none \
    --as-genera \
    --grid \
    --remove-zeroes \
    --width 5 \
    --dimensions 8 4 \
    --output test-data/barplot3.svg

humann2_barplot \
    --input 'test-data/hmp_pathabund.txt' \
    --focal-feature ANAGLYCOLYSIS-PWY \
    --top-strata 12 \
    --sort sum \
    --last-metadatum STSite \
    --focal-metadatum STSite \
    --exclude-unclassified \
    --scaling normalize \
    --width 5 \
    --dimensions 8 4 \
    --output test-data/barplot4.eps

# humann2_regroup_table
humann2_regroup_table \
    --input 'test-data/demo_genefamilies.tsv' \
    -o 'test-data/regrouped_gene_families_to_infogo1000.tsv' \
    --function 'sum' \

    --function 'max' \
    --sort-by 'name'

# humann2_rename_table
humann2_rename_table \
    --input 'test-data/regrouped_gene_families_to_infogo1000.tsv' \
    --names 'infogo1000' \
    -o 'test-data/renamed_info1000.tsv'
    
humann2_rename_table \
    --input 'test-data/regrouped_gene_families_to_infogo1000.tsv' \
    --custom 'test-data/rename_mapping' \
    -o 'test-data/renamed_mapping_info1000.tsv'

# humann2_rna_dna_norm
humann2_rna_dna_norm \
    --input_dna 'test-data/rna_dna_norm-dna.txt' \
    --input_rna 'test-data/rna_dna_norm-rna.txt' \
    --output_basename 'results1' \
    --method 'laplace'

humann2_rna_dna_norm \
    --input_dna 'test-data/rna_dna_norm-dna.txt' \
    --input_rna 'test-data/rna_dna_norm-rna.txt' \
    --output_basename 'results2' \
    --method 'witten_bell' \
    --log_transform \
    --log_base '2'

# humann2_split_stratified_table
humann2_split_stratified_table \
    --input 'test-data/demo_genefamilies.tsv' \
    --output 'split_stratified_table'

# humann2_split_table
mkdir split_table_tmp_dir
humann2_split_table \
    --input 'test-data/joined_pathway_coverage_abundance.tsv' \
    -o 'split_table_tmp_dir'
cp 'split_table_tmp_dir/demo_Abundance.tsv' 'test-data/'
cp 'split_table_tmp_dir/demo_Coverage.tsv' 'test-data/'

# humann2_strain_profiler
humann2_strain_profiler \
    --input 'test-data/strain_profiler-input.txt' \
    --critical_mean '1' \
    --critical_count '2' \
    --pinterval '1e-10' '1' \
    --critical_samples '2'

# humann2_unpack_pathways
humann2_unpack_pathways \
    --input-genes 'test-data/demo_genefamilies.tsv' \
    --input-pathways 'test-data/demo_pathabundance.tsv' \
    --output unpack1

humann2_unpack_pathways \
    --input-genes 'test-data/demo_genefamilies.tsv' \
    --input-pathways 'test-data/demo_pathabundance.tsv' \
    --remove-taxonomy \
    --output unpack2