Burden test error in ukb23158_500k_OQFE.sets.txt.gz

Elham Alhathli

I got the same error that was reported earlier here

I did the following

for chr in {1..21}; do

run_regenie_burden="regenie --step 2 --pred ap38_results_pred.list --out assoc.c${chr} \

--bed ${data_field}_c${chr}_b0_v1 --ref-first \

--phenoFile Ex1.phe --covarFile Ex1.covar \

--phenoCol Exercise_cc --covarCol age --covarCol sex \

--covarCol genotyping --covarCol PCA{1:10} \

--covarCol Centre --covarCol Month \

--set-list ukb23158_500k_OQFE.sets.txt.gz --anno-file ukb23158_500k_OQFE.annotations.txt.gz --mask-def ldl_custom_masks.txt

--aaf-bins 0.01,0.001 --nauto 23 --bsize 200 --extract-sets extract_gene_names.txt --check-burden-files --threads 16 --gz\

--write-mask-snplist --vc-tests skato,acato-full --out Exercise_skato_acato_chr${i}"

 

dx run swiss-army-knife -iin="${exome_file_dir}/${data_field}_c${chr}_b0_v1.bed" \

-iin="${exome_file_dir}/${data_field}_c${chr}_b0_v1.bim" \

-iin="${exome_file_dir}/${data_field}_c${chr}_b0_v1.fam" \

-iin="${txt_file_dir}/ldl_custom_masks.txt" \

-iin="${txt_file_dir}/Ex1.phe" \

-iin="${txt_file_dir}/Ex1.covar" \

-iin="${txt_file_dir}/extract_gene_names.txt" \

-iin="${data_file_dir}/ap38_results_pred.list" \

-iin="${data_file_dir}/ap38_results_1.loco.gz" \

-iin="${txt_file_dir}/ukb23158_500k_OQFE.sets.txt.gz" \

-iin="${txt_file_dir}/ukb23158_500k_OQFE.annotations.txt.gz" \

-icmd="${run_regenie_burden}" --tag="burden_EX" --instance-type "mem2_ssd1_v2_x96" \

--destination="${project}:/Data/ap_wes_gwas/" --brief --yes


 

done

and I got this error 

* set file         : [ukb23158_500k_OQFE.sets.txt.gz] ERROR: unknown chromosome code in set list file.

I tried to modify the sets file but it's read only file and I can't modify it or make it tab delimited

 

I even tried to run the same script as bash script on Jupyter notebook 

./regenie_v2.2.4.gz_x86_64_Linux_mkl \
   --step 2 \
   --pred "Exercise_cc /mnt/project/Data/ap_wes_gwas/ap38_results_1.loco" \
   --ref-first \
   --bed "/mnt/project/Bulk/Exome sequences/Population level exome OQFE variants, PLINK format - final release/ukb23158_c21_b0_v1" \
   --phenoFile "/mnt/project/gwas_cohort_textfiles/Ex1.phe" \
   --covarFile "/mnt/project/gwas_cohort_textfiles/Ex1.covar" \
   --phenoCol Exercise_cc --covarCol age --covarCol sex \
   --covarCol genotyping --covarCol PCA{1:10} \
   --covarCol Centre --covarCol Month \
   --set-list "/mnt/project/gwas_cohort_textfiles/ukb23158_500k_OQFE.sets.txt.gz" \
   --anno-file "/mnt/project/gwas_cohort_textfiles/ukb23158_500k_OQFE.annotations.txt.gz" \
   --mask-def "/mnt/project/gwas_cohort_textfiles/ukb23158_500k_OQFE.masks" \
   --nauto 23 \
   --aaf-bins 0.01,0.001 \
   --bsize 200 \
   --extract-sets "/mnt/project/gwas_cohort_textfiles/extract_gene_names.txt" \
   --out ldl_pcsk9 \
   --gz

 

and I got the following error * phenotypes       : [/mnt/project/gwas_cohort_textfiles/Ex1.phe] ERROR: cannot read file : Exercise_cc /mnt/project/Data/ap_wes_gwas/ap38_results_1.loco though the file is present in the same directory.

 

 

Comments

0 comments

Please sign in to leave a comment.