diff --git a/docs/_images/example_01.png b/docs/_images/example_01.png index 1ee485d..f4ccda1 100644 Binary files a/docs/_images/example_01.png and b/docs/_images/example_01.png differ diff --git a/docs/_images/example_02.png b/docs/_images/example_02.png index 2eb8ccf..22a2bfe 100644 Binary files a/docs/_images/example_02.png and b/docs/_images/example_02.png differ diff --git a/docs/_images/example_06.png b/docs/_images/example_06.png index 5dfcdf5..0ef2243 100644 Binary files a/docs/_images/example_06.png and b/docs/_images/example_06.png differ diff --git a/docs/_images/example_06b.png b/docs/_images/example_06b.png index e86c45e..a87c7e2 100644 Binary files a/docs/_images/example_06b.png and b/docs/_images/example_06b.png differ diff --git a/docs/_images/example_08.png b/docs/_images/example_08.png index cc0d9df..ba46c29 100644 Binary files a/docs/_images/example_08.png and b/docs/_images/example_08.png differ diff --git a/docs/_images/example_13.png b/docs/_images/example_13.png index 918c679..9172d06 100644 Binary files a/docs/_images/example_13.png and b/docs/_images/example_13.png differ diff --git a/docs/_static/documentation_options.js b/docs/_static/documentation_options.js index f29af04..f1bd55b 100644 --- a/docs/_static/documentation_options.js +++ b/docs/_static/documentation_options.js @@ -1,6 +1,6 @@ var DOCUMENTATION_OPTIONS = { URL_ROOT: document.getElementById("documentation_options").getAttribute('data-url_root'), - VERSION: '1.2.9', + VERSION: '1.3.0', LANGUAGE: 'en', COLLAPSE_INDEX: false, BUILDER: 'html', diff --git a/docs/_static/example_01.png b/docs/_static/example_01.png index f4ccda1..2da49f3 100644 Binary files a/docs/_static/example_01.png and b/docs/_static/example_01.png differ diff --git a/docs/_static/example_02.png b/docs/_static/example_02.png index 22a2bfe..cb060dc 100644 Binary files a/docs/_static/example_02.png and b/docs/_static/example_02.png differ diff --git a/docs/_static/example_06.png b/docs/_static/example_06.png index 0ef2243..a4f2592 100644 Binary files a/docs/_static/example_06.png and b/docs/_static/example_06.png differ diff --git a/docs/_static/example_06b.png b/docs/_static/example_06b.png index a87c7e2..129da2a 100644 Binary files a/docs/_static/example_06b.png and b/docs/_static/example_06b.png differ diff --git a/docs/_static/example_08.png b/docs/_static/example_08.png index ba46c29..e44412b 100644 Binary files a/docs/_static/example_08.png and b/docs/_static/example_08.png differ diff --git a/docs/_static/example_13.png b/docs/_static/example_13.png index 9172d06..3397a81 100644 Binary files a/docs/_static/example_13.png and b/docs/_static/example_13.png differ diff --git a/docs/_static/example_pa_01.pdf b/docs/_static/example_pa_01.pdf index db3b912..272d7fb 100644 Binary files a/docs/_static/example_pa_01.pdf and b/docs/_static/example_pa_01.pdf differ diff --git a/docs/_static/example_pa_02.pdf b/docs/_static/example_pa_02.pdf index 8b8e868..ca579c3 100644 Binary files a/docs/_static/example_pa_02.pdf and b/docs/_static/example_pa_02.pdf differ diff --git a/docs/_static/example_pa_03.pdf b/docs/_static/example_pa_03.pdf index 9303b3a..f7f16db 100644 Binary files a/docs/_static/example_pa_03.pdf and b/docs/_static/example_pa_03.pdf differ diff --git a/docs/_static/example_pa_04.pdf b/docs/_static/example_pa_04.pdf index 7b501d3..efab0c6 100644 Binary files a/docs/_static/example_pa_04.pdf and b/docs/_static/example_pa_04.pdf differ diff --git a/docs/_static/merge_ologram_stats_01.pdf b/docs/_static/merge_ologram_stats_01.pdf index 86454df..3291f80 100644 Binary files a/docs/_static/merge_ologram_stats_01.pdf and b/docs/_static/merge_ologram_stats_01.pdf differ diff --git a/docs/_static/treeified.pdf b/docs/_static/treeified.pdf index e69de29..1478791 100644 Binary files a/docs/_static/treeified.pdf and b/docs/_static/treeified.pdf differ diff --git a/docs/about.html b/docs/about.html index 4ad1a0e..174791d 100644 --- a/docs/about.html +++ b/docs/about.html @@ -17,7 +17,7 @@ var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s); })(); -
Print the nb last lines of the GTF object.
nb – The number of line to display.
+nb – The number of line to display.
returned – If True, don’t print but returns the message.
>>> from pygtftk.utils import get_example_file
@@ -3681,7 +3684,7 @@ Navigation
$ gtftk profile -D -i mini_real_promoter.zip -o profile_prom -pf png -if example_01.png
- |-- 11:02-WARNING-profile : --group-by not set. Choosing 'bwig'.
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:21-WARNING-profile : --group-by not set. Choosing 'bwig'.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
Changing colors and applying color order can be done using the following syntax:
$ gtftk profile -D -i mini_real_promoter.zip -c 'red,blue,violet' -d H3K79me,H3K4me3,H3K36me3 -o profile_prom -pf png -if example_01b.png
- |-- 11:02-WARNING-profile : --group-by not set. Choosing 'bwig'.
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:21-WARNING-profile : --group-by not set. Choosing 'bwig'.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
Transcript coverage is obtained using the mini_real_tx.zip matrix. This provides a simple overlayed profile of all epigenetic marks along the transcript body extended in 5’ and 3’ regions:
$ gtftk profile -D -i mini_real_tx.zip -o profile_tx -pf png -if example_02.png
- |-- 11:02-WARNING-profile : --group-by not set. Choosing 'bwig'.
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:21-WARNING-profile : --group-by not set. Choosing 'bwig'.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
$ gtftk profile -D -i mini_real_promoter.zip -f tx_classes -g bwig -t tx_classes.txt -o profile_prom -pf png -if example_05.png -e -V 2 -fc 2
- |-- 11:02-DEBUG-profile : Using pandas version 1.2.4
- |-- 11:02-DEBUG-profile : Pandas location /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pandas/__init__.py
- |-- 11:02-DEBUG-profile : Using numpy version 1.20.3
- |-- 11:02-DEBUG-profile : Pandas numpy /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/numpy/__init__.py
- |-- 11:02-DEBUG-profile : Using plotnine version 0.8.0
- |-- 11:02-DEBUG-profile : Pandas plotnine /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/plotnine/__init__.py
- |-- 11:02-DEBUG-profile : Creating directory : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_cxvl97ol
- |-- 11:02-DEBUG-profile : Uncompressing : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_cxvl97ol
- |-- 11:02-DEBUG-profile : Reading : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_cxvl97ol/mini_real_promoter
- |-- 11:02-INFO-profile : Getting configuration info from input file.
- |-- 11:02-INFO-profile : Reading transcript file.
- |-- 11:02-INFO-profile : Deleting duplicates in transcript-file.
- |-- 11:02-INFO-profile : Checking how many transcripts where found in the transcript list.
- |-- 11:02-INFO-profile : Keeping 804 transcript out of 833 in input transcript list.
- |-- 11:02-DEBUG-profile : Color order : ['H3K36me3', 'H3K79me', 'H3K4me3']
- |-- 11:02-DEBUG-profile : Profile color : ['#000000', '#00bb00', '#cccccc']
- |-- 11:02-INFO-profile : Searching coverage columns.
- |-- 11:02-INFO-profile : Melting.
- |-- 11:02-INFO-profile : Ceiling
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
- |-- 11:02-INFO-profile : Computing column ordering.
- |-- 11:02-INFO-profile : Preparing diagram
- |-- 11:02-INFO-profile : Theming and ordering. Please be patient...
- |-- 11:02-INFO-profile : Preparing x axis
- |-- 11:02-INFO-profile : facet_col 2
- |-- 11:02-INFO-profile : Page width set to 6
- |-- 11:02-INFO-profile : Page height set to 5.0
- |-- 11:02-INFO-profile : Saving diagram to file : example_05.png
- |-- 11:02-INFO-profile : Be patient. This may be long for large datasets.
- |-- 11:02-DEBUG-profile : Deleting temp file : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_cxvl97ol
+ |-- 13:21-DEBUG-profile : Using pandas version 1.2.4
+ |-- 13:21-DEBUG-profile : Pandas location /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pandas/__init__.py
+ |-- 13:21-DEBUG-profile : Using numpy version 1.20.3
+ |-- 13:21-DEBUG-profile : Pandas numpy /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/numpy/__init__.py
+ |-- 13:21-DEBUG-profile : Using plotnine version 0.8.0
+ |-- 13:21-DEBUG-profile : Pandas plotnine /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/plotnine/__init__.py
+ |-- 13:21-DEBUG-profile : Creating directory : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_jqjsd6n_
+ |-- 13:21-DEBUG-profile : Uncompressing : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_jqjsd6n_
+ |-- 13:21-DEBUG-profile : Reading : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_jqjsd6n_/mini_real_promoter
+ |-- 13:21-INFO-profile : Getting configuration info from input file.
+ |-- 13:21-INFO-profile : Reading transcript file.
+ |-- 13:21-INFO-profile : Deleting duplicates in transcript-file.
+ |-- 13:21-INFO-profile : Checking how many transcripts where found in the transcript list.
+ |-- 13:21-INFO-profile : Keeping 804 transcript out of 833 in input transcript list.
+ |-- 13:21-DEBUG-profile : Color order : ['H3K36me3', 'H3K79me', 'H3K4me3']
+ |-- 13:21-DEBUG-profile : Profile color : ['#000000', '#00bb00', '#cccccc']
+ |-- 13:21-INFO-profile : Searching coverage columns.
+ |-- 13:21-INFO-profile : Melting.
+ |-- 13:21-INFO-profile : Ceiling
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:21-INFO-profile : Computing column ordering.
+ |-- 13:21-INFO-profile : Preparing diagram
+ |-- 13:21-INFO-profile : Theming and ordering. Please be patient...
+ |-- 13:21-INFO-profile : Preparing x axis
+ |-- 13:21-INFO-profile : facet_col 2
+ |-- 13:21-INFO-profile : Page width set to 6
+ |-- 13:21-INFO-profile : Page height set to 5.0
+ |-- 13:21-INFO-profile : Saving diagram to file : example_05.png
+ |-- 13:21-INFO-profile : Be patient. This may be long for large datasets.
+ |-- 13:21-DEBUG-profile : Deleting temp file : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_jqjsd6n_
Alternatively, the groups can be set to chromosomes or transcript classes:
$ gtftk profile -D -i mini_real_promoter.zip -g tx_classes -f bwig -t tx_classes.txt -o profile_prom -pf png -if example_06.png -V 2 -nm ranging
- |-- 11:02-DEBUG-profile : Using pandas version 1.2.4
- |-- 11:02-DEBUG-profile : Pandas location /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pandas/__init__.py
- |-- 11:02-DEBUG-profile : Using numpy version 1.20.3
- |-- 11:02-DEBUG-profile : Pandas numpy /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/numpy/__init__.py
- |-- 11:02-DEBUG-profile : Using plotnine version 0.8.0
- |-- 11:02-DEBUG-profile : Pandas plotnine /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/plotnine/__init__.py
- |-- 11:02-DEBUG-profile : Creating directory : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_f198it4c
- |-- 11:02-DEBUG-profile : Uncompressing : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_f198it4c
- |-- 11:02-DEBUG-profile : Reading : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_f198it4c/mini_real_promoter
- |-- 11:02-INFO-profile : Getting configuration info from input file.
- |-- 11:02-INFO-profile : Reading transcript file.
- |-- 11:02-INFO-profile : Deleting duplicates in transcript-file.
- |-- 11:02-INFO-profile : Checking how many transcripts where found in the transcript list.
- |-- 11:02-INFO-profile : Keeping 804 transcript out of 833 in input transcript list.
- |-- 11:02-DEBUG-profile : Color order : ['lincRNA', 'antisense', 'protein_coding']
- |-- 11:02-DEBUG-profile : Profile color : ['#000000', '#00bb00', '#cccccc']
- |-- 11:02-INFO-profile : Searching coverage columns.
- |-- 11:02-INFO-profile : Melting.
- |-- 11:02-INFO-profile : Ceiling
- |-- 11:02-INFO-profile : Normalizing (ranging)
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
- |-- 11:02-INFO-profile : Computing column ordering.
- |-- 11:02-INFO-profile : Preparing diagram
- |-- 11:02-INFO-profile : Theming and ordering. Please be patient...
- |-- 11:02-INFO-profile : Preparing x axis
- |-- 11:02-INFO-profile : facet_col 3
- |-- 11:02-INFO-profile : Page width set to 9
- |-- 11:02-INFO-profile : Page height set to 2.0
- |-- 11:02-INFO-profile : Saving diagram to file : example_06.png
- |-- 11:02-INFO-profile : Be patient. This may be long for large datasets.
- |-- 11:02-DEBUG-profile : Deleting temp file : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_f198it4c
+ |-- 13:21-DEBUG-profile : Using pandas version 1.2.4
+ |-- 13:21-DEBUG-profile : Pandas location /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pandas/__init__.py
+ |-- 13:21-DEBUG-profile : Using numpy version 1.20.3
+ |-- 13:21-DEBUG-profile : Pandas numpy /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/numpy/__init__.py
+ |-- 13:21-DEBUG-profile : Using plotnine version 0.8.0
+ |-- 13:21-DEBUG-profile : Pandas plotnine /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/plotnine/__init__.py
+ |-- 13:21-DEBUG-profile : Creating directory : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_xb7sgm3e
+ |-- 13:21-DEBUG-profile : Uncompressing : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_xb7sgm3e
+ |-- 13:21-DEBUG-profile : Reading : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_xb7sgm3e/mini_real_promoter
+ |-- 13:21-INFO-profile : Getting configuration info from input file.
+ |-- 13:21-INFO-profile : Reading transcript file.
+ |-- 13:21-INFO-profile : Deleting duplicates in transcript-file.
+ |-- 13:21-INFO-profile : Checking how many transcripts where found in the transcript list.
+ |-- 13:21-INFO-profile : Keeping 804 transcript out of 833 in input transcript list.
+ |-- 13:21-DEBUG-profile : Color order : ['lincRNA', 'protein_coding', 'antisense']
+ |-- 13:21-DEBUG-profile : Profile color : ['#000000', '#00bb00', '#cccccc']
+ |-- 13:21-INFO-profile : Searching coverage columns.
+ |-- 13:21-INFO-profile : Melting.
+ |-- 13:21-INFO-profile : Ceiling
+ |-- 13:21-INFO-profile : Normalizing (ranging)
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:22-INFO-profile : Computing column ordering.
+ |-- 13:22-INFO-profile : Preparing diagram
+ |-- 13:22-INFO-profile : Theming and ordering. Please be patient...
+ |-- 13:22-INFO-profile : Preparing x axis
+ |-- 13:22-INFO-profile : facet_col 3
+ |-- 13:22-INFO-profile : Page width set to 9
+ |-- 13:22-INFO-profile : Page height set to 2.0
+ |-- 13:22-INFO-profile : Saving diagram to file : example_06.png
+ |-- 13:22-INFO-profile : Be patient. This may be long for large datasets.
+ |-- 13:22-DEBUG-profile : Deleting temp file : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_xb7sgm3e
$ gtftk profile -D -i mini_real_promoter.zip -g chrom -f bwig -t tx_classes.txt -o profile_prom -pf png -if example_06b.png -V 2 -nm ranging
- |-- 11:02-DEBUG-profile : Using pandas version 1.2.4
- |-- 11:02-DEBUG-profile : Pandas location /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pandas/__init__.py
- |-- 11:02-DEBUG-profile : Using numpy version 1.20.3
- |-- 11:02-DEBUG-profile : Pandas numpy /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/numpy/__init__.py
- |-- 11:02-DEBUG-profile : Using plotnine version 0.8.0
- |-- 11:02-DEBUG-profile : Pandas plotnine /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/plotnine/__init__.py
- |-- 11:02-DEBUG-profile : Creating directory : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_woyx5gcz
- |-- 11:02-DEBUG-profile : Uncompressing : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_woyx5gcz
- |-- 11:02-DEBUG-profile : Reading : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_woyx5gcz/mini_real_promoter
- |-- 11:02-INFO-profile : Getting configuration info from input file.
- |-- 11:02-DEBUG-profile : Color order : ['chr1', 'chr6', 'chr15', 'chr7', 'chr9', 'chr5', 'chr8', 'chr14', 'chr11', 'chr19', 'chr3', 'chr4', 'chr10', 'chr16', 'chr2', 'chrX', 'chr17', 'chr21', 'chr12', 'chr18', 'chr22', 'chr13', 'chr20']
- |-- 11:02-DEBUG-profile : Profile color : ['#000000', '#6c007c', '#850096', '#2500a5', '#0000ca', '#0041dd', '#0086dd', '#009fca', '#00aaa1', '#00a76f', '#009c00', '#00bb00', '#00da00', '#00f900', '#88ff00', '#dbf400', '#f7db00', '#ffb500', '#ff6100', '#f60000', '#da0000', '#cc1313', '#cccccc']
- |-- 11:02-INFO-profile : Searching coverage columns.
- |-- 11:02-INFO-profile : Melting.
- |-- 11:02-INFO-profile : Ceiling
- |-- 11:02-INFO-profile : Normalizing (ranging)
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
- |-- 11:02-INFO-profile : Computing column ordering.
- |-- 11:02-INFO-profile : Preparing diagram
- |-- 11:02-INFO-profile : Theming and ordering. Please be patient...
- |-- 11:02-INFO-profile : Preparing x axis
- |-- 11:02-INFO-profile : facet_col 3
- |-- 11:02-INFO-profile : Page width set to 9
- |-- 11:02-INFO-profile : Page height set to 2.0
- |-- 11:02-INFO-profile : Saving diagram to file : example_06b.png
- |-- 11:02-INFO-profile : Be patient. This may be long for large datasets.
- |-- 11:02-DEBUG-profile : Deleting temp file : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_woyx5gcz
+ |-- 13:22-DEBUG-profile : Using pandas version 1.2.4
+ |-- 13:22-DEBUG-profile : Pandas location /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pandas/__init__.py
+ |-- 13:22-DEBUG-profile : Using numpy version 1.20.3
+ |-- 13:22-DEBUG-profile : Pandas numpy /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/numpy/__init__.py
+ |-- 13:22-DEBUG-profile : Using plotnine version 0.8.0
+ |-- 13:22-DEBUG-profile : Pandas plotnine /Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/plotnine/__init__.py
+ |-- 13:22-DEBUG-profile : Creating directory : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_umru5mcw
+ |-- 13:22-DEBUG-profile : Uncompressing : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_umru5mcw
+ |-- 13:22-DEBUG-profile : Reading : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_umru5mcw/mini_real_promoter
+ |-- 13:22-INFO-profile : Getting configuration info from input file.
+ |-- 13:22-DEBUG-profile : Color order : ['chr9', 'chr1', 'chr22', 'chr6', 'chr13', 'chrX', 'chr21', 'chr2', 'chr12', 'chr5', 'chr18', 'chr14', 'chr8', 'chr10', 'chr17', 'chr15', 'chr20', 'chr16', 'chr4', 'chr7', 'chr11', 'chr19', 'chr3']
+ |-- 13:22-DEBUG-profile : Profile color : ['#000000', '#6c007c', '#850096', '#2500a5', '#0000ca', '#0041dd', '#0086dd', '#009fca', '#00aaa1', '#00a76f', '#009c00', '#00bb00', '#00da00', '#00f900', '#88ff00', '#dbf400', '#f7db00', '#ffb500', '#ff6100', '#f60000', '#da0000', '#cc1313', '#cccccc']
+ |-- 13:22-INFO-profile : Searching coverage columns.
+ |-- 13:22-INFO-profile : Melting.
+ |-- 13:22-INFO-profile : Ceiling
+ |-- 13:22-INFO-profile : Normalizing (ranging)
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:22-INFO-profile : Computing column ordering.
+ |-- 13:22-INFO-profile : Preparing diagram
+ |-- 13:22-INFO-profile : Theming and ordering. Please be patient...
+ |-- 13:22-INFO-profile : Preparing x axis
+ |-- 13:22-INFO-profile : facet_col 3
+ |-- 13:22-INFO-profile : Page width set to 9
+ |-- 13:22-INFO-profile : Page height set to 2.0
+ |-- 13:22-INFO-profile : Saving diagram to file : example_06b.png
+ |-- 13:22-INFO-profile : Be patient. This may be long for large datasets.
+ |-- 13:22-DEBUG-profile : Deleting temp file : /var/folders/zy/wl3dj2_n76zfc8sdvny1q06c0000gn/T/profile_matrix__pygtftk_umru5mcw
Note that facets may also be associated to epigenetic marks. In this case each the –group-by can be set to tx_classes or chrom.
$ gtftk profile -D -i mini_real_tx.zip -g tx_classes -t tx_classes.txt -f bwig -o profile_tx -pf png -if example_07.png -w -nm ranging
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:890: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
- |-- 11:02-WARNING-profile : PlotnineError error, can not show group number: Aesthetics {'ha'} specified two times.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:890: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:22-WARNING-profile : PlotnineError error, can not show group number: Aesthetics {'ha'} specified two times.
$ gtftk profile -D -i mini_real_tx.zip -g chrom -f bwig -o profile_tx -pf png -if example_08.png -w -nm ranging
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
-/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.2.9-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:890: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
- |-- 11:03-WARNING-profile : PlotnineError error, can not show group number: Aesthetics {'ha'} specified two times.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:815: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+/Users/puthier/anaconda3/envs/python_dev_388/lib/python3.8/site-packages/pygtftk-1.3.0-py3.8-macosx-10.9-x86_64.egg/pygtftk/plugins/profile.py:890: FutureWarning: Index.ravel returning ndarray is deprecated; in a future version this will return a view on self.
+ |-- 13:22-WARNING-profile : PlotnineError error, can not show group number: Aesthetics {'ha'} specified two times.
$ gtftk join_attr -i simple.gtf -j simple.join_mat -k gene_id -m | gtftk discretize_key -k S1 -d S1_d -n 2 -l A,B | gtftk select_by_key -k feature -v gene
- |-- 11:03-INFO-discretize_key : Categories: ['A', 'B']
+ |-- 13:23-INFO-discretize_key : Categories: ['A', 'B']
chr1 gtftk gene 125 138 . + . gene_id "G0001";
chr1 gtftk gene 180 189 . + . gene_id "G0002";
chr1 gtftk gene 50 61 . - . gene_id "G0003"; S1 "0.2322"; S2 "0.4"; S1_d "A";
@@ -450,7 +450,7 @@ discretize_key$ gtftk join_attr -i mini_real.gtf.gz -H -j mini_real_counts_ENCFF630HEX.tsv -k gene_name -n exprs -t gene | gtftk discretize_key -k exprs -p -d exprs_class -n 10 -l A,B,C,D,E,F,G,H,I,J | gtftk tabulate -k exprs_class -Hn | sort | uniq -c
- |-- 11:03-INFO-discretize_key : Categories: ['A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J']
+ |-- 13:23-INFO-discretize_key : Categories: ['A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J']
96 A
83 B
89 C
@@ -567,7 +567,7 @@ Navigation
previous |
- gtftk 1.2.9 documentation »
+ gtftk 1.3.0 documentation »
Commands from section ‘Editing’
diff --git a/docs/example_01.png b/docs/example_01.png
index 1ee485d..f4ccda1 100644
Binary files a/docs/example_01.png and b/docs/example_01.png differ
diff --git a/docs/example_02.png b/docs/example_02.png
index 2eb8ccf..22a2bfe 100644
Binary files a/docs/example_02.png and b/docs/example_02.png differ
diff --git a/docs/example_06.png b/docs/example_06.png
index 5dfcdf5..0ef2243 100644
Binary files a/docs/example_06.png and b/docs/example_06.png differ
diff --git a/docs/example_06b.png b/docs/example_06b.png
index e86c45e..a87c7e2 100644
Binary files a/docs/example_06b.png and b/docs/example_06b.png differ
diff --git a/docs/example_08.png b/docs/example_08.png
index cc0d9df..ba46c29 100644
Binary files a/docs/example_08.png and b/docs/example_08.png differ
diff --git a/docs/example_13.png b/docs/example_13.png
index 918c679..9172d06 100644
Binary files a/docs/example_13.png and b/docs/example_13.png differ
diff --git a/docs/example_pa_01.pdf b/docs/example_pa_01.pdf
index 9c9c33e..db3b912 100644
Binary files a/docs/example_pa_01.pdf and b/docs/example_pa_01.pdf differ
diff --git a/docs/example_pa_02.pdf b/docs/example_pa_02.pdf
index 8226f96..8b8e868 100644
Binary files a/docs/example_pa_02.pdf and b/docs/example_pa_02.pdf differ
diff --git a/docs/example_pa_03.pdf b/docs/example_pa_03.pdf
index a8b8906..9303b3a 100644
Binary files a/docs/example_pa_03.pdf and b/docs/example_pa_03.pdf differ
diff --git a/docs/example_pa_04.pdf b/docs/example_pa_04.pdf
index 0d33b73..7b501d3 100644
Binary files a/docs/example_pa_04.pdf and b/docs/example_pa_04.pdf differ
diff --git a/docs/genindex.html b/docs/genindex.html
index 3ccb722..fe19693 100644
--- a/docs/genindex.html
+++ b/docs/genindex.html
@@ -16,7 +16,7 @@
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s);
})();
- Index — gtftk 1.2.9 documentation
+ Index — gtftk 1.3.0 documentation
@@ -36,7 +36,7 @@ Navigation
modules |
- gtftk 1.2.9 documentation »
+ gtftk 1.3.0 documentation »
Index
$ gtftk apropos -k promoter
- |-- 11:03-INFO-apropos : >> Keyword 'promoter' was found in the following command:
+ |-- 13:23-INFO-apropos : >> Keyword 'promoter' was found in the following command:
+ - divergent.
- ologram.
- coverage.
- - divergent.
Arguments:
@@ -724,7 +724,7 @@$ gtftk control_list -i mini_real_counts_ENCFF630HEX.tsv -r mini_real_control_1.txt -D -V 2 -s -l -p 1 -ju -if example_13.png -pf png
- |-- 11:04-INFO-control_list : 0 duplicate lines have been deleted in reference file.
- |-- 11:04-INFO-control_list : Found 50 genes of the reference in the provided signal file
- |-- 11:04-INFO-control_list : All reference genes were found.
- |-- 11:04-INFO-control_list : Searching for genes with matched signal.
- |-- 11:04-INFO-control_list : Preparing a dataframe for plotting.
- |-- 11:04-INFO-control_list : Saving diagram to file : example_13.png
- |-- 11:04-INFO-control_list : Be patient. This may be long for large datasets.
+ |-- 13:23-INFO-control_list : 0 duplicate lines have been deleted in reference file.
+ |-- 13:23-INFO-control_list : Found 50 genes of the reference in the provided signal file
+ |-- 13:23-INFO-control_list : All reference genes were found.
+ |-- 13:23-INFO-control_list : Searching for genes with matched signal.
+ |-- 13:23-INFO-control_list : Preparing a dataframe for plotting.
+ |-- 13:23-INFO-control_list : Saving diagram to file : example_13.png
+ |-- 13:23-INFO-control_list : Be patient. This may be long for large datasets.
$ gtftk ologram -i hg38_chr1.gtf.gz -p ENCFF112BHN_H3K4me3_chr1.bed -c hg38_chr1.genome -u 1500 -d 1500 -D -pf example_pa_01.pdf -k 8 -j summed_bp_overlaps_pvalue
- |-- 11:04-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
- |-- 11:05-WARNING-ologram : Computing log(p-val) for a Neg Binom with mean >= var ; var was set to mean+1 (start_codon)
+ |-- 13:23-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
+ |-- 13:24-WARNING-ologram : Computing log(p-val) for a Neg Binom with mean >= var ; var was set to mean+1 (start_codon)
$ gtftk select_by_key -i mini_real.gtf.gz -k gene_biotype -v protein_coding,lincRNA,antisense,processed_transcript | gtftk ologram -m gene_biotype -p ENCFF112BHN_H3K4me3_K562_sub.bed -c hg38 -D -n -pf example_pa_02.pdf -k 8 -j summed_bp_overlaps_pvalue
- |-- 11:06-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
+ |-- 13:24-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
$ gtftk join_attr -i mini_real.gtf.gz -H -j mini_real_counts_ENCFF630HEX.tsv -k gene_name -n exprs -t exon | gtftk discretize_key -k exprs -p -d exprs_class -n 6 -u | gtftk ologram -p ENCFF119BYM_H3K36me3_K562_sub.bed -c hg38 -D -n -m exprs_class -pf example_pa_03.pdf -k 8 -j summed_bp_overlaps_pvalue
- |-- 11:06-INFO-discretize_key : Categories: ['[0.0_183.0]', '(183.0_549.0]', '(549.0_1018.0]', '(1018.0_1631.0]', '(1631.0_3139.0]', '(3139.0_41703.0]']
- |-- 11:06-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
+ |-- 13:25-INFO-discretize_key : Categories: ['[0.0_183.0]', '(183.0_549.0]', '(549.0_1018.0]', '(1018.0_1631.0]', '(1631.0_3139.0]', '(3139.0_41703.0]']
+ |-- 13:25-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
$ gtftk add_exon_nb -k exon_nbr -i mini_real.gtf.gz | gtftk discretize_key -p -d exon_nbr_cat -n 5 -k exon_nbr | gtftk ologram -p ENCFF112BHN_H3K4me3_K562_sub.bed -c hg38 -D -n -m exon_nbr_cat -pf example_pa_04.pdf -k 8 -j summed_bp_overlaps_pvalue
- |-- 11:07-INFO-discretize_key : Categories: ['[1.0_2.0]', '(2.0_4.0]', '(4.0_6.0]', '(6.0_12.0]', '(12.0_107.0]']
- |-- 11:07-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
+ |-- 13:26-INFO-discretize_key : Categories: ['[1.0_2.0]', '(2.0_4.0]', '(4.0_6.0]', '(6.0_12.0]', '(12.0_107.0]']
+ |-- 13:26-WARNING-ologram : Using only 8 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
$ gtftk ologram -z -w -q -c simple_07.chromInfo -p simple_07_peaks.bed --more-bed simple_07_peaks.1.bed simple_07_peaks.2.bed --more-bed-multiple-overlap
- |-- 11:08-WARNING : Converting to bed6 format (simple_07_peaks.bed).
- |-- 11:08-WARNING : Converting to bed6 format (simple_07_peaks.1.bed).
- |-- 11:08-WARNING : Converting to bed6 format (simple_07_peaks.2.bed).
- |-- 11:08-WARNING-ologram : Using only 1 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
- |-- 11:08-WARNING-ologram : --more-bed-labels was not set, automatically defaulting to --more-bed file names.
- |-- 11:08-WARNING-ologram : [Query + simple_07_peaks_1 + ... ]: there may be a poor fit for this feature. Check fit quality in the results. This is likely due to there being too few regions.
- |-- 11:08-WARNING-ologram : [Query + simple_07_peaks_1 + simple_07_peaks_2 + ... ]: there may be a poor fit for this feature. Check fit quality in the results. This is likely due to there being too few regions.
- |-- 11:08-WARNING-ologram : Computing log(p-val) for a Neg Binom with mean >= var ; var was set to mean+1 ([Query + simple_07_peaks_1 + simple_07_peaks_2 + ... ])
+ |-- 13:26-WARNING : Converting to bed6 format (simple_07_peaks.bed).
+ |-- 13:26-WARNING : Converting to bed6 format (simple_07_peaks.1.bed).
+ |-- 13:26-WARNING : Converting to bed6 format (simple_07_peaks.2.bed).
+ |-- 13:26-WARNING-ologram : Using only 1 threads, but 16 cores are available. Consider changing the --nb-threads parameter.
+ |-- 13:26-WARNING-ologram : --more-bed-labels was not set, automatically defaulting to --more-bed file names.
+ |-- 13:26-WARNING-ologram : [Query + simple_07_peaks_1 + ... ]: there may be a poor fit for this feature. Check fit quality in the results. This is likely due to there being too few regions.
+ |-- 13:26-WARNING-ologram : [Query + simple_07_peaks_1 + simple_07_peaks_2 + ... ]: there may be a poor fit for this feature. Check fit quality in the results. This is likely due to there being too few regions.
+ |-- 13:26-WARNING-ologram : Computing log(p-val) for a Neg Binom with mean >= var ; var was set to mean+1 ([Query + simple_07_peaks_1 + simple_07_peaks_2 + ... ])
Detailed example:
@@ -680,7 +680,7 @@$ gtftk random_list -n 3 -i simple.gtf | gtftk count
transcript 3
-exon 6
-CDS 5
+exon 3
+CDS 3
Arguments:
@@ -479,12 +479,12 @@$ gtftk random_tx -m 1 -i simple.gtf | gtftk select_by_key -k feature -v gene,transcript| gtftk tabulate -k gene_id,transcript_id
gene_id transcript_id
-G0001 G0001T001
+G0001 G0001T002
G0002 G0002T001
G0003 G0003T001
-G0004 G0004T001
+G0004 G0004T002
G0005 G0005T001
-G0006 G0006T002
+G0006 G0006T001
G0007 G0007T001
G0008 G0008T001
G0009 G0009T001
@@ -805,7 +805,7 @@ Navigation