File tree Expand file tree Collapse file tree 3 files changed +47
-13
lines changed Expand file tree Collapse file tree 3 files changed +47
-13
lines changed Original file line number Diff line number Diff line change 139
139
140
140
echo " Successfully generated correlation matrices for each phenotype!"
141
141
echo " Generating PRS-phenotype associations for each subset"
142
- Rscript resources/phenotypes/score.r ${phenotype_processed_dir} /phenolist ${genotype_processed_dir} /scratch/tophits ${results_dir} /04
142
+ Rscript resources/phenotypes/score.r ${phenotype_processed_dir} /phenolist ${genotype_processed_dir} /scratch/tophits ${results_dir} /03
143
143
echo " Successfully generated scores!"
Original file line number Diff line number Diff line change 2
2
3
3
source config.env
4
4
5
+
6
+ echo " Checking R packages..."
5
7
Rscript -e " renv::status()"
6
8
7
- ./bin/plink2 --version
8
- ./bin/flashpca --version
9
- ./bin/gcta-1.94.1
10
- ./bin/king
9
+ tf=$( mktemp)
10
+
11
+ echo " Checking plink..."
12
+ if ./bin/plink2 --version > $tf 2>&1 ; then
13
+ echo " All good!"
14
+ else
15
+ cat $tf
16
+ fi
17
+
18
+ echo " Checking flashpca..."
19
+ if ./bin/flashpca --version > $tf 2>&1 ; then
20
+ echo " All good!"
21
+ else
22
+ cat $tf
23
+ fi
24
+
25
+ echo " Checking gcta..."
26
+
27
+ ./bin/gcta-1.94.1 > $tf 2>&1
28
+ if cat $tf | grep -q " Genome" ; then
29
+ echo " All good!"
30
+ else
31
+ cat $tf
32
+ fi
33
+
34
+ echo " Checking gcta..."
35
+
36
+ ./bin/king > $tf 2>&1
37
+ if cat $tf | grep -q " KING" ; then
38
+ echo " All good!"
39
+ else
40
+ cat $tf
41
+ fi
42
+
43
+
Original file line number Diff line number Diff line change 2
2
3
3
source config.env
4
4
5
- apptainer pull -F docker://MRCIEU /lifecourse-gwas:latest
5
+ # apptainer pull -F docker://mrcieu /lifecourse-gwas:latest
6
6
7
7
# Run without pwd binding
8
8
9
- d=${PWD}
10
- echo $d
11
- td=$( mktemp -d)
12
- cd $td
9
+ # d=${PWD}
10
+ # echo $d
11
+ # td=$(mktemp -d)
12
+ # cd $td
13
13
14
14
apptainer run \
15
- --bind " ${d} /config.env:/project/config.env" \
15
+ --containall \
16
+ --bind " ${PWD} /config.env:/project/config.env" \
16
17
--bind " ${phenotype_input_dir} :${phenotype_input_dir} " \
17
18
--bind " ${phenotype_processed_dir} :${phenotype_processed_dir} " \
18
19
--bind " ${genotype_input_dir} :${genotype_input_dir} " \
19
20
--bind " ${genotype_processed_dir} :${genotype_processed_dir} " \
20
21
--bind " ${results_dir} :${results_dir} " \
21
22
--cwd /project \
22
- ${d } /lifecourse-gwas_latest.sif \
23
+ ${PWD } /lifecourse-gwas_latest.sif \
23
24
" $@ "
24
25
25
- cd $d
26
+ # cd $d
You can’t perform that action at this time.
0 commit comments