gwaslab 3.6.6__py3-none-any.whl → 3.6.8__py3-none-any.whl
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.
Potentially problematic release.
This version of gwaslab might be problematic. Click here for more details.
- gwaslab/__init__.py +57 -47
- gwaslab/{bd_common_data.py → bd/bd_common_data.py} +10 -9
- gwaslab/bd/bd_config.py +28 -0
- gwaslab/{bd_download.py → bd/bd_download.py} +1 -1
- gwaslab/{bd_get_hapmap3.py → bd/bd_get_hapmap3.py} +9 -6
- gwaslab/bd/bd_path_manager.py +110 -0
- gwaslab/data/formatbook.json +805 -9
- gwaslab/{ldsc_irwls.py → extension/ldsc/ldsc_irwls.py} +1 -1
- gwaslab/{ldsc_regressions.py → extension/ldsc/ldsc_regressions.py} +2 -2
- gwaslab/{ldsc_sumstats.py → extension/ldsc/ldsc_sumstats.py} +2 -2
- gwaslab/{prscs_mcmc_gtb.py → extension/prscs/prscs_mcmc_gtb.py} +1 -1
- gwaslab/g_Sumstats.py +130 -96
- gwaslab/g_SumstatsMulti.py +69 -40
- gwaslab/g_SumstatsPair.py +54 -37
- gwaslab/g_SumstatsSet.py +88 -81
- gwaslab/g_SumstatsT.py +6 -6
- gwaslab/g_Sumstats_polars.py +84 -84
- gwaslab/g_meta_update.py +1 -1
- gwaslab/g_vchange_status.py +4 -4
- gwaslab/g_version.py +2 -2
- gwaslab/{hm_casting.py → hm/hm_casting.py} +4 -4
- gwaslab/{hm_casting_polars.py → hm/hm_casting_polars.py} +4 -4
- gwaslab/hm/hm_harmonize_sumstats.py +1635 -0
- gwaslab/hm_harmonize_sumstats.py +3 -8
- gwaslab/{io_load_ld.py → io/io_load_ld.py} +16 -13
- gwaslab/{io_preformat_input.py → io/io_preformat_input.py} +152 -73
- gwaslab/{io_preformat_input_polars.py → io/io_preformat_input_polars.py} +7 -7
- gwaslab/{io_read_pipcs.py → io/io_read_pipcs.py} +2 -2
- gwaslab/{io_read_tabular.py → io/io_read_tabular.py} +2 -2
- gwaslab/{io_to_formats.py → io/io_to_formats.py} +11 -8
- gwaslab/{io_to_pickle.py → io/io_to_pickle.py} +16 -1
- gwaslab/{qc_check_datatype_polars.py → qc/qc_check_datatype_polars.py} +2 -2
- gwaslab/{qc_fix_sumstats.py → qc/qc_fix_sumstats.py} +60 -33
- gwaslab/{qc_fix_sumstats_polars.py → qc/qc_fix_sumstats_polars.py} +15 -11
- gwaslab/{util_abf_finemapping.py → util/util_abf_finemapping.py} +2 -2
- gwaslab/{util_ex_calculate_ldmatrix.py → util/util_ex_calculate_ldmatrix.py} +18 -8
- gwaslab/{util_ex_calculate_prs.py → util/util_ex_calculate_prs.py} +2 -2
- gwaslab/{util_ex_ldproxyfinder.py → util/util_ex_ldproxyfinder.py} +6 -6
- gwaslab/{util_ex_ldsc.py → util/util_ex_ldsc.py} +18 -13
- gwaslab/{util_ex_match_ldmatrix.py → util/util_ex_match_ldmatrix.py} +8 -7
- gwaslab/util/util_ex_phewwas.py +127 -0
- gwaslab/{util_ex_process_h5.py → util/util_ex_process_h5.py} +2 -2
- gwaslab/{util_ex_process_ref.py → util/util_ex_process_ref.py} +2 -2
- gwaslab/{util_ex_run_2samplemr.py → util/util_ex_run_2samplemr.py} +18 -7
- gwaslab/{util_ex_run_ccgwas.py → util/util_ex_run_ccgwas.py} +4 -4
- gwaslab/{util_ex_run_clumping.py → util/util_ex_run_clumping.py} +28 -13
- gwaslab/{util_ex_run_coloc.py → util/util_ex_run_coloc.py} +22 -10
- gwaslab/{util_ex_run_hyprcoloc.py → util/util_ex_run_hyprcoloc.py} +4 -4
- gwaslab/{util_ex_run_magma.py → util/util_ex_run_magma.py} +21 -11
- gwaslab/{util_ex_run_mesusie.py → util/util_ex_run_mesusie.py} +3 -3
- gwaslab/{util_ex_run_mtag.py → util/util_ex_run_mtag.py} +50 -18
- gwaslab/{util_ex_run_prscs.py → util/util_ex_run_prscs.py} +3 -3
- gwaslab/{util_ex_run_scdrs.py → util/util_ex_run_scdrs.py} +10 -4
- gwaslab/{util_ex_run_susie.py → util/util_ex_run_susie.py} +49 -26
- gwaslab/{util_in_fill_data.py → util/util_in_fill_data.py} +1 -1
- gwaslab/{util_in_filter_value.py → util/util_in_filter_value.py} +18 -11
- gwaslab/{util_in_get_sig.py → util/util_in_get_sig.py} +15 -13
- gwaslab/{util_in_meta.py → util/util_in_meta.py} +1 -1
- gwaslab/{util_in_meta_polars.py → util/util_in_meta_polars.py} +1 -1
- gwaslab/{viz_aux_annotate_plot.py → viz/viz_aux_annotate_plot.py} +1 -1
- gwaslab/{viz_aux_quickfix.py → viz/viz_aux_quickfix.py} +2 -2
- gwaslab/{viz_plot_compare_af.py → viz/viz_plot_compare_af.py} +1 -1
- gwaslab/{viz_plot_compare_effect.py → viz/viz_plot_compare_effect.py} +16 -8
- gwaslab/{viz_plot_credible_sets.py → viz/viz_plot_credible_sets.py} +6 -6
- gwaslab/{viz_plot_effect.py → viz/viz_plot_effect.py} +37 -69
- gwaslab/{viz_plot_miamiplot.py → viz/viz_plot_miamiplot.py} +28 -20
- gwaslab/{viz_plot_miamiplot2.py → viz/viz_plot_miamiplot2.py} +27 -22
- gwaslab/{viz_plot_mqqplot.py → viz/viz_plot_mqqplot.py} +100 -46
- gwaslab/{viz_plot_phe_heatmap.py → viz/viz_plot_phe_heatmap.py} +18 -15
- gwaslab/{viz_plot_qqplot.py → viz/viz_plot_qqplot.py} +12 -28
- gwaslab/{viz_plot_regional2.py → viz/viz_plot_regional2.py} +11 -9
- gwaslab/{viz_plot_regionalplot.py → viz/viz_plot_regionalplot.py} +5 -4
- gwaslab/{viz_plot_rg_heatmap.py → viz/viz_plot_rg_heatmap.py} +1 -1
- gwaslab/{viz_plot_scatter_with_reg.py → viz/viz_plot_scatter_with_reg.py} +10 -7
- gwaslab/{viz_plot_stackedregional.py → viz/viz_plot_stackedregional.py} +67 -33
- gwaslab/{viz_plot_trumpetplot.py → viz/viz_plot_trumpetplot.py} +15 -9
- {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/METADATA +1 -1
- gwaslab-3.6.8.dist-info/RECORD +123 -0
- gwaslab/bd_config.py +0 -18
- gwaslab-3.6.6.dist-info/RECORD +0 -120
- /gwaslab/{ldsc_jackknife.py → extension/ldsc/ldsc_jackknife.py} +0 -0
- /gwaslab/{ldsc_ldscore.py → extension/ldsc/ldsc_ldscore.py} +0 -0
- /gwaslab/{ldsc_parse.py → extension/ldsc/ldsc_parse.py} +0 -0
- /gwaslab/{prscs_gigrnd.py → extension/prscs/prscs_gigrnd.py} +0 -0
- /gwaslab/{prscs_parse_genet.py → extension/prscs/prscs_parse_genet.py} +0 -0
- /gwaslab/{hm_rsid_to_chrpos.py → hm/hm_rsid_to_chrpos.py} +0 -0
- /gwaslab/{io_process_args.py → io/io_process_args.py} +0 -0
- /gwaslab/{io_read_ldsc.py → io/io_read_ldsc.py} +0 -0
- /gwaslab/{qc_build.py → qc/qc_build.py} +0 -0
- /gwaslab/{qc_check_datatype.py → qc/qc_check_datatype.py} +0 -0
- /gwaslab/{util_ex_gwascatalog.py → util/util_ex_gwascatalog.py} +0 -0
- /gwaslab/{util_ex_infer_ancestry.py → util/util_ex_infer_ancestry.py} +0 -0
- /gwaslab/{util_ex_plink_filter.py → util/util_ex_plink_filter.py} +0 -0
- /gwaslab/{util_in_calculate_gc.py → util/util_in_calculate_gc.py} +0 -0
- /gwaslab/{util_in_calculate_power.py → util/util_in_calculate_power.py} +0 -0
- /gwaslab/{util_in_convert_h2.py → util/util_in_convert_h2.py} +0 -0
- /gwaslab/{util_in_correct_winnerscurse.py → util/util_in_correct_winnerscurse.py} +0 -0
- /gwaslab/{util_in_estimate_ess.py → util/util_in_estimate_ess.py} +0 -0
- /gwaslab/{util_in_get_density.py → util/util_in_get_density.py} +0 -0
- /gwaslab/{util_in_merge.py → util/util_in_merge.py} +0 -0
- /gwaslab/{util_in_snphwe.py → util/util_in_snphwe.py} +0 -0
- /gwaslab/{viz_aux_chromatin.py → viz/viz_aux_chromatin.py} +0 -0
- /gwaslab/{viz_aux_property.py → viz/viz_aux_property.py} +0 -0
- /gwaslab/{viz_aux_reposition_text.py → viz/viz_aux_reposition_text.py} +0 -0
- /gwaslab/{viz_aux_save_figure.py → viz/viz_aux_save_figure.py} +0 -0
- /gwaslab/{viz_plot_forestplot.py → viz/viz_plot_forestplot.py} +0 -0
- {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/WHEEL +0 -0
- {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/licenses/LICENSE +0 -0
- {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/licenses/LICENSE_before_v3.4.39 +0 -0
- {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/top_level.txt +0 -0
gwaslab/g_SumstatsSet.py
CHANGED
|
@@ -3,91 +3,98 @@ import numpy as np
|
|
|
3
3
|
import time
|
|
4
4
|
import re
|
|
5
5
|
import copy
|
|
6
|
-
|
|
7
|
-
|
|
8
|
-
from gwaslab.io_preformat_input import preformat
|
|
9
|
-
from gwaslab.io_to_formats import _to_format
|
|
10
|
-
from gwaslab.g_Log import Log
|
|
11
|
-
from gwaslab.qc_fix_sumstats import fixID
|
|
12
|
-
from gwaslab.qc_fix_sumstats import flipSNPID
|
|
13
|
-
from gwaslab.qc_fix_sumstats import stripSNPID
|
|
14
|
-
from gwaslab.qc_fix_sumstats import removedup
|
|
15
|
-
from gwaslab.qc_fix_sumstats import fixchr
|
|
16
|
-
from gwaslab.qc_fix_sumstats import fixpos
|
|
17
|
-
from gwaslab.qc_fix_sumstats import fixallele
|
|
18
|
-
from gwaslab.qc_fix_sumstats import parallelnormalizeallele
|
|
19
|
-
from gwaslab.qc_fix_sumstats import sanitycheckstats
|
|
20
|
-
from gwaslab.qc_fix_sumstats import parallelizeliftovervariant
|
|
21
|
-
from gwaslab.qc_fix_sumstats import flipallelestats
|
|
22
|
-
from gwaslab.qc_fix_sumstats import sortcoordinate
|
|
23
|
-
from gwaslab.qc_fix_sumstats import sortcolumn
|
|
24
|
-
from gwaslab.qc_fix_sumstats import _set_build
|
|
25
|
-
from gwaslab.qc_fix_sumstats import _process_build
|
|
26
|
-
from gwaslab.hm_harmonize_sumstats import parallelecheckaf
|
|
27
|
-
from gwaslab.hm_harmonize_sumstats import paralleleinferaf
|
|
28
|
-
from gwaslab.hm_harmonize_sumstats import checkref
|
|
29
|
-
from gwaslab.hm_harmonize_sumstats import oldcheckref
|
|
30
|
-
from gwaslab.hm_harmonize_sumstats import rsidtochrpos
|
|
31
|
-
from gwaslab.hm_harmonize_sumstats import parallelizeassignrsid
|
|
32
|
-
from gwaslab.hm_harmonize_sumstats import parallelinferstrand
|
|
33
|
-
from gwaslab.hm_harmonize_sumstats import parallelrsidtochrpos
|
|
34
|
-
from gwaslab.hm_harmonize_sumstats import _paralleleinferafwithmaf
|
|
35
|
-
from gwaslab.util_in_filter_value import filtervalues
|
|
36
|
-
from gwaslab.util_in_filter_value import filterout
|
|
37
|
-
from gwaslab.util_in_filter_value import filterin
|
|
38
|
-
from gwaslab.util_in_filter_value import filterregionin
|
|
39
|
-
from gwaslab.util_in_filter_value import filterregionout
|
|
40
|
-
from gwaslab.util_in_filter_value import _filter_indel
|
|
41
|
-
from gwaslab.util_in_filter_value import _filter_palindromic
|
|
42
|
-
from gwaslab.util_in_filter_value import _filter_snp
|
|
43
|
-
from gwaslab.util_in_filter_value import _exclude_hla
|
|
44
|
-
from gwaslab.util_in_filter_value import inferbuild
|
|
45
|
-
from gwaslab.util_in_filter_value import sampling
|
|
46
|
-
from gwaslab.util_in_filter_value import _get_flanking
|
|
47
|
-
from gwaslab.util_in_filter_value import _get_flanking_by_chrpos
|
|
48
|
-
from gwaslab.util_in_filter_value import _get_flanking_by_id
|
|
49
|
-
from gwaslab.util_in_calculate_gc import lambdaGC
|
|
50
|
-
from gwaslab.util_in_convert_h2 import _get_per_snp_r2
|
|
51
|
-
from gwaslab.util_in_get_sig import getsig
|
|
52
|
-
from gwaslab.util_in_get_density import getsignaldensity
|
|
53
|
-
from gwaslab.util_in_get_density import assigndensity
|
|
54
|
-
from gwaslab.util_in_get_sig import annogene
|
|
55
|
-
from gwaslab.util_in_get_sig import getnovel
|
|
56
|
-
from gwaslab.util_in_get_sig import _check_cis
|
|
57
|
-
from gwaslab.util_in_get_sig import _check_novel_set
|
|
58
|
-
from gwaslab.util_in_fill_data import filldata
|
|
59
|
-
from gwaslab.bd_get_hapmap3 import gethapmap3
|
|
60
|
-
from gwaslab.bd_common_data import get_chr_list
|
|
61
|
-
from gwaslab.bd_common_data import get_number_to_chr
|
|
62
|
-
from gwaslab.bd_common_data import get_chr_to_number
|
|
63
|
-
from gwaslab.bd_common_data import get_high_ld
|
|
64
|
-
from gwaslab.bd_common_data import get_format_dict
|
|
65
|
-
from gwaslab.bd_common_data import get_formats_list
|
|
6
|
+
import gc
|
|
7
|
+
|
|
66
8
|
from gwaslab.g_version import _show_version
|
|
67
9
|
from gwaslab.g_version import gwaslab_info
|
|
68
10
|
from gwaslab.g_meta import _init_meta
|
|
69
11
|
from gwaslab.g_meta import _append_meta_record
|
|
70
|
-
from gwaslab.
|
|
71
|
-
from gwaslab.
|
|
72
|
-
from gwaslab.
|
|
73
|
-
|
|
74
|
-
from gwaslab.
|
|
75
|
-
from gwaslab.
|
|
76
|
-
from gwaslab.
|
|
77
|
-
from gwaslab.qc_fix_sumstats import
|
|
78
|
-
from gwaslab.
|
|
79
|
-
from gwaslab.
|
|
80
|
-
from gwaslab.
|
|
81
|
-
from gwaslab.
|
|
82
|
-
from gwaslab.
|
|
83
|
-
from gwaslab.
|
|
84
|
-
from gwaslab.
|
|
85
|
-
from gwaslab.
|
|
86
|
-
from gwaslab.
|
|
87
|
-
import
|
|
88
|
-
from gwaslab.
|
|
89
|
-
from gwaslab.
|
|
90
|
-
|
|
12
|
+
from gwaslab.g_Sumstats_summary import summarize
|
|
13
|
+
from gwaslab.g_Sumstats_summary import lookupstatus
|
|
14
|
+
from gwaslab.g_Log import Log
|
|
15
|
+
|
|
16
|
+
from gwaslab.qc.qc_fix_sumstats import fixID
|
|
17
|
+
from gwaslab.qc.qc_fix_sumstats import flipSNPID
|
|
18
|
+
from gwaslab.qc.qc_fix_sumstats import stripSNPID
|
|
19
|
+
from gwaslab.qc.qc_fix_sumstats import removedup
|
|
20
|
+
from gwaslab.qc.qc_fix_sumstats import fixchr
|
|
21
|
+
from gwaslab.qc.qc_fix_sumstats import fixpos
|
|
22
|
+
from gwaslab.qc.qc_fix_sumstats import fixallele
|
|
23
|
+
from gwaslab.qc.qc_fix_sumstats import parallelnormalizeallele
|
|
24
|
+
from gwaslab.qc.qc_fix_sumstats import sanitycheckstats
|
|
25
|
+
from gwaslab.qc.qc_fix_sumstats import parallelizeliftovervariant
|
|
26
|
+
from gwaslab.qc.qc_fix_sumstats import flipallelestats
|
|
27
|
+
from gwaslab.qc.qc_fix_sumstats import sortcoordinate
|
|
28
|
+
from gwaslab.qc.qc_fix_sumstats import sortcolumn
|
|
29
|
+
from gwaslab.qc.qc_fix_sumstats import _set_build
|
|
30
|
+
from gwaslab.qc.qc_fix_sumstats import _process_build
|
|
31
|
+
from gwaslab.qc.qc_fix_sumstats import _check_data_consistency
|
|
32
|
+
|
|
33
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelecheckaf
|
|
34
|
+
from gwaslab.hm.hm_harmonize_sumstats import paralleleinferaf
|
|
35
|
+
from gwaslab.hm.hm_harmonize_sumstats import checkref
|
|
36
|
+
from gwaslab.hm.hm_harmonize_sumstats import oldcheckref
|
|
37
|
+
from gwaslab.hm.hm_harmonize_sumstats import rsidtochrpos
|
|
38
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelizeassignrsid
|
|
39
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelinferstrand
|
|
40
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelrsidtochrpos
|
|
41
|
+
from gwaslab.hm.hm_harmonize_sumstats import _paralleleinferafwithmaf
|
|
42
|
+
|
|
43
|
+
from gwaslab.util.util_in_filter_value import filtervalues
|
|
44
|
+
from gwaslab.util.util_in_filter_value import filterout
|
|
45
|
+
from gwaslab.util.util_in_filter_value import filterin
|
|
46
|
+
from gwaslab.util.util_in_filter_value import filterregionin
|
|
47
|
+
from gwaslab.util.util_in_filter_value import filterregionout
|
|
48
|
+
from gwaslab.util.util_in_filter_value import _filter_indel
|
|
49
|
+
from gwaslab.util.util_in_filter_value import _filter_palindromic
|
|
50
|
+
from gwaslab.util.util_in_filter_value import _filter_snp
|
|
51
|
+
from gwaslab.util.util_in_filter_value import _exclude_hla
|
|
52
|
+
from gwaslab.util.util_in_filter_value import inferbuild
|
|
53
|
+
from gwaslab.util.util_in_filter_value import sampling
|
|
54
|
+
from gwaslab.util.util_in_filter_value import _get_flanking
|
|
55
|
+
from gwaslab.util.util_in_filter_value import _get_flanking_by_chrpos
|
|
56
|
+
from gwaslab.util.util_in_filter_value import _get_flanking_by_id
|
|
57
|
+
from gwaslab.util.util_in_calculate_gc import lambdaGC
|
|
58
|
+
from gwaslab.util.util_in_convert_h2 import _get_per_snp_r2
|
|
59
|
+
from gwaslab.util.util_in_get_sig import getsig
|
|
60
|
+
from gwaslab.util.util_in_get_density import getsignaldensity
|
|
61
|
+
from gwaslab.util.util_in_get_density import assigndensity
|
|
62
|
+
from gwaslab.util.util_in_get_sig import annogene
|
|
63
|
+
from gwaslab.util.util_in_get_sig import getnovel
|
|
64
|
+
from gwaslab.util.util_in_get_sig import _check_cis
|
|
65
|
+
from gwaslab.util.util_in_get_sig import _check_novel_set
|
|
66
|
+
from gwaslab.util.util_in_fill_data import filldata
|
|
67
|
+
from gwaslab.util.util_ex_run_clumping import _clump
|
|
68
|
+
from gwaslab.util.util_ex_calculate_ldmatrix import tofinemapping
|
|
69
|
+
from gwaslab.util.util_ex_calculate_prs import _calculate_prs
|
|
70
|
+
from gwaslab.util.util_ex_ldsc import _estimate_h2_by_ldsc
|
|
71
|
+
from gwaslab.util.util_ex_ldsc import _estimate_rg_by_ldsc
|
|
72
|
+
from gwaslab.util.util_ex_ldsc import _estimate_h2_cts_by_ldsc
|
|
73
|
+
from gwaslab.util.util_ex_ldsc import _estimate_partitioned_h2_by_ldsc
|
|
74
|
+
from gwaslab.util.util_abf_finemapping import abf_finemapping
|
|
75
|
+
from gwaslab.util.util_abf_finemapping import make_cs
|
|
76
|
+
from gwaslab.util.util_in_merge import _extract_variant
|
|
77
|
+
from gwaslab.util.util_ex_run_susie import _run_susie_rss
|
|
78
|
+
|
|
79
|
+
from gwaslab.bd.bd_get_hapmap3 import gethapmap3
|
|
80
|
+
from gwaslab.bd.bd_common_data import get_chr_list
|
|
81
|
+
from gwaslab.bd.bd_common_data import get_number_to_chr
|
|
82
|
+
from gwaslab.bd.bd_common_data import get_chr_to_number
|
|
83
|
+
from gwaslab.bd.bd_common_data import get_high_ld
|
|
84
|
+
from gwaslab.bd.bd_common_data import get_format_dict
|
|
85
|
+
from gwaslab.bd.bd_common_data import get_formats_list
|
|
86
|
+
from gwaslab.bd.bd_get_hapmap3 import gethapmap3
|
|
87
|
+
|
|
88
|
+
from gwaslab.viz.viz_plot_mqqplot import mqqplot
|
|
89
|
+
from gwaslab.viz.viz_plot_trumpetplot import plottrumpet
|
|
90
|
+
from gwaslab.viz.viz_plot_compare_af import plotdaf
|
|
91
|
+
from gwaslab.viz.viz_plot_credible_sets import _plot_cs
|
|
92
|
+
from gwaslab.viz.viz_plot_phe_heatmap import _gwheatmap
|
|
93
|
+
from gwaslab.viz.viz_plot_effect import _plot_effect
|
|
94
|
+
|
|
95
|
+
from gwaslab.io.io_read_pipcs import _read_pipcs
|
|
96
|
+
from gwaslab.io.io_preformat_input import preformat
|
|
97
|
+
from gwaslab.io.io_to_formats import _to_format
|
|
91
98
|
|
|
92
99
|
#20250215
|
|
93
100
|
class SumstatsSet():
|
gwaslab/g_SumstatsT.py
CHANGED
|
@@ -1,13 +1,13 @@
|
|
|
1
1
|
import pandas as pd
|
|
2
2
|
import numpy as np
|
|
3
3
|
from gwaslab.g_Sumstats import Sumstats
|
|
4
|
-
from gwaslab.hm_casting import _merge_mold_with_sumstats_by_chrpos
|
|
5
|
-
from gwaslab.hm_casting import _align_with_mold
|
|
6
|
-
from gwaslab.hm_casting import _fill_missing_columns
|
|
7
|
-
from gwaslab.hm_casting import _check_daf
|
|
4
|
+
from gwaslab.hm.hm_casting import _merge_mold_with_sumstats_by_chrpos
|
|
5
|
+
from gwaslab.hm.hm_casting import _align_with_mold
|
|
6
|
+
from gwaslab.hm.hm_casting import _fill_missing_columns
|
|
7
|
+
from gwaslab.hm.hm_casting import _check_daf
|
|
8
|
+
from gwaslab.hm.hm_casting import _assign_warning_code
|
|
8
9
|
|
|
9
|
-
from gwaslab.
|
|
10
|
-
from gwaslab.qc_fix_sumstats import flipallelestats
|
|
10
|
+
from gwaslab.qc.qc_fix_sumstats import flipallelestats
|
|
11
11
|
|
|
12
12
|
|
|
13
13
|
class SumstatsT( ):
|
gwaslab/g_Sumstats_polars.py
CHANGED
|
@@ -4,97 +4,97 @@ import time
|
|
|
4
4
|
import copy
|
|
5
5
|
from gwaslab.g_Sumstats_summary import summarize
|
|
6
6
|
from gwaslab.g_Sumstats_summary import lookupstatus
|
|
7
|
-
from gwaslab.io_preformat_input_polars import preformatp
|
|
8
|
-
from gwaslab.io_to_formats import _to_format
|
|
7
|
+
from gwaslab.io.io_preformat_input_polars import preformatp
|
|
8
|
+
from gwaslab.io.io_to_formats import _to_format
|
|
9
9
|
from gwaslab.g_Log import Log
|
|
10
|
-
from gwaslab.qc_fix_sumstats import fixID
|
|
11
|
-
from gwaslab.qc_fix_sumstats import flipSNPID
|
|
12
|
-
from gwaslab.qc_fix_sumstats import stripSNPID
|
|
13
|
-
from gwaslab.qc_fix_sumstats import removedup
|
|
14
|
-
from gwaslab.qc_fix_sumstats import fixchr
|
|
15
|
-
from gwaslab.qc_fix_sumstats import fixpos
|
|
16
|
-
from gwaslab.qc_fix_sumstats import fixallele
|
|
17
|
-
from gwaslab.qc_fix_sumstats import parallelnormalizeallele
|
|
18
|
-
from gwaslab.qc_fix_sumstats import sanitycheckstats
|
|
19
|
-
from gwaslab.qc_fix_sumstats import parallelizeliftovervariant
|
|
20
|
-
from gwaslab.qc_fix_sumstats import flipallelestats
|
|
21
|
-
from gwaslab.qc_fix_sumstats import sortcoordinate
|
|
22
|
-
from gwaslab.qc_fix_sumstats import sortcolumn
|
|
23
|
-
from gwaslab.qc_fix_sumstats import _set_build
|
|
24
|
-
from gwaslab.qc_fix_sumstats import _process_build
|
|
25
|
-
from gwaslab.hm_harmonize_sumstats import parallelecheckaf
|
|
26
|
-
from gwaslab.hm_harmonize_sumstats import paralleleinferaf
|
|
27
|
-
from gwaslab.hm_harmonize_sumstats import checkref
|
|
28
|
-
from gwaslab.hm_harmonize_sumstats import oldcheckref
|
|
29
|
-
from gwaslab.hm_harmonize_sumstats import rsidtochrpos
|
|
30
|
-
from gwaslab.hm_harmonize_sumstats import parallelizeassignrsid
|
|
31
|
-
from gwaslab.hm_harmonize_sumstats import parallelinferstrand
|
|
32
|
-
from gwaslab.hm_harmonize_sumstats import parallelrsidtochrpos
|
|
33
|
-
from gwaslab.hm_harmonize_sumstats import _paralleleinferafwithmaf
|
|
34
|
-
from gwaslab.util_in_filter_value import filtervalues
|
|
35
|
-
from gwaslab.util_in_filter_value import filterout
|
|
36
|
-
from gwaslab.util_in_filter_value import filterin
|
|
37
|
-
from gwaslab.util_in_filter_value import filterregionin
|
|
38
|
-
from gwaslab.util_in_filter_value import filterregionout
|
|
39
|
-
from gwaslab.util_in_filter_value import _filter_indel
|
|
40
|
-
from gwaslab.util_in_filter_value import _filter_palindromic
|
|
41
|
-
from gwaslab.util_in_filter_value import _filter_snp
|
|
42
|
-
from gwaslab.util_in_filter_value import _filter_region
|
|
43
|
-
from gwaslab.util_in_filter_value import _exclude_hla
|
|
44
|
-
from gwaslab.util_in_filter_value import _search_variants
|
|
45
|
-
from gwaslab.util_in_filter_value import inferbuild
|
|
46
|
-
from gwaslab.util_in_filter_value import sampling
|
|
47
|
-
from gwaslab.util_in_filter_value import _get_flanking
|
|
48
|
-
from gwaslab.util_in_filter_value import _get_flanking_by_chrpos
|
|
49
|
-
from gwaslab.util_in_filter_value import _get_flanking_by_id
|
|
50
|
-
from gwaslab.util_in_calculate_gc import lambdaGC
|
|
51
|
-
from gwaslab.util_in_convert_h2 import _get_per_snp_r2
|
|
52
|
-
from gwaslab.util_in_get_sig import getsig
|
|
53
|
-
from gwaslab.util_in_get_density import getsignaldensity
|
|
54
|
-
from gwaslab.util_in_get_density import assigndensity
|
|
55
|
-
from gwaslab.util_in_get_sig import annogene
|
|
56
|
-
from gwaslab.util_in_get_sig import getnovel
|
|
57
|
-
from gwaslab.util_in_get_sig import _check_cis
|
|
58
|
-
from gwaslab.util_in_get_sig import _check_novel_set
|
|
59
|
-
from gwaslab.util_in_fill_data import filldata
|
|
60
|
-
from gwaslab.bd_get_hapmap3 import gethapmap3
|
|
61
|
-
from gwaslab.bd_common_data import get_chr_list
|
|
62
|
-
from gwaslab.bd_common_data import get_number_to_chr
|
|
63
|
-
from gwaslab.bd_common_data import get_chr_to_number
|
|
64
|
-
from gwaslab.bd_common_data import get_high_ld
|
|
65
|
-
from gwaslab.bd_common_data import get_format_dict
|
|
66
|
-
from gwaslab.bd_common_data import get_formats_list
|
|
10
|
+
from gwaslab.qc.qc_fix_sumstats import fixID
|
|
11
|
+
from gwaslab.qc.qc_fix_sumstats import flipSNPID
|
|
12
|
+
from gwaslab.qc.qc_fix_sumstats import stripSNPID
|
|
13
|
+
from gwaslab.qc.qc_fix_sumstats import removedup
|
|
14
|
+
from gwaslab.qc.qc_fix_sumstats import fixchr
|
|
15
|
+
from gwaslab.qc.qc_fix_sumstats import fixpos
|
|
16
|
+
from gwaslab.qc.qc_fix_sumstats import fixallele
|
|
17
|
+
from gwaslab.qc.qc_fix_sumstats import parallelnormalizeallele
|
|
18
|
+
from gwaslab.qc.qc_fix_sumstats import sanitycheckstats
|
|
19
|
+
from gwaslab.qc.qc_fix_sumstats import parallelizeliftovervariant
|
|
20
|
+
from gwaslab.qc.qc_fix_sumstats import flipallelestats
|
|
21
|
+
from gwaslab.qc.qc_fix_sumstats import sortcoordinate
|
|
22
|
+
from gwaslab.qc.qc_fix_sumstats import sortcolumn
|
|
23
|
+
from gwaslab.qc.qc_fix_sumstats import _set_build
|
|
24
|
+
from gwaslab.qc.qc_fix_sumstats import _process_build
|
|
25
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelecheckaf
|
|
26
|
+
from gwaslab.hm.hm_harmonize_sumstats import paralleleinferaf
|
|
27
|
+
from gwaslab.hm.hm_harmonize_sumstats import checkref
|
|
28
|
+
from gwaslab.hm.hm_harmonize_sumstats import oldcheckref
|
|
29
|
+
from gwaslab.hm.hm_harmonize_sumstats import rsidtochrpos
|
|
30
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelizeassignrsid
|
|
31
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelinferstrand
|
|
32
|
+
from gwaslab.hm.hm_harmonize_sumstats import parallelrsidtochrpos
|
|
33
|
+
from gwaslab.hm.hm_harmonize_sumstats import _paralleleinferafwithmaf
|
|
34
|
+
from gwaslab.util.util_in_filter_value import filtervalues
|
|
35
|
+
from gwaslab.util.util_in_filter_value import filterout
|
|
36
|
+
from gwaslab.util.util_in_filter_value import filterin
|
|
37
|
+
from gwaslab.util.util_in_filter_value import filterregionin
|
|
38
|
+
from gwaslab.util.util_in_filter_value import filterregionout
|
|
39
|
+
from gwaslab.util.util_in_filter_value import _filter_indel
|
|
40
|
+
from gwaslab.util.util_in_filter_value import _filter_palindromic
|
|
41
|
+
from gwaslab.util.util_in_filter_value import _filter_snp
|
|
42
|
+
from gwaslab.util.util_in_filter_value import _filter_region
|
|
43
|
+
from gwaslab.util.util_in_filter_value import _exclude_hla
|
|
44
|
+
from gwaslab.util.util_in_filter_value import _search_variants
|
|
45
|
+
from gwaslab.util.util_in_filter_value import inferbuild
|
|
46
|
+
from gwaslab.util.util_in_filter_value import sampling
|
|
47
|
+
from gwaslab.util.util_in_filter_value import _get_flanking
|
|
48
|
+
from gwaslab.util.util_in_filter_value import _get_flanking_by_chrpos
|
|
49
|
+
from gwaslab.util.util_in_filter_value import _get_flanking_by_id
|
|
50
|
+
from gwaslab.util.util_in_calculate_gc import lambdaGC
|
|
51
|
+
from gwaslab.util.util_in_convert_h2 import _get_per_snp_r2
|
|
52
|
+
from gwaslab.util.util_in_get_sig import getsig
|
|
53
|
+
from gwaslab.util.util_in_get_density import getsignaldensity
|
|
54
|
+
from gwaslab.util.util_in_get_density import assigndensity
|
|
55
|
+
from gwaslab.util.util_in_get_sig import annogene
|
|
56
|
+
from gwaslab.util.util_in_get_sig import getnovel
|
|
57
|
+
from gwaslab.util.util_in_get_sig import _check_cis
|
|
58
|
+
from gwaslab.util.util_in_get_sig import _check_novel_set
|
|
59
|
+
from gwaslab.util.util_in_fill_data import filldata
|
|
60
|
+
from gwaslab.bd.bd_get_hapmap3 import gethapmap3
|
|
61
|
+
from gwaslab.bd.bd_common_data import get_chr_list
|
|
62
|
+
from gwaslab.bd.bd_common_data import get_number_to_chr
|
|
63
|
+
from gwaslab.bd.bd_common_data import get_chr_to_number
|
|
64
|
+
from gwaslab.bd.bd_common_data import get_high_ld
|
|
65
|
+
from gwaslab.bd.bd_common_data import get_format_dict
|
|
66
|
+
from gwaslab.bd.bd_common_data import get_formats_list
|
|
67
67
|
from gwaslab.g_version import _show_version
|
|
68
68
|
from gwaslab.g_version import gwaslab_info
|
|
69
69
|
from gwaslab.g_meta import _init_meta
|
|
70
70
|
from gwaslab.g_meta import _append_meta_record
|
|
71
71
|
from gwaslab.g_meta_update import _update_meta
|
|
72
|
-
from gwaslab.util_ex_run_clumping import _clump
|
|
73
|
-
from gwaslab.util_ex_calculate_ldmatrix import tofinemapping
|
|
74
|
-
from gwaslab.io_load_ld import tofinemapping_using_ld
|
|
75
|
-
from gwaslab.util_ex_calculate_prs import _calculate_prs
|
|
76
|
-
from gwaslab.viz_plot_mqqplot import mqqplot
|
|
77
|
-
from gwaslab.viz_plot_trumpetplot import plottrumpet
|
|
78
|
-
from gwaslab.viz_plot_compare_af import plotdaf
|
|
79
|
-
from gwaslab.util_ex_run_susie import _run_susie_rss
|
|
80
|
-
from gwaslab.util_ex_run_susie import _get_cs_lead
|
|
81
|
-
from gwaslab.qc_fix_sumstats import _check_data_consistency
|
|
82
|
-
from gwaslab.util_ex_ldsc import _estimate_h2_by_ldsc
|
|
83
|
-
from gwaslab.util_ex_ldsc import _estimate_rg_by_ldsc
|
|
84
|
-
from gwaslab.util_ex_ldsc import _estimate_h2_cts_by_ldsc
|
|
85
|
-
from gwaslab.util_ex_ldsc import _estimate_partitioned_h2_by_ldsc
|
|
86
|
-
from gwaslab.util_ex_ldproxyfinder import _extract_ld_proxy
|
|
87
|
-
from gwaslab.bd_get_hapmap3 import gethapmap3
|
|
88
|
-
from gwaslab.util_abf_finemapping import abf_finemapping
|
|
89
|
-
from gwaslab.util_abf_finemapping import make_cs
|
|
90
|
-
from gwaslab.io_read_pipcs import _read_pipcs
|
|
91
|
-
from gwaslab.util_in_estimate_ess import _get_ess
|
|
92
|
-
from gwaslab.viz_plot_credible_sets import _plot_cs
|
|
93
|
-
from gwaslab.hm_casting import _align_with_mold
|
|
94
|
-
from gwaslab.hm_casting import _merge_mold_with_sumstats_by_chrpos
|
|
72
|
+
from gwaslab.util.util_ex_run_clumping import _clump
|
|
73
|
+
from gwaslab.util.util_ex_calculate_ldmatrix import tofinemapping
|
|
74
|
+
from gwaslab.io.io_load_ld import tofinemapping_using_ld
|
|
75
|
+
from gwaslab.util.util_ex_calculate_prs import _calculate_prs
|
|
76
|
+
from gwaslab.viz.viz_plot_mqqplot import mqqplot
|
|
77
|
+
from gwaslab.viz.viz_plot_trumpetplot import plottrumpet
|
|
78
|
+
from gwaslab.viz.viz_plot_compare_af import plotdaf
|
|
79
|
+
from gwaslab.util.util_ex_run_susie import _run_susie_rss
|
|
80
|
+
from gwaslab.util.util_ex_run_susie import _get_cs_lead
|
|
81
|
+
from gwaslab.qc.qc_fix_sumstats import _check_data_consistency
|
|
82
|
+
from gwaslab.util.util_ex_ldsc import _estimate_h2_by_ldsc
|
|
83
|
+
from gwaslab.util.util_ex_ldsc import _estimate_rg_by_ldsc
|
|
84
|
+
from gwaslab.util.util_ex_ldsc import _estimate_h2_cts_by_ldsc
|
|
85
|
+
from gwaslab.util.util_ex_ldsc import _estimate_partitioned_h2_by_ldsc
|
|
86
|
+
from gwaslab.util.util_ex_ldproxyfinder import _extract_ld_proxy
|
|
87
|
+
from gwaslab.bd.bd_get_hapmap3 import gethapmap3
|
|
88
|
+
from gwaslab.util.util_abf_finemapping import abf_finemapping
|
|
89
|
+
from gwaslab.util.util_abf_finemapping import make_cs
|
|
90
|
+
from gwaslab.io.io_read_pipcs import _read_pipcs
|
|
91
|
+
from gwaslab.util.util_in_estimate_ess import _get_ess
|
|
92
|
+
from gwaslab.viz.viz_plot_credible_sets import _plot_cs
|
|
93
|
+
from gwaslab.hm.hm_casting import _align_with_mold
|
|
94
|
+
from gwaslab.hm.hm_casting import _merge_mold_with_sumstats_by_chrpos
|
|
95
95
|
import gc
|
|
96
|
-
from gwaslab.viz_plot_phe_heatmap import _gwheatmap
|
|
97
|
-
from gwaslab.util_ex_run_prscs import _run_prscs
|
|
96
|
+
from gwaslab.viz.viz_plot_phe_heatmap import _gwheatmap
|
|
97
|
+
from gwaslab.util.util_ex_run_prscs import _run_prscs
|
|
98
98
|
|
|
99
99
|
#20220309
|
|
100
100
|
class Sumstatsp():
|
gwaslab/g_meta_update.py
CHANGED
gwaslab/g_vchange_status.py
CHANGED
|
@@ -1,15 +1,15 @@
|
|
|
1
1
|
import pandas as pd
|
|
2
2
|
|
|
3
|
-
|
|
3
|
+
STATUS_CATEGORIES = [str(j+i) for j in [1300000,1800000,1900000,3800000,3900000,9700000,9800000,9900000] for i in range(0,100000)]
|
|
4
4
|
|
|
5
5
|
def vchange_status(status,digit,before,after):
|
|
6
6
|
dic={}
|
|
7
7
|
for i in range(len(before)):
|
|
8
|
-
dic[before[i]]=after[i]
|
|
8
|
+
dic[str(before[i])] = str(after[i])
|
|
9
9
|
if digit>1:
|
|
10
|
-
return pd.Categorical(status.str[:digit-1]+status.str[digit-1].replace(dic)+status.str[digit:],categories=
|
|
10
|
+
return pd.Categorical(status.str[:digit-1] + status.str[digit-1].replace(dic) + status.str[digit:],categories=STATUS_CATEGORIES)
|
|
11
11
|
else:
|
|
12
|
-
return pd.Categorical(status.str[digit-1].replace(dic)+status.str[digit:],categories=
|
|
12
|
+
return pd.Categorical(status.str[digit-1].replace(dic) + status.str[digit:],categories=STATUS_CATEGORIES)
|
|
13
13
|
|
|
14
14
|
def copy_status(from_status,to_status, digit):
|
|
15
15
|
if digit>1:
|
gwaslab/g_version.py
CHANGED
|
@@ -4,10 +4,10 @@ from gwaslab.g_Log import Log
|
|
|
4
4
|
from pandas.api.types import CategoricalDtype
|
|
5
5
|
from gwaslab.g_vchange_status import copy_status
|
|
6
6
|
from gwaslab.g_vchange_status import vchange_status
|
|
7
|
-
from gwaslab.qc_fix_sumstats import flipallelestats
|
|
8
|
-
from gwaslab.qc_check_datatype import check_datatype
|
|
9
|
-
from gwaslab.qc_fix_sumstats import start_to
|
|
10
|
-
from gwaslab.util_in_fill_data import filldata
|
|
7
|
+
from gwaslab.qc.qc_fix_sumstats import flipallelestats
|
|
8
|
+
from gwaslab.qc.qc_check_datatype import check_datatype
|
|
9
|
+
from gwaslab.qc.qc_fix_sumstats import start_to
|
|
10
|
+
from gwaslab.util.util_in_fill_data import filldata
|
|
11
11
|
from Bio import SeqIO
|
|
12
12
|
from itertools import combinations
|
|
13
13
|
|
|
@@ -5,10 +5,10 @@ from pandas.api.types import CategoricalDtype
|
|
|
5
5
|
from gwaslab.g_vchange_status import copy_status
|
|
6
6
|
from gwaslab.g_vchange_status_polars import vchange_statusp
|
|
7
7
|
from gwaslab.g_vchange_status_polars import copy_statusp
|
|
8
|
-
from gwaslab.qc_fix_sumstats import flipallelestats
|
|
9
|
-
from gwaslab.qc_check_datatype import check_datatype
|
|
10
|
-
from gwaslab.qc_fix_sumstats import start_to
|
|
11
|
-
from gwaslab.util_in_fill_data import filldata
|
|
8
|
+
from gwaslab.qc.qc_fix_sumstats import flipallelestats
|
|
9
|
+
from gwaslab.qc.qc_check_datatype import check_datatype
|
|
10
|
+
from gwaslab.qc.qc_fix_sumstats import start_to
|
|
11
|
+
from gwaslab.util.util_in_fill_data import filldata
|
|
12
12
|
from Bio import SeqIO
|
|
13
13
|
from itertools import combinations
|
|
14
14
|
import polars as pl
|