@@ -331,9 +331,9 @@ workflow_mode="complete" # Select any of the modes described here: http
331331lai_tool=" rfmix2" # not required if "phasing_only" or "tractor_only"
332332
333333# # SHAPEIT5 (Mandatory Arguments)
334- shapeit5_input_vcf=" /path/to/input_QC .vcf.gz"
335- shapeit5_chunkfile=" /path/to/chunks_chr${chr} .txt"
336- shapeit5_genetic_map=" /path/to/maps/b37/chr${chr} .b37.gmap.gz"
334+ shapeit5_input_vcf=" /path/to/test_data/admixed_cohort/ASW.unphased .vcf.gz"
335+ shapeit5_chunkfile=" /path/to/TractorWorkflow/resources/genomic_chunks/chunks_fullchromosome/ chunks_chr${chr} .txt"
336+ shapeit5_genetic_map=" /path/to/shapeit5/resources/ maps/b37/chr${chr} .b37.gmap.gz"
337337
338338# # SHAPEIT5 (Optional Arguments)
339339# shapeit5_ref_vcf=""
@@ -347,9 +347,9 @@ shapeit5_genetic_map="/path/to/maps/b37/chr${chr}.b37.gmap.gz"
347347# NOTE: If you need additional (optional) arguments, add them in Step 5 as well.
348348
349349# # RFMix2 (Mandatory Arguments)
350- rfmix2_ref_vcf=" /path/to/reference/reference_chr ${chr} _phased .vcf.gz"
351- rfmix2_sample_map=" /path/to/reference/TGP_HGDP_2way_AFR_EUR .txt"
352- rfmix2_genetic_map=" /path/to/maps/b37/shapeit5_genetic_map_LAIformat1 .txt"
350+ rfmix2_ref_vcf=" /path/to/test_data/references/TGP_HGDP_QC_hg19_chr ${chr} .vcf.gz"
351+ rfmix2_sample_map=" /path/to/test_data/references/YRI_GBR_samplemap .txt"
352+ rfmix2_genetic_map=" /path/to/TractorWorkflow/resources/genetic_maps/shapeit5_genetic_map_b37_LAIformat1 .txt"
353353
354354# # RFMix2 (Optional Arguments)
355355# rfmix2_reanalyze_ref=""
@@ -362,15 +362,15 @@ rfmix2_genetic_map="/path/to/maps/b37/shapeit5_genetic_map_LAIformat1.txt"
362362# # GNomix (Mandatory Arguments)
363363# gnomix_dir="/path/software/gnomix"
364364# gnomix_config="${gnomix_dir}/config.yaml"
365- # gnomix_ref_vcf="/path/reference/reference_chr ${chr}_phased .vcf.gz"
366- # gnomix_sample_map="/path/reference/TGP_HGDP_2way_AFR_EUR .txt"
367- # gnomix_genetic_map="/path/maps/b37/shapeit5_genetic_map_LAIformat1 .txt"
365+ # gnomix_ref_vcf="/path/to/test_data/references/TGP_HGDP_QC_hg19_chr ${chr}.vcf.gz"
366+ # gnomix_sample_map="/path/to/test_data/references/YRI_GBR_samplemap .txt"
367+ # gnomix_genetic_map="/path/to/TractorWorkflow/resources/genetic_maps/shapeit5_genetic_map_b37_LAIformat1 .txt"
368368# gnomix_phase="false"
369369
370370# # FLARE (Mandatory Arguments)
371371# flare_dir="/path/softwares/flare"
372- # flare_ref_vcf="/path/reference /TGP_HGDP_QC_hg19_chr${chr}.vcf.gz"
373- # flare_sample_map="/path/reference /YRI_GBR_samplemap.txt"
372+ # flare_ref_vcf="/path/to/test_data/references /TGP_HGDP_QC_hg19_chr${chr}.vcf.gz"
373+ # flare_sample_map="/path/to/test_data/references /YRI_GBR_samplemap.txt"
374374# flare_genetic_map="/path/beagle_genetic_map_files/plink.chr${chr}.GRCh37.map"
375375
376376# FLARE (Optional Arguments)
@@ -395,13 +395,13 @@ extracts_num_ancs=2
395395# extracts_compress=""
396396
397397# # Tractor (Mandatory Arguments)
398- tractor_phenotype=" /path/to/phenotype_files/pheno_irnt .txt"
399- tractor_covarcollist=" age,sex,PC1,PC2,PC3,PC4,PC5 "
398+ tractor_phenotype=" /path/to/test_data/phenotype/Phe_linear_covars_mod1 .txt"
399+ tractor_covarcollist=" age,sex"
400400tractor_regression_method=" linear"
401401
402402# # Tractor (one of the following arguments is Mandatory, unless "y" is used as phenocol bc that's default)
403403# tractor_phenocol=""
404- tractor_phenolist_file=" /path/to/pheno_irnt_phenolist .txt"
404+ tractor_phenolist_file=" /path/to/test_data/phenotype/Phe_linear_covars_mod1_phenolist .txt"
405405
406406# # Tractor (Optional Arguments)
407407# tractor_sampleidcol=""
@@ -447,7 +447,7 @@ Note that the SLURM directives are used to launch the Tractor workflow, which in
447447| ` --partition=<ATTN: add-partition-name> ` | Specifies the partition or queue for job submission. |
448448| ` --job-name=launch_nxf_tractor ` | Custom name for the job. |
449449| ` --output=launch_nxf_tractor_%a_%j.out ` | File to store job logs, ` %a ` : array index (chr no.) and ` %j ` : job ID.|
450- | ` --array=1-22 ` | Runs the job as an array i.e. 22 jobs are launched, with one task per chromosome (1–22). |
450+ | ` --array=1-22 ` | Runs the job as an array i.e. 22 jobs are launched, with one task per chromosome (1–22). |
451451
452452** Explanation of Nextflow directives**
453453
0 commit comments