prstools 0.0.76__tar.gz → 0.0.78__tar.gz
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- {prstools-0.0.76/prstools.egg-info → prstools-0.0.78}/PKG-INFO +26 -48
- {prstools-0.0.76 → prstools-0.0.78}/README.md +25 -47
- {prstools-0.0.76 → prstools-0.0.78}/prstools/__init__.py +1 -1
- {prstools-0.0.76 → prstools-0.0.78}/prstools/_parser_vars.py +2 -2
- {prstools-0.0.76 → prstools-0.0.78}/prstools/models/_base.py +7 -7
- {prstools-0.0.76 → prstools-0.0.78/prstools.egg-info}/PKG-INFO +26 -48
- {prstools-0.0.76 → prstools-0.0.78}/settings.ini +1 -1
- {prstools-0.0.76 → prstools-0.0.78}/LICENSE +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/MANIFEST.in +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/_cmd.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/_ext_utils.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/_modidx.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/EUR/1kg_chr22_22004675_23374984.EUR.edgelist +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/EUR/1kg_chr22_30667654_32269392.EUR.edgelist +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/EUR/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/snplist/1kg_chr22_22004675_23374984.snplist +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/snplist/1kg_chr22_30667654_32269392.snplist +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldgm_1kg_pop/snplist/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldref_1kg_pop/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldref_1kg_pop/ldblk_1kg_chr22.hdf5 +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldref_1kg_pop/snpinfo_1kg_hm3 +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/sumstats.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/target.bed +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/target.bim +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/target.fam +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_1blk_shift=0.regdef.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_2blk_shift=0.regdef.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_2blk_shift=1.regdef.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_3blk_shift=0.regdef.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_3blk_shift=1.regdef.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_3blk_shift=2.regdef.tsv +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/errors.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/io.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/linkage/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/linkage/_base.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/models/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/models/_compute.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/parse_genet.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/scores.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/tests/__init__.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/tests/test_cli_model_examples.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools/utils.py +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools.egg-info/SOURCES.txt +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools.egg-info/dependency_links.txt +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools.egg-info/entry_points.txt +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools.egg-info/not-zip-safe +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools.egg-info/requires.txt +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/prstools.egg-info/top_level.txt +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/setup.cfg +0 -0
- {prstools-0.0.76 → prstools-0.0.78}/setup.py +0 -0
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
Metadata-Version: 2.1
|
|
2
2
|
Name: prstools
|
|
3
|
-
Version: 0.0.
|
|
3
|
+
Version: 0.0.78
|
|
4
4
|
Summary: Convenient and powerfull Polygenic Risk Score creation.
|
|
5
5
|
Home-page: https://github.com/mennojw/prstools-release
|
|
6
6
|
Author: Menno Witteveen et al.
|
|
@@ -43,7 +43,7 @@ Requires-Dist: seaborn; extra == "full"
|
|
|
43
43
|
<!-- WARNING: THIS FILE WAS AUTOGENERATED! DO NOT EDIT! -->
|
|
44
44
|
|
|
45
45
|
`prstools` is software to create Polygenic Risk Scores (PRS) directly
|
|
46
|
-
from the commandline <br>
|
|
46
|
+
from the commandline. <br>
|
|
47
47
|
|
|
48
48
|
It makes PRS generation easier, compared to previous tools, by:
|
|
49
49
|
|
|
@@ -94,62 +94,40 @@ prstools prscs2 --ref ldref_1kg_pop --target target \
|
|
|
94
94
|
This will run PRS-CS2 on the example data, using the new implementation
|
|
95
95
|
to demonstrate the capabilities of `prstools` and makes PRS predictions
|
|
96
96
|
for the example dataset. The best and fastest way to get a PRS for your
|
|
97
|
-
is to try the **Tutorial** below. <br>
|
|
98
|
-
documentation, which resides inside of the command-line interface, which
|
|
99
|
-
you can see by typing `prstools`.
|
|
100
|
-
|
|
101
|
-
``` console
|
|
102
|
-
Usage:
|
|
103
|
-
prstools <command> ...
|
|
104
|
-
|
|
105
|
-
Convenient and powerfull Polygenic Risk Score creation.
|
|
106
|
-
'prst' is a commandline shorthand for 'prstools'
|
|
107
|
-
|
|
108
|
-
Models & Utility Commands:
|
|
109
|
-
<command>
|
|
110
|
-
downloadutil Download and unpack LD reference panels and other data.
|
|
111
|
-
combine A tool to combine genetics-related text files.
|
|
112
|
-
prscs2 PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes
|
|
113
|
-
under continuous shrinkage (CS) priors.
|
|
114
|
-
```
|
|
97
|
+
case is to try the **Tutorial** below. <br>
|
|
115
98
|
|
|
116
|
-
|
|
117
|
-
|
|
118
|
-
`prstools prscs2` will output the
|
|
119
|
-
subcommand:
|
|
99
|
+
There is also the `prstools` documentation, residing inside of the
|
|
100
|
+
command-line interface, which you can see by typing `prstools` or a
|
|
101
|
+
subcommand. <br> Forinstance typing `prstools prscs2` will output the
|
|
102
|
+
**documentation** for the `prscs2` subcommand:
|
|
120
103
|
|
|
121
|
-
|
|
122
|
-
|
|
123
|
-
prstools prscs2 [-h --cpus <number-of-cpus>] --ref <dir/refcode> --target <bim-prefix>
|
|
104
|
+
Usage:
|
|
105
|
+
prst prscs2 [-h --cpus <num-of-cpus>] --ref <dir/refcode> --target <bim-prefix>
|
|
124
106
|
--sst <file> --out <dir+prefix> [--n_gwas <num> --chrom <chroms>]
|
|
125
|
-
[--colmap <
|
|
126
|
-
[--
|
|
127
|
-
[--b <b> --phi <phi> --clip <clip>
|
|
128
|
-
|
|
129
|
-
PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under continuous shrinkage (CS) priors.
|
|
107
|
+
[--colmap <colnames> --rsidmode <yes/no> --pred <yes/no>]
|
|
108
|
+
[--n_iter <n_iter> --n_burnin <n_burnin> --n_slice <n_slice>]
|
|
109
|
+
[--seed <seed> --a <a> --b <b> --phi <phi> --clip <clip>]
|
|
110
|
+
[--sampler <sampler> --n_jobs <n_jobs>]
|
|
130
111
|
|
|
131
|
-
|
|
132
|
-
|
|
133
|
-
-c, --cpus <number-of-cpus> The number of cpus to use. Generally most efficient if
|
|
134
|
-
chosen to be between 1 and 5. Functionality can be
|
|
112
|
+
PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under
|
|
113
|
+
continuous shrinkage (CS) priors.
|
|
135
114
|
|
|
136
|
-
|
|
115
|
+
General Options:
|
|
116
|
+
-h, --help Show this help message and exit.
|
|
137
117
|
|
|
118
|
+
... [omitted for readability] ...
|
|
138
119
|
|
|
139
|
-
|
|
140
|
-
|
|
141
|
-
--sampler <sampler> Sampler algorithm. Rue sampling is the original sampler,
|
|
142
|
-
which gives good results. (default: Rue)
|
|
120
|
+
--n_jobs <n_jobs> This sets the number of jobs for parallel processing. (default:
|
|
121
|
+
8)
|
|
143
122
|
|
|
144
|
-
Examples --> can be directly copy-pasted (:
|
|
145
|
-
|
|
146
|
-
|
|
147
|
-
prst prscs2 -r ldref_1kg_pop -t target -s sumstats.tsv -n 2565 -o ./result-prscs2 --pred # A shorter version of previous that also does the predictions.
|
|
148
|
-
```
|
|
123
|
+
# Examples (get data, run model) --> can be directly copy-pasted (:
|
|
124
|
+
prstools downloadutil --pattern example --destdir ./; cd example
|
|
125
|
+
prst prscs2 --ref ldref_1kg_pop -t target -s sumstats.tsv --n_gwas 2565 --out result-prscs2
|
|
149
126
|
|
|
150
127
|
As can be seen, there are examples at the end of the help output to
|
|
151
|
-
illustrate usage, which should work with a simple copy-paste.
|
|
152
|
-
|
|
128
|
+
illustrate usage, which should work with a simple copy-paste.
|
|
129
|
+
|
|
130
|
+
There is now also an online version of all this documentation
|
|
153
131
|
(https://prstools.readthedocs.io/). <br>
|
|
154
132
|
|
|
155
133
|
## Tutorial + Video
|
|
@@ -4,7 +4,7 @@
|
|
|
4
4
|
<!-- WARNING: THIS FILE WAS AUTOGENERATED! DO NOT EDIT! -->
|
|
5
5
|
|
|
6
6
|
`prstools` is software to create Polygenic Risk Scores (PRS) directly
|
|
7
|
-
from the commandline <br>
|
|
7
|
+
from the commandline. <br>
|
|
8
8
|
|
|
9
9
|
It makes PRS generation easier, compared to previous tools, by:
|
|
10
10
|
|
|
@@ -55,62 +55,40 @@ prstools prscs2 --ref ldref_1kg_pop --target target \
|
|
|
55
55
|
This will run PRS-CS2 on the example data, using the new implementation
|
|
56
56
|
to demonstrate the capabilities of `prstools` and makes PRS predictions
|
|
57
57
|
for the example dataset. The best and fastest way to get a PRS for your
|
|
58
|
-
is to try the **Tutorial** below. <br>
|
|
59
|
-
documentation, which resides inside of the command-line interface, which
|
|
60
|
-
you can see by typing `prstools`.
|
|
61
|
-
|
|
62
|
-
``` console
|
|
63
|
-
Usage:
|
|
64
|
-
prstools <command> ...
|
|
65
|
-
|
|
66
|
-
Convenient and powerfull Polygenic Risk Score creation.
|
|
67
|
-
'prst' is a commandline shorthand for 'prstools'
|
|
68
|
-
|
|
69
|
-
Models & Utility Commands:
|
|
70
|
-
<command>
|
|
71
|
-
downloadutil Download and unpack LD reference panels and other data.
|
|
72
|
-
combine A tool to combine genetics-related text files.
|
|
73
|
-
prscs2 PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes
|
|
74
|
-
under continuous shrinkage (CS) priors.
|
|
75
|
-
```
|
|
58
|
+
case is to try the **Tutorial** below. <br>
|
|
76
59
|
|
|
77
|
-
|
|
78
|
-
|
|
79
|
-
`prstools prscs2` will output the
|
|
80
|
-
subcommand:
|
|
60
|
+
There is also the `prstools` documentation, residing inside of the
|
|
61
|
+
command-line interface, which you can see by typing `prstools` or a
|
|
62
|
+
subcommand. <br> Forinstance typing `prstools prscs2` will output the
|
|
63
|
+
**documentation** for the `prscs2` subcommand:
|
|
81
64
|
|
|
82
|
-
|
|
83
|
-
|
|
84
|
-
prstools prscs2 [-h --cpus <number-of-cpus>] --ref <dir/refcode> --target <bim-prefix>
|
|
65
|
+
Usage:
|
|
66
|
+
prst prscs2 [-h --cpus <num-of-cpus>] --ref <dir/refcode> --target <bim-prefix>
|
|
85
67
|
--sst <file> --out <dir+prefix> [--n_gwas <num> --chrom <chroms>]
|
|
86
|
-
[--colmap <
|
|
87
|
-
[--
|
|
88
|
-
[--b <b> --phi <phi> --clip <clip>
|
|
89
|
-
|
|
90
|
-
PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under continuous shrinkage (CS) priors.
|
|
68
|
+
[--colmap <colnames> --rsidmode <yes/no> --pred <yes/no>]
|
|
69
|
+
[--n_iter <n_iter> --n_burnin <n_burnin> --n_slice <n_slice>]
|
|
70
|
+
[--seed <seed> --a <a> --b <b> --phi <phi> --clip <clip>]
|
|
71
|
+
[--sampler <sampler> --n_jobs <n_jobs>]
|
|
91
72
|
|
|
92
|
-
|
|
93
|
-
|
|
94
|
-
-c, --cpus <number-of-cpus> The number of cpus to use. Generally most efficient if
|
|
95
|
-
chosen to be between 1 and 5. Functionality can be
|
|
73
|
+
PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under
|
|
74
|
+
continuous shrinkage (CS) priors.
|
|
96
75
|
|
|
97
|
-
|
|
76
|
+
General Options:
|
|
77
|
+
-h, --help Show this help message and exit.
|
|
98
78
|
|
|
79
|
+
... [omitted for readability] ...
|
|
99
80
|
|
|
100
|
-
|
|
101
|
-
|
|
102
|
-
--sampler <sampler> Sampler algorithm. Rue sampling is the original sampler,
|
|
103
|
-
which gives good results. (default: Rue)
|
|
81
|
+
--n_jobs <n_jobs> This sets the number of jobs for parallel processing. (default:
|
|
82
|
+
8)
|
|
104
83
|
|
|
105
|
-
Examples --> can be directly copy-pasted (:
|
|
106
|
-
|
|
107
|
-
|
|
108
|
-
prst prscs2 -r ldref_1kg_pop -t target -s sumstats.tsv -n 2565 -o ./result-prscs2 --pred # A shorter version of previous that also does the predictions.
|
|
109
|
-
```
|
|
84
|
+
# Examples (get data, run model) --> can be directly copy-pasted (:
|
|
85
|
+
prstools downloadutil --pattern example --destdir ./; cd example
|
|
86
|
+
prst prscs2 --ref ldref_1kg_pop -t target -s sumstats.tsv --n_gwas 2565 --out result-prscs2
|
|
110
87
|
|
|
111
88
|
As can be seen, there are examples at the end of the help output to
|
|
112
|
-
illustrate usage, which should work with a simple copy-paste.
|
|
113
|
-
|
|
89
|
+
illustrate usage, which should work with a simple copy-paste.
|
|
90
|
+
|
|
91
|
+
There is now also an online version of all this documentation
|
|
114
92
|
(https://prstools.readthedocs.io/). <br>
|
|
115
93
|
|
|
116
94
|
## Tutorial + Video
|
|
@@ -150,8 +150,8 @@ def get_subparserkwg_lst():
|
|
|
150
150
|
'display_info': True,
|
|
151
151
|
'help': 'PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under continuous shrinkage (CS) priors.',
|
|
152
152
|
'epilog': '\x1b[32m# Examples (get data, run model) --> can be directly copy-pasted (:\x1b[0m\n'
|
|
153
|
-
'
|
|
154
|
-
'
|
|
153
|
+
'prstools downloadutil --pattern example --destdir ./; cd example \n'
|
|
154
|
+
'prst prscs2 --ref ldref_1kg_pop -t target -s sumstats.tsv --n_gwas 2565 --out result-prscs2\n',
|
|
155
155
|
'modulename': 'prstools.models._base',
|
|
156
156
|
'groups': {'general': {'grpheader': 'General Options',
|
|
157
157
|
'pkwargs': {'basics': {'args': ['-h', '--help'], 'kwargs': {'action': 'help', 'help': 'Show this help message and exit.'}},
|
|
@@ -87,8 +87,8 @@ class BasePred(ABC):
|
|
|
87
87
|
|
|
88
88
|
epilog=f'''\
|
|
89
89
|
# Examples (get data, run model) --> can be directly copy-pasted (:
|
|
90
|
-
|
|
91
|
-
|
|
90
|
+
prstools downloadutil --pattern example --destdir ./; cd example {insert}
|
|
91
|
+
prst {cmdname} --ref {ldrefname} -t target -s sumstats.tsv {chromopt}--n_gwas 2565 --out result-{cmdname}
|
|
92
92
|
'''
|
|
93
93
|
|
|
94
94
|
# prst {cmdname} -r {ldrefname} -t target -s sumstats.tsv -n 2565 {chromopt}-o ./result-{cmdname} --pred # A shorter version of previous that also does the predictions.
|
|
@@ -209,7 +209,7 @@ class BasePred(ABC):
|
|
|
209
209
|
return pkwargs
|
|
210
210
|
|
|
211
211
|
@classmethod
|
|
212
|
-
def _get_linkageclass(pkwargs=None, ref=None):
|
|
212
|
+
def _get_linkageclass(cls, pkwargs=None, ref=None):
|
|
213
213
|
if pkwargs is None: pkwargs = cls._get_pkwargs_for_class(cls)
|
|
214
214
|
try: from prstools.linkage import AutoLinkageData as linkcls
|
|
215
215
|
except: from prstools.linkage import RefLinkageData as linkcls
|
|
@@ -246,11 +246,11 @@ class BasePred(ABC):
|
|
|
246
246
|
if pred and pred != 'no': # Prediction
|
|
247
247
|
try:
|
|
248
248
|
bed = prst.io.load_bed(target, verbose=verbose);
|
|
249
|
-
yhat = model.predict(bed, rsidmode=rsidmode)
|
|
249
|
+
yhat = model.predict(bed, rsidmode=rsidmode)
|
|
250
250
|
prst.io.save_prs(yhat, fn=out_fnfmt, verbose=verbose); ysv=True # Store prediction result (ysv is helper var, to see if step finished)
|
|
251
|
-
pheno = prst.io.load_pheno(target, verbose=verbose)
|
|
252
|
-
scores = prst.scores.eval(pheno, yhat, metrics=['R2','AUC','etc'], verbose=verbose)
|
|
253
|
-
prst.io.save_scores(scores, fn=out_fnfmt, verbose=verbose)
|
|
251
|
+
#pheno = prst.io.load_pheno(target, verbose=verbose)
|
|
252
|
+
#scores = prst.scores.eval(pheno, yhat, metrics=['R2','AUC','etc'], verbose=verbose)
|
|
253
|
+
#prst.io.save_scores(scores, fn=out_fnfmt, verbose=verbose)
|
|
254
254
|
except Exception as e: # One could have some remarks about the logic of this section, but Menno did not want an if/else jungle here.
|
|
255
255
|
inject = 'evaluation' if ysv else 'prediction'
|
|
256
256
|
msg = (f"Could not generate {inject} (e.g. plink/pheno file missing)"
|
|
@@ -1,6 +1,6 @@
|
|
|
1
1
|
Metadata-Version: 2.1
|
|
2
2
|
Name: prstools
|
|
3
|
-
Version: 0.0.
|
|
3
|
+
Version: 0.0.78
|
|
4
4
|
Summary: Convenient and powerfull Polygenic Risk Score creation.
|
|
5
5
|
Home-page: https://github.com/mennojw/prstools-release
|
|
6
6
|
Author: Menno Witteveen et al.
|
|
@@ -43,7 +43,7 @@ Requires-Dist: seaborn; extra == "full"
|
|
|
43
43
|
<!-- WARNING: THIS FILE WAS AUTOGENERATED! DO NOT EDIT! -->
|
|
44
44
|
|
|
45
45
|
`prstools` is software to create Polygenic Risk Scores (PRS) directly
|
|
46
|
-
from the commandline <br>
|
|
46
|
+
from the commandline. <br>
|
|
47
47
|
|
|
48
48
|
It makes PRS generation easier, compared to previous tools, by:
|
|
49
49
|
|
|
@@ -94,62 +94,40 @@ prstools prscs2 --ref ldref_1kg_pop --target target \
|
|
|
94
94
|
This will run PRS-CS2 on the example data, using the new implementation
|
|
95
95
|
to demonstrate the capabilities of `prstools` and makes PRS predictions
|
|
96
96
|
for the example dataset. The best and fastest way to get a PRS for your
|
|
97
|
-
is to try the **Tutorial** below. <br>
|
|
98
|
-
documentation, which resides inside of the command-line interface, which
|
|
99
|
-
you can see by typing `prstools`.
|
|
100
|
-
|
|
101
|
-
``` console
|
|
102
|
-
Usage:
|
|
103
|
-
prstools <command> ...
|
|
104
|
-
|
|
105
|
-
Convenient and powerfull Polygenic Risk Score creation.
|
|
106
|
-
'prst' is a commandline shorthand for 'prstools'
|
|
107
|
-
|
|
108
|
-
Models & Utility Commands:
|
|
109
|
-
<command>
|
|
110
|
-
downloadutil Download and unpack LD reference panels and other data.
|
|
111
|
-
combine A tool to combine genetics-related text files.
|
|
112
|
-
prscs2 PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes
|
|
113
|
-
under continuous shrinkage (CS) priors.
|
|
114
|
-
```
|
|
97
|
+
case is to try the **Tutorial** below. <br>
|
|
115
98
|
|
|
116
|
-
|
|
117
|
-
|
|
118
|
-
`prstools prscs2` will output the
|
|
119
|
-
subcommand:
|
|
99
|
+
There is also the `prstools` documentation, residing inside of the
|
|
100
|
+
command-line interface, which you can see by typing `prstools` or a
|
|
101
|
+
subcommand. <br> Forinstance typing `prstools prscs2` will output the
|
|
102
|
+
**documentation** for the `prscs2` subcommand:
|
|
120
103
|
|
|
121
|
-
|
|
122
|
-
|
|
123
|
-
prstools prscs2 [-h --cpus <number-of-cpus>] --ref <dir/refcode> --target <bim-prefix>
|
|
104
|
+
Usage:
|
|
105
|
+
prst prscs2 [-h --cpus <num-of-cpus>] --ref <dir/refcode> --target <bim-prefix>
|
|
124
106
|
--sst <file> --out <dir+prefix> [--n_gwas <num> --chrom <chroms>]
|
|
125
|
-
[--colmap <
|
|
126
|
-
[--
|
|
127
|
-
[--b <b> --phi <phi> --clip <clip>
|
|
128
|
-
|
|
129
|
-
PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under continuous shrinkage (CS) priors.
|
|
107
|
+
[--colmap <colnames> --rsidmode <yes/no> --pred <yes/no>]
|
|
108
|
+
[--n_iter <n_iter> --n_burnin <n_burnin> --n_slice <n_slice>]
|
|
109
|
+
[--seed <seed> --a <a> --b <b> --phi <phi> --clip <clip>]
|
|
110
|
+
[--sampler <sampler> --n_jobs <n_jobs>]
|
|
130
111
|
|
|
131
|
-
|
|
132
|
-
|
|
133
|
-
-c, --cpus <number-of-cpus> The number of cpus to use. Generally most efficient if
|
|
134
|
-
chosen to be between 1 and 5. Functionality can be
|
|
112
|
+
PRS-CS v2: A polygenic prediction method that infers posterior SNP effect sizes under
|
|
113
|
+
continuous shrinkage (CS) priors.
|
|
135
114
|
|
|
136
|
-
|
|
115
|
+
General Options:
|
|
116
|
+
-h, --help Show this help message and exit.
|
|
137
117
|
|
|
118
|
+
... [omitted for readability] ...
|
|
138
119
|
|
|
139
|
-
|
|
140
|
-
|
|
141
|
-
--sampler <sampler> Sampler algorithm. Rue sampling is the original sampler,
|
|
142
|
-
which gives good results. (default: Rue)
|
|
120
|
+
--n_jobs <n_jobs> This sets the number of jobs for parallel processing. (default:
|
|
121
|
+
8)
|
|
143
122
|
|
|
144
|
-
Examples --> can be directly copy-pasted (:
|
|
145
|
-
|
|
146
|
-
|
|
147
|
-
prst prscs2 -r ldref_1kg_pop -t target -s sumstats.tsv -n 2565 -o ./result-prscs2 --pred # A shorter version of previous that also does the predictions.
|
|
148
|
-
```
|
|
123
|
+
# Examples (get data, run model) --> can be directly copy-pasted (:
|
|
124
|
+
prstools downloadutil --pattern example --destdir ./; cd example
|
|
125
|
+
prst prscs2 --ref ldref_1kg_pop -t target -s sumstats.tsv --n_gwas 2565 --out result-prscs2
|
|
149
126
|
|
|
150
127
|
As can be seen, there are examples at the end of the help output to
|
|
151
|
-
illustrate usage, which should work with a simple copy-paste.
|
|
152
|
-
|
|
128
|
+
illustrate usage, which should work with a simple copy-paste.
|
|
129
|
+
|
|
130
|
+
There is now also an online version of all this documentation
|
|
153
131
|
(https://prstools.readthedocs.io/). <br>
|
|
154
132
|
|
|
155
133
|
## Tutorial + Video
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/_example/ldref_1kg_pop/ldblk_1kg_chr22.hdf5
RENAMED
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_1blk_shift=0.regdef.tsv
RENAMED
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_2blk_shift=0.regdef.tsv
RENAMED
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_2blk_shift=1.regdef.tsv
RENAMED
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_3blk_shift=0.regdef.tsv
RENAMED
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_3blk_shift=1.regdef.tsv
RENAMED
|
File without changes
|
{prstools-0.0.76 → prstools-0.0.78}/prstools/data/defs/regdef/regions_3blk_shift=2.regdef.tsv
RENAMED
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|
|
File without changes
|