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.

Files changed (110) hide show
  1. gwaslab/__init__.py +57 -47
  2. gwaslab/{bd_common_data.py → bd/bd_common_data.py} +10 -9
  3. gwaslab/bd/bd_config.py +28 -0
  4. gwaslab/{bd_download.py → bd/bd_download.py} +1 -1
  5. gwaslab/{bd_get_hapmap3.py → bd/bd_get_hapmap3.py} +9 -6
  6. gwaslab/bd/bd_path_manager.py +110 -0
  7. gwaslab/data/formatbook.json +805 -9
  8. gwaslab/{ldsc_irwls.py → extension/ldsc/ldsc_irwls.py} +1 -1
  9. gwaslab/{ldsc_regressions.py → extension/ldsc/ldsc_regressions.py} +2 -2
  10. gwaslab/{ldsc_sumstats.py → extension/ldsc/ldsc_sumstats.py} +2 -2
  11. gwaslab/{prscs_mcmc_gtb.py → extension/prscs/prscs_mcmc_gtb.py} +1 -1
  12. gwaslab/g_Sumstats.py +130 -96
  13. gwaslab/g_SumstatsMulti.py +69 -40
  14. gwaslab/g_SumstatsPair.py +54 -37
  15. gwaslab/g_SumstatsSet.py +88 -81
  16. gwaslab/g_SumstatsT.py +6 -6
  17. gwaslab/g_Sumstats_polars.py +84 -84
  18. gwaslab/g_meta_update.py +1 -1
  19. gwaslab/g_vchange_status.py +4 -4
  20. gwaslab/g_version.py +2 -2
  21. gwaslab/{hm_casting.py → hm/hm_casting.py} +4 -4
  22. gwaslab/{hm_casting_polars.py → hm/hm_casting_polars.py} +4 -4
  23. gwaslab/hm/hm_harmonize_sumstats.py +1635 -0
  24. gwaslab/hm_harmonize_sumstats.py +3 -8
  25. gwaslab/{io_load_ld.py → io/io_load_ld.py} +16 -13
  26. gwaslab/{io_preformat_input.py → io/io_preformat_input.py} +152 -73
  27. gwaslab/{io_preformat_input_polars.py → io/io_preformat_input_polars.py} +7 -7
  28. gwaslab/{io_read_pipcs.py → io/io_read_pipcs.py} +2 -2
  29. gwaslab/{io_read_tabular.py → io/io_read_tabular.py} +2 -2
  30. gwaslab/{io_to_formats.py → io/io_to_formats.py} +11 -8
  31. gwaslab/{io_to_pickle.py → io/io_to_pickle.py} +16 -1
  32. gwaslab/{qc_check_datatype_polars.py → qc/qc_check_datatype_polars.py} +2 -2
  33. gwaslab/{qc_fix_sumstats.py → qc/qc_fix_sumstats.py} +60 -33
  34. gwaslab/{qc_fix_sumstats_polars.py → qc/qc_fix_sumstats_polars.py} +15 -11
  35. gwaslab/{util_abf_finemapping.py → util/util_abf_finemapping.py} +2 -2
  36. gwaslab/{util_ex_calculate_ldmatrix.py → util/util_ex_calculate_ldmatrix.py} +18 -8
  37. gwaslab/{util_ex_calculate_prs.py → util/util_ex_calculate_prs.py} +2 -2
  38. gwaslab/{util_ex_ldproxyfinder.py → util/util_ex_ldproxyfinder.py} +6 -6
  39. gwaslab/{util_ex_ldsc.py → util/util_ex_ldsc.py} +18 -13
  40. gwaslab/{util_ex_match_ldmatrix.py → util/util_ex_match_ldmatrix.py} +8 -7
  41. gwaslab/util/util_ex_phewwas.py +127 -0
  42. gwaslab/{util_ex_process_h5.py → util/util_ex_process_h5.py} +2 -2
  43. gwaslab/{util_ex_process_ref.py → util/util_ex_process_ref.py} +2 -2
  44. gwaslab/{util_ex_run_2samplemr.py → util/util_ex_run_2samplemr.py} +18 -7
  45. gwaslab/{util_ex_run_ccgwas.py → util/util_ex_run_ccgwas.py} +4 -4
  46. gwaslab/{util_ex_run_clumping.py → util/util_ex_run_clumping.py} +28 -13
  47. gwaslab/{util_ex_run_coloc.py → util/util_ex_run_coloc.py} +22 -10
  48. gwaslab/{util_ex_run_hyprcoloc.py → util/util_ex_run_hyprcoloc.py} +4 -4
  49. gwaslab/{util_ex_run_magma.py → util/util_ex_run_magma.py} +21 -11
  50. gwaslab/{util_ex_run_mesusie.py → util/util_ex_run_mesusie.py} +3 -3
  51. gwaslab/{util_ex_run_mtag.py → util/util_ex_run_mtag.py} +50 -18
  52. gwaslab/{util_ex_run_prscs.py → util/util_ex_run_prscs.py} +3 -3
  53. gwaslab/{util_ex_run_scdrs.py → util/util_ex_run_scdrs.py} +10 -4
  54. gwaslab/{util_ex_run_susie.py → util/util_ex_run_susie.py} +49 -26
  55. gwaslab/{util_in_fill_data.py → util/util_in_fill_data.py} +1 -1
  56. gwaslab/{util_in_filter_value.py → util/util_in_filter_value.py} +18 -11
  57. gwaslab/{util_in_get_sig.py → util/util_in_get_sig.py} +15 -13
  58. gwaslab/{util_in_meta.py → util/util_in_meta.py} +1 -1
  59. gwaslab/{util_in_meta_polars.py → util/util_in_meta_polars.py} +1 -1
  60. gwaslab/{viz_aux_annotate_plot.py → viz/viz_aux_annotate_plot.py} +1 -1
  61. gwaslab/{viz_aux_quickfix.py → viz/viz_aux_quickfix.py} +2 -2
  62. gwaslab/{viz_plot_compare_af.py → viz/viz_plot_compare_af.py} +1 -1
  63. gwaslab/{viz_plot_compare_effect.py → viz/viz_plot_compare_effect.py} +16 -8
  64. gwaslab/{viz_plot_credible_sets.py → viz/viz_plot_credible_sets.py} +6 -6
  65. gwaslab/{viz_plot_effect.py → viz/viz_plot_effect.py} +37 -69
  66. gwaslab/{viz_plot_miamiplot.py → viz/viz_plot_miamiplot.py} +28 -20
  67. gwaslab/{viz_plot_miamiplot2.py → viz/viz_plot_miamiplot2.py} +27 -22
  68. gwaslab/{viz_plot_mqqplot.py → viz/viz_plot_mqqplot.py} +100 -46
  69. gwaslab/{viz_plot_phe_heatmap.py → viz/viz_plot_phe_heatmap.py} +18 -15
  70. gwaslab/{viz_plot_qqplot.py → viz/viz_plot_qqplot.py} +12 -28
  71. gwaslab/{viz_plot_regional2.py → viz/viz_plot_regional2.py} +11 -9
  72. gwaslab/{viz_plot_regionalplot.py → viz/viz_plot_regionalplot.py} +5 -4
  73. gwaslab/{viz_plot_rg_heatmap.py → viz/viz_plot_rg_heatmap.py} +1 -1
  74. gwaslab/{viz_plot_scatter_with_reg.py → viz/viz_plot_scatter_with_reg.py} +10 -7
  75. gwaslab/{viz_plot_stackedregional.py → viz/viz_plot_stackedregional.py} +67 -33
  76. gwaslab/{viz_plot_trumpetplot.py → viz/viz_plot_trumpetplot.py} +15 -9
  77. {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/METADATA +1 -1
  78. gwaslab-3.6.8.dist-info/RECORD +123 -0
  79. gwaslab/bd_config.py +0 -18
  80. gwaslab-3.6.6.dist-info/RECORD +0 -120
  81. /gwaslab/{ldsc_jackknife.py → extension/ldsc/ldsc_jackknife.py} +0 -0
  82. /gwaslab/{ldsc_ldscore.py → extension/ldsc/ldsc_ldscore.py} +0 -0
  83. /gwaslab/{ldsc_parse.py → extension/ldsc/ldsc_parse.py} +0 -0
  84. /gwaslab/{prscs_gigrnd.py → extension/prscs/prscs_gigrnd.py} +0 -0
  85. /gwaslab/{prscs_parse_genet.py → extension/prscs/prscs_parse_genet.py} +0 -0
  86. /gwaslab/{hm_rsid_to_chrpos.py → hm/hm_rsid_to_chrpos.py} +0 -0
  87. /gwaslab/{io_process_args.py → io/io_process_args.py} +0 -0
  88. /gwaslab/{io_read_ldsc.py → io/io_read_ldsc.py} +0 -0
  89. /gwaslab/{qc_build.py → qc/qc_build.py} +0 -0
  90. /gwaslab/{qc_check_datatype.py → qc/qc_check_datatype.py} +0 -0
  91. /gwaslab/{util_ex_gwascatalog.py → util/util_ex_gwascatalog.py} +0 -0
  92. /gwaslab/{util_ex_infer_ancestry.py → util/util_ex_infer_ancestry.py} +0 -0
  93. /gwaslab/{util_ex_plink_filter.py → util/util_ex_plink_filter.py} +0 -0
  94. /gwaslab/{util_in_calculate_gc.py → util/util_in_calculate_gc.py} +0 -0
  95. /gwaslab/{util_in_calculate_power.py → util/util_in_calculate_power.py} +0 -0
  96. /gwaslab/{util_in_convert_h2.py → util/util_in_convert_h2.py} +0 -0
  97. /gwaslab/{util_in_correct_winnerscurse.py → util/util_in_correct_winnerscurse.py} +0 -0
  98. /gwaslab/{util_in_estimate_ess.py → util/util_in_estimate_ess.py} +0 -0
  99. /gwaslab/{util_in_get_density.py → util/util_in_get_density.py} +0 -0
  100. /gwaslab/{util_in_merge.py → util/util_in_merge.py} +0 -0
  101. /gwaslab/{util_in_snphwe.py → util/util_in_snphwe.py} +0 -0
  102. /gwaslab/{viz_aux_chromatin.py → viz/viz_aux_chromatin.py} +0 -0
  103. /gwaslab/{viz_aux_property.py → viz/viz_aux_property.py} +0 -0
  104. /gwaslab/{viz_aux_reposition_text.py → viz/viz_aux_reposition_text.py} +0 -0
  105. /gwaslab/{viz_aux_save_figure.py → viz/viz_aux_save_figure.py} +0 -0
  106. /gwaslab/{viz_plot_forestplot.py → viz/viz_plot_forestplot.py} +0 -0
  107. {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/WHEEL +0 -0
  108. {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/licenses/LICENSE +0 -0
  109. {gwaslab-3.6.6.dist-info → gwaslab-3.6.8.dist-info}/licenses/LICENSE_before_v3.4.39 +0 -0
  110. {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
- from gwaslab.g_Sumstats_summary import summarize
7
- from gwaslab.g_Sumstats_summary import lookupstatus
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.util_ex_run_clumping import _clump
71
- from gwaslab.util_ex_calculate_ldmatrix import tofinemapping
72
- from gwaslab.util_ex_calculate_prs import _calculate_prs
73
- from gwaslab.viz_plot_mqqplot import mqqplot
74
- from gwaslab.viz_plot_trumpetplot import plottrumpet
75
- from gwaslab.viz_plot_compare_af import plotdaf
76
- from gwaslab.util_ex_run_susie import _run_susie_rss
77
- from gwaslab.qc_fix_sumstats import _check_data_consistency
78
- from gwaslab.util_ex_ldsc import _estimate_h2_by_ldsc
79
- from gwaslab.util_ex_ldsc import _estimate_rg_by_ldsc
80
- from gwaslab.util_ex_ldsc import _estimate_h2_cts_by_ldsc
81
- from gwaslab.util_ex_ldsc import _estimate_partitioned_h2_by_ldsc
82
- from gwaslab.bd_get_hapmap3 import gethapmap3
83
- from gwaslab.util_abf_finemapping import abf_finemapping
84
- from gwaslab.util_abf_finemapping import make_cs
85
- from gwaslab.io_read_pipcs import _read_pipcs
86
- from gwaslab.viz_plot_credible_sets import _plot_cs
87
- import gc
88
- from gwaslab.viz_plot_phe_heatmap import _gwheatmap
89
- from gwaslab.viz_plot_effect import _plot_effect
90
- from gwaslab.util_in_merge import _extract_variant
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.hm_casting import _assign_warning_code
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( ):
@@ -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
@@ -1,5 +1,5 @@
1
1
  import numpy as np
2
- from gwaslab.util_in_filter_value import inferbuild
2
+ from gwaslab.util.util_in_filter_value import inferbuild
3
3
  from gwaslab.g_Log import Log
4
4
  import time
5
5
 
@@ -1,15 +1,15 @@
1
1
  import pandas as pd
2
2
 
3
- CATEGORIES = {str(j+i) for j in [1300000,1800000,1900000,3800000,9700000,9800000,9900000] for i in range(0,100000)}
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=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=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
@@ -17,8 +17,8 @@ def _get_version():
17
17
  def gwaslab_info():
18
18
  # version meta information
19
19
  dic={
20
- "version":"3.6.6",
21
- "release_date":"20250704"
20
+ "version":"3.6.0",
21
+ "release_date":"20250902"
22
22
  }
23
23
  return dic
24
24
 
@@ -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