Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Help needed for processing #2

Open
zara03 opened this issue May 2, 2024 · 0 comments
Open

Help needed for processing #2

zara03 opened this issue May 2, 2024 · 0 comments

Comments

@zara03
Copy link

zara03 commented May 2, 2024

I have check.sh in ldsc directory where i am running this job script.
My bed file is in directory /home/saz22/Chiou.CRE_bed/GSE163160_T1D_snATAC.cCREs.bed

while 1000G_EUR_Phase3_plink is in directory with 1000G.EUR.QC.1.bim ,these kind of files so I wrote as
/home/saz22/LDSC/ochroGWAS/1000G.EUR.hg38.$val.bim

I am submitting this job script but its getting exit code 2.
I added activate ldsc ( ldsc environment) and added that in script because ldsc has update from python 2 to 3 i have wrote python 3 in script.
Even when i remove activate ldsc job is getting failed. I actually tried with only script provided by you even then its getting failed. I am writing both scripts which are not working,
Can you please look into it and tell me whats the error and how to make it run.

module load python/2.7
module load bedtools/2.30.0
source ~/ENV/bin/activate

vals=($(seq 1 1 22))

for val in ${vals[@]}; do
echo $val
declare -a files= (/home/saz22/Chiou.CRE_bed/GSE163160_T1D_snATAC.cCREs.bed)

for t in "${files[@]}"; do
    declare noPathName="${t##*/}"
    declare noExt="${noPathName%.bed}"
    echo "$noExt"
    make_annot.py \              				
           --bed-file $t \
           --bimfile /home/saz22/LDSC/ochroGWAS/1000G_EUR_Phase3_plink/1000G.EUR.$val.bim \			
           --annot-file /home/saz22/Chiou.CRE_bed$noExt".$val.annot

OR THIS AS WELL:

module purge
module load miniconda
conda activate test

vals=($(seq 1 1 22))

for val in ${vals[@]}; do
echo $val
declare -a files= (/home/saz22/Chiou.CRE_bed/GSE163160_T1D_snATAC.cCREs.bed)

for t in "${files[@]}"; do
    declare noPathName="${t##*/}"
    declare noExt="${noPathName%.bed}"
    echo "$noExt"
    make_annot.py \              				
           --bed-file $t \
           --bimfile /home/saz22/LDSC/ochroGWAS/1000G_EUR_Phase3_plink/1000G.EUR.$val.bim \			
           --annot-file /home/saz22/Chiou.CRE_bed$noExt".$val.annot
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant