pyLoopSage 0.0.6__tar.gz → 0.1.0__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.
@@ -1,6 +1,6 @@
1
1
  Metadata-Version: 2.1
2
2
  Name: pyLoopSage
3
- Version: 0.0.6
3
+ Version: 0.1.0
4
4
  Summary: An energy-based stochastic model of loop extrusion in chromatin.
5
5
  Home-page: https://github.com/SFGLab/pyLoopSage
6
6
  Author: Sebastian Korsak
@@ -29,6 +29,7 @@ Requires-Dist: OpenMM
29
29
  Requires-Dist: OpenMM-cuda
30
30
  Requires-Dist: statsmodels
31
31
  Requires-Dist: imageio
32
+ Requires-Dist: imageio[ffmpeg]
32
33
  Requires-Dist: pillow
33
34
 
34
35
  # pyLoopSage
@@ -50,14 +51,14 @@ Please cite the method paper in case that you would like to use this model for y
50
51
 
51
52
  ## The model
52
53
 
53
- We have a polymer chain with $$N_{\text{beads}}$$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
54
+ We have a polymer chain with $N_{\text{beads}}$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
54
55
 
55
56
  Let's assume that each cohesin $i$ can be represented of two coordinates $(m_{i},n_{i})$ we allow three moves in our simulation:
56
57
 
57
58
  * Slide both locations randomly (as 1D random walk) or
58
59
  * Rebind somewhere else.
59
60
 
60
- In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $$N_{\text{lef}}=2N_{\text{CTCF}}$$.
61
+ In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $N_{\text{lef}}=2N_{\text{CTCF}}$.
61
62
 
62
63
  The main idea of the algorithm is to ensemble loop extrusion from a Boltzmann probability distribution, with Hamiltonian,
63
64
 
@@ -70,7 +71,7 @@ These $L(\cdot), R(\cdot)$ functions are two functions that define the binding p
70
71
 
71
72
  $$c_{\text{fold}}=-\dfrac{N_{\text{beads}}f}{N_{\text{lef}}\log(N_{\text{beads}}/N_{\text{lef}})},\quad c_{\text{bind}}=-\dfrac{N_{\text{beads}}b}{\sum_i \left(L(m_i)+R(n_i)\right)},\quad c_{\text{cross}}=\kappa \times 10^4.$$
72
73
 
73
- The parameters are defined in such a way that when $$f=b=\kappa=1$$, the three terms of the stochastic energy are balanced.
74
+ The parameters are defined in such a way that when $f=b=\kappa=1$, the three terms of the stochastic energy are balanced.
74
75
 
75
76
  And the energy difference can be expressed as the energy difference of each term,
76
77
 
@@ -95,7 +96,7 @@ The main script is `LoopSage.py`. The implementation of the code is very easy an
95
96
 
96
97
  ```python
97
98
  # Definition of Monte Carlo parameters
98
- import Loopsage.stochastic_simulation as lps
99
+ import loopsage.stochastic_simulation as lps
99
100
 
100
101
  N_steps, MC_step, burnin, T, T_min = int(4e4), int(5e2), 1000, 2.5, 1.0
101
102
  mode = 'Metropolis'
@@ -17,14 +17,14 @@ Please cite the method paper in case that you would like to use this model for y
17
17
 
18
18
  ## The model
19
19
 
20
- We have a polymer chain with $$N_{\text{beads}}$$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
20
+ We have a polymer chain with $N_{\text{beads}}$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
21
21
 
22
22
  Let's assume that each cohesin $i$ can be represented of two coordinates $(m_{i},n_{i})$ we allow three moves in our simulation:
23
23
 
24
24
  * Slide both locations randomly (as 1D random walk) or
25
25
  * Rebind somewhere else.
26
26
 
27
- In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $$N_{\text{lef}}=2N_{\text{CTCF}}$$.
27
+ In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $N_{\text{lef}}=2N_{\text{CTCF}}$.
28
28
 
29
29
  The main idea of the algorithm is to ensemble loop extrusion from a Boltzmann probability distribution, with Hamiltonian,
30
30
 
@@ -37,7 +37,7 @@ These $L(\cdot), R(\cdot)$ functions are two functions that define the binding p
37
37
 
38
38
  $$c_{\text{fold}}=-\dfrac{N_{\text{beads}}f}{N_{\text{lef}}\log(N_{\text{beads}}/N_{\text{lef}})},\quad c_{\text{bind}}=-\dfrac{N_{\text{beads}}b}{\sum_i \left(L(m_i)+R(n_i)\right)},\quad c_{\text{cross}}=\kappa \times 10^4.$$
39
39
 
40
- The parameters are defined in such a way that when $$f=b=\kappa=1$$, the three terms of the stochastic energy are balanced.
40
+ The parameters are defined in such a way that when $f=b=\kappa=1$, the three terms of the stochastic energy are balanced.
41
41
 
42
42
  And the energy difference can be expressed as the energy difference of each term,
43
43
 
@@ -62,7 +62,7 @@ The main script is `LoopSage.py`. The implementation of the code is very easy an
62
62
 
63
63
  ```python
64
64
  # Definition of Monte Carlo parameters
65
- import Loopsage.stochastic_simulation as lps
65
+ import loopsage.stochastic_simulation as lps
66
66
 
67
67
  N_steps, MC_step, burnin, T, T_min = int(4e4), int(5e2), 1000, 2.5, 1.0
68
68
  mode = 'Metropolis'
@@ -1,6 +1,6 @@
1
1
  Metadata-Version: 2.1
2
2
  Name: pyLoopSage
3
- Version: 0.0.6
3
+ Version: 0.1.0
4
4
  Summary: An energy-based stochastic model of loop extrusion in chromatin.
5
5
  Home-page: https://github.com/SFGLab/pyLoopSage
6
6
  Author: Sebastian Korsak
@@ -29,6 +29,7 @@ Requires-Dist: OpenMM
29
29
  Requires-Dist: OpenMM-cuda
30
30
  Requires-Dist: statsmodels
31
31
  Requires-Dist: imageio
32
+ Requires-Dist: imageio[ffmpeg]
32
33
  Requires-Dist: pillow
33
34
 
34
35
  # pyLoopSage
@@ -50,14 +51,14 @@ Please cite the method paper in case that you would like to use this model for y
50
51
 
51
52
  ## The model
52
53
 
53
- We have a polymer chain with $$N_{\text{beads}}$$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
54
+ We have a polymer chain with $N_{\text{beads}}$ number of monomers. In general we can scale by deault the granularity of the simulation so as to give reasonable results. Therefore if we have a `region` counted in genomic coordinates, we can assume that `N_beads=(region[1]-region[0])//2000`.
54
55
 
55
56
  Let's assume that each cohesin $i$ can be represented of two coordinates $(m_{i},n_{i})$ we allow three moves in our simulation:
56
57
 
57
58
  * Slide both locations randomly (as 1D random walk) or
58
59
  * Rebind somewhere else.
59
60
 
60
- In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $$N_{\text{lef}}=2N_{\text{CTCF}}$$.
61
+ In general a good working assumption is that the number of cohesins (or loop extrusion factors LEFs) is $N_{\text{lef}}=2N_{\text{CTCF}}$.
61
62
 
62
63
  The main idea of the algorithm is to ensemble loop extrusion from a Boltzmann probability distribution, with Hamiltonian,
63
64
 
@@ -70,7 +71,7 @@ These $L(\cdot), R(\cdot)$ functions are two functions that define the binding p
70
71
 
71
72
  $$c_{\text{fold}}=-\dfrac{N_{\text{beads}}f}{N_{\text{lef}}\log(N_{\text{beads}}/N_{\text{lef}})},\quad c_{\text{bind}}=-\dfrac{N_{\text{beads}}b}{\sum_i \left(L(m_i)+R(n_i)\right)},\quad c_{\text{cross}}=\kappa \times 10^4.$$
72
73
 
73
- The parameters are defined in such a way that when $$f=b=\kappa=1$$, the three terms of the stochastic energy are balanced.
74
+ The parameters are defined in such a way that when $f=b=\kappa=1$, the three terms of the stochastic energy are balanced.
74
75
 
75
76
  And the energy difference can be expressed as the energy difference of each term,
76
77
 
@@ -95,7 +96,7 @@ The main script is `LoopSage.py`. The implementation of the code is very easy an
95
96
 
96
97
  ```python
97
98
  # Definition of Monte Carlo parameters
98
- import Loopsage.stochastic_simulation as lps
99
+ import loopsage.stochastic_simulation as lps
99
100
 
100
101
  N_steps, MC_step, burnin, T, T_min = int(4e4), int(5e2), 1000, 2.5, 1.0
101
102
  mode = 'Metropolis'
@@ -15,4 +15,5 @@ OpenMM
15
15
  OpenMM-cuda
16
16
  statsmodels
17
17
  imageio
18
+ imageio[ffmpeg]
18
19
  pillow
@@ -6,7 +6,7 @@ with open('README.md', 'r', encoding='utf-8') as f:
6
6
 
7
7
  setup(
8
8
  name='pyLoopSage', # Package name
9
- version='0.0.6', # Version of the software
9
+ version='0.1.0', # Version of the software
10
10
  description='An energy-based stochastic model of loop extrusion in chromatin.',
11
11
  long_description=long_description,
12
12
  long_description_content_type='text/markdown',
@@ -33,11 +33,12 @@ setup(
33
33
  'OpenMM-cuda',
34
34
  'statsmodels',
35
35
  'imageio',
36
+ 'imageio[ffmpeg]',
36
37
  'pillow',
37
38
  ],
38
39
  entry_points={
39
40
  'console_scripts': [
40
- 'loopsage=LoopSage.run:main', # MultiMM command points to run.py's main function
41
+ 'loopsage=LoopSage.run:main', # loopsage command points to run.py's main function
41
42
  ],
42
43
  },
43
44
  classifiers=[
File without changes
File without changes