lets-plot 4.5.3a1__cp312-cp312-win_amd64.whl → 4.6.0a2__cp312-cp312-win_amd64.whl

lets_plot/plot/__init__.py

@@ -13,6 +13,7 @@ from .geom_extras import *
 from .geom_function_ import *
 from .geom_imshow_ import *
 from .geom_livemap_ import *
+from .ggbunch_ import *
 from .gggrid_ import *
 from .ggtb_ import *
 from .guide import *
@@ -56,6 +57,7 @@ __all__ = (coord.__all__ +
            annotation.__all__ +
            marginal_layer.__all__ +
            font_features.__all__ +
+           ggbunch_.__all__ +
            gggrid_.__all__ +
            ggtb_.__all__ +
            expand_limits_.__all__

lets_plot/plot/core.py

@@ -706,7 +706,7 @@ class LayerSpec(FeatureSpec):
                 map_data_meta = {'georeference': {}}
             else:
                 # Fetch proper GeoDataFrame. Further processing is the same as if map was a GDF.
-                if name in ['point', 'pie', 'text', 'livemap']:
+                if name in ['point', 'pie', 'text', 'label', 'livemap']:
                     map = map.get_centroids()
                 elif name in ['map', 'polygon']:
                     map = map.get_boundaries()

lets_plot/plot/geom.py

@@ -12,7 +12,7 @@ from .util import as_annotated_data, is_geo_data_frame, geo_data_frame_to_crs, g
 #
 __all__ = ['geom_point', 'geom_path', 'geom_line',
            'geom_smooth', 'geom_bar',
-           'geom_histogram', 'geom_dotplot', 'geom_bin2d',
+           'geom_histogram', 'geom_dotplot', 'geom_bin2d', 'geom_hex',
            'geom_tile', 'geom_raster',
            'geom_errorbar', 'geom_crossbar', 'geom_linerange', 'geom_pointrange',
            'geom_contour',
@@ -128,7 +128,7 @@ def geom_point(mapping=None, *, data=None, stat=None, position=None, show_legend
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Point` and `MultiPoint`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `centroids()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_centroids()` function.
 
     ----
 
@@ -193,7 +193,7 @@ def geom_point(mapping=None, *, data=None, stat=None, position=None, show_legend
 
     .. jupyter-execute::
         :linenos:
-        :emphasize-lines: 8-11
+        :emphasize-lines: 8-12
 
         from lets_plot import *
         from lets_plot.geo_data import *
@@ -205,7 +205,8 @@ def geom_point(mapping=None, *, data=None, stat=None, position=None, show_legend
             geom_point(aes(size="est_pop_2019"), color="red", show_legend=False, \\
                        data=data, map=centroids, map_join="city", \\
                        tooltips=layer_tooltips().title("@city")
-                                                .line("population|@est_pop_2019"))
+                                                .line("population|@est_pop_2019")
+                                                .format("@est_pop_2019", ".2~s"))
 
     """
     return _geom('point',
@@ -1420,6 +1421,173 @@ def geom_bin2d(mapping=None, *, data=None, stat=None, position=None, show_legend
                  **other_args)
 
 
+def geom_hex(mapping=None, *, data=None, stat=None, position=None, show_legend=None, inherit_aes=None,
+             manual_key=None, sampling=None,
+             tooltips=None,
+             bins=None,
+             binwidth=None,
+             drop=None,
+             color_by=None, fill_by=None,
+             **other_args):
+    """
+    Apply a hexagonal grid to the plane, count observations in each cell (hexagonal bin) of the grid,
+    and map the count to the fill color of the cell (hexagonal tile).
+
+    By default, this geom uses `coord_fixed()`.
+    However, this may not be the best choice when the values on the X/Y axis have significantly different magnitudes.
+    In such cases, try using `coord_cartesian()`.
+
+    Parameters
+    ----------
+    mapping : `FeatureSpec`
+        Set of aesthetic mappings created by `aes()` function.
+        Aesthetic mappings describe the way that variables in the data are
+        mapped to plot "aesthetics".
+    data : dict or Pandas or Polars `DataFrame`
+        The data to be displayed in this layer. If None, the default, the data
+        is inherited from the plot data as specified in the call to ggplot.
+    stat : str, default='binhex'
+        The statistical transformation to use on the data for this layer, as a string.
+    position : str or `FeatureSpec`, default='identity'
+        Position adjustment.
+        Either a position adjustment name: 'dodge', 'dodgev', 'jitter', 'nudge', 'jitterdodge', 'fill',
+        'stack' or 'identity', or the result of calling a position adjustment function (e.g., `position_dodge()` etc.).
+    show_legend : bool, default=True
+        False - do not show legend for this layer.
+    inherit_aes : bool, default=True
+        False - do not combine the layer aesthetic mappings with the plot shared mappings.
+    manual_key : str or `layer_key`
+        The key to show in the manual legend.
+        Specify text for the legend label or advanced settings using the `layer_key()` function.
+    sampling : `FeatureSpec`
+        Result of the call to the `sampling_xxx()` function.
+        To prevent any sampling for this layer pass value "none" (string "none").
+    tooltips : `layer_tooltips`
+        Result of the call to the `layer_tooltips()` function.
+        Specify appearance, style and content.
+        Set tooltips='none' to hide tooltips from the layer.
+    bins : list of int, default=[30, 30]
+        Number of hexagonal bins in both directions, vertical and horizontal. Overridden by `binwidth`.
+    binwidth : list of float
+        The width of the hexagonal bins in both directions, vertical and horizontal.
+        Override `bins`. The default is to use bin widths that cover the entire range of the data.
+    drop : bool, default=True
+        Specify whether to remove all hexagonal bins with 0 counts.
+    color_by : {'fill', 'color', 'paint_a', 'paint_b', 'paint_c'}, default='color'
+        Define the color aesthetic for the geometry.
+    fill_by : {'fill', 'color', 'paint_a', 'paint_b', 'paint_c'}, default='fill'
+        Define the fill aesthetic for the geometry.
+    other_args
+        Other arguments passed on to the layer.
+        These are often aesthetics settings used to set an aesthetic to a fixed value,
+        like color='red', fill='blue', size=3 or shape=21.
+        They may also be parameters to the paired geom/stat.
+
+    Returns
+    -------
+    `LayerSpec`
+        Geom object specification.
+
+    Notes
+    -----
+
+    Computed variables:
+
+    - ..count.. : number of points with coordinates in the same hexagonal bin.
+
+    `geom_hex()` understands the following aesthetics mappings:
+
+    - x : x-axis value.
+    - y : y-axis value.
+    - alpha : transparency level of a layer. Accept values between 0 and 1.
+    - color (colour) : color of the geometry lines. For more info see `Color and Fill <https://lets-plot.org/python/pages/aesthetics.html#color-and-fill>`__.
+    - fill : fill color. For more info see `Color and Fill <https://lets-plot.org/python/pages/aesthetics.html#color-and-fill>`__.
+    - size : line width, default=0 (i.e. tiles outline initially is not visible).
+    - weight : used by 'binhex' stat to compute weighted sum instead of simple count.
+    - width : width of the hexagon.
+    - height : the real height of the hexagon will be 2/sqrt(3) times this value, so with width=height the hexagon will be the regular.
+
+    ----
+
+    To hide axis tooltips, set 'blank' or the result of `element_blank()`
+    to the `axis_tooltip`, `axis_tooltip_x` or `axis_tooltip_y` parameter of the `theme()`.
+
+    Examples
+    --------
+    .. jupyter-execute::
+        :linenos:
+        :emphasize-lines: 8
+
+        import numpy as np
+        from lets_plot import *
+        LetsPlot.setup_html()
+        np.random.seed(42)
+        mean = np.zeros(2)
+        cov = np.eye(2)
+        x, y = np.random.multivariate_normal(mean, cov, 1000).T
+        ggplot({'x': x, 'y': y}, aes(x='x', y='y')) + geom_hex()
+
+    |
+
+    .. jupyter-execute::
+        :linenos:
+        :emphasize-lines: 9-14
+
+        import numpy as np
+        from lets_plot import *
+        LetsPlot.setup_html()
+        np.random.seed(42)
+        n = 5000
+        x = np.random.uniform(-2, 2, size=n)
+        y = np.random.normal(scale=.5, size=n)
+        ggplot({'x': x, 'y': y}, aes(x='x', y='y')) + \\
+            geom_hex(aes(fill='..density..'), binwidth=[.25, .24], \\
+                     tooltips=layer_tooltips().format('@x', '.2f')
+                             .format('@y', '.2f').line('(@x, @y)')
+                             .line('count|@..count..')
+                             .format('@..density..', '.3f')
+                             .line('density|@..density..')) + \\
+            scale_fill_gradient(low='black', high='red')
+
+    |
+
+    .. jupyter-execute::
+        :linenos:
+        :emphasize-lines: 10-11
+
+        import numpy as np
+        from lets_plot import *
+        LetsPlot.setup_html()
+        np.random.seed(42)
+        mean = np.zeros(2)
+        cov = [[1, .5],
+               [.5, 1]]
+        x, y = np.random.multivariate_normal(mean, cov, 500).T
+        ggplot({'x': x, 'y': y}, aes(x='x', y='y')) + \\
+            geom_hex(aes(alpha='..count..'), bins=[20, 20], \\
+                     fill='darkgreen') + \\
+            geom_point(size=1.5, shape=21, color='white', \\
+                       fill='darkgreen') + \\
+            ggsize(600, 450)
+
+    """
+    return _geom('hex',
+                 mapping=mapping,
+                 data=data,
+                 stat=stat,
+                 position=position,
+                 show_legend=show_legend,
+                 inherit_aes=inherit_aes,
+                 manual_key=manual_key,
+                 sampling=sampling,
+                 tooltips=tooltips,
+                 bins=bins,
+                 binwidth=binwidth,
+                 drop=drop,
+                 color_by=color_by, fill_by=fill_by,
+                 **other_args)
+
+
 def geom_tile(mapping=None, *, data=None, stat=None, position=None, show_legend=None, inherit_aes=None,
               manual_key=None, sampling=None,
               tooltips=None,
@@ -2617,7 +2785,7 @@ def geom_polygon(mapping=None, *, data=None, stat=None, position=None, show_lege
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Polygon` and `MultiPolygon`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `boundaries()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_boundaries()` function.
 
     ----
 
@@ -2807,7 +2975,7 @@ def geom_map(mapping=None, *, data=None, stat=None, position=None, show_legend=N
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Polygon` and `MultiPolygon`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `boundaries()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_boundaries()` function.
 
     ----
 
@@ -4652,7 +4820,7 @@ def geom_density(mapping=None, *, data=None, stat=None, position=None, show_lege
 
     .. jupyter-execute::
         :linenos:
-        :emphasize-lines: 10
+        :emphasize-lines: 11
 
         import numpy as np
         from lets_plot import *
@@ -4660,13 +4828,17 @@ def geom_density(mapping=None, *, data=None, stat=None, position=None, show_lege
         np.random.seed(42)
         x = np.random.normal(size=1000)
         p = ggplot({'x': x}, aes(x='x'))
-        bunch = GGBunch()
-        for i, bw in enumerate([.1, .2, .4]):
-            for j, n in enumerate([16, 64, 256]):
-                bunch.add_plot(p + geom_density(kernel='epanechikov', bw=bw, n=n) + \\
-                                   ggtitle('bw={0}, n={1}'.format(bw, n)),
-                               j * 300, i * 200, 300, 200)
-        bunch.show()
+        bandwidths = [0.1, 0.2, 0.4]
+        sample_sizes = [16, 64, 256]
+        plots = [
+            p + geom_density(
+                kernel='epanechikov', bw=bw, n=n
+            ) + ggtitle(f'bw={bw}, n={n}')
+            for bw in bandwidths
+            for n in sample_sizes
+        ]
+
+        gggrid(plots, ncol=3) + ggsize(900, 600)
 
     |
 
@@ -4681,13 +4853,15 @@ def geom_density(mapping=None, *, data=None, stat=None, position=None, show_lege
         x = np.random.normal(size=1000)
         y = np.sign(x)
         p = ggplot({'x': x, 'y': y}, aes(x='x'))
-        bunch = GGBunch()
-        for i, adjust in [(i, .5 * (1 + i)) for i in range(3)]:
-            bunch.add_plot(p + geom_density(aes(weight='y'), kernel='cosine', \\
-                                            adjust=adjust) + \\
-                               ggtitle('adjust={0}'.format(adjust)),
-                           i * 300, 0, 300, 200)
-        bunch.show()
+        adjustments = [0.5 * (1 + i) for i in range(3)]
+        plots = [
+            p + geom_density(aes(weight='y'),
+                kernel='cosine', adjust=adjust
+            ) + ggtitle(f'adjust={adjust}')
+            for adjust in adjustments
+        ]
+
+        gggrid(plots) + ggsize(800, 200)
 
     """
     return _geom('density',
@@ -4847,7 +5021,7 @@ def geom_density2d(mapping=None, *, data=None, stat=None, position=None, show_le
 
     .. jupyter-execute::
         :linenos:
-        :emphasize-lines: 12
+        :emphasize-lines: 13
 
         import numpy as np
         from lets_plot import *
@@ -4857,19 +5031,23 @@ def geom_density2d(mapping=None, *, data=None, stat=None, position=None, show_le
         x = np.random.normal(size=n)
         y = np.random.normal(size=n)
         p = ggplot({'x': x, 'y': y}, aes('x', 'y'))
-        bunch = GGBunch()
-        for i, bw in enumerate([.2, .4]):
-            for j, n in enumerate([16, 256]):
-                bunch.add_plot(p + geom_density2d(kernel='epanechikov', bw=bw, n=n) + \\
-                                   ggtitle('bw={0}, n={1}'.format(bw, n)),
-                               j * 400, i * 400, 400, 400)
-        bunch.show()
+        bandwidths = [0.2, 0.4]
+        sample_sizes = [16, 256]
+        plots = [
+            p + geom_density2d(kernel='epanechikov',
+                               bw=bw, n=n
+            ) + ggtitle(f'bw={bw}, n={n}')
+            for bw in bandwidths
+            for n in sample_sizes
+        ]
+
+        gggrid(plots, ncol=2) + ggsize(600, 650)
 
     |
 
     .. jupyter-execute::
         :linenos:
-        :emphasize-lines: 12-13
+        :emphasize-lines: 13-14
 
         import numpy as np
         from lets_plot import *
@@ -4879,14 +5057,18 @@ def geom_density2d(mapping=None, *, data=None, stat=None, position=None, show_le
         x = np.random.normal(size=n)
         y = np.random.normal(size=n)
         p = ggplot({'x': x, 'y': y}, aes('x', 'y'))
-        bunch = GGBunch()
-        for i, adjust in enumerate([1.5, 2.5]):
-            for j, bins in enumerate([5, 15]):
-                bunch.add_plot(p + geom_density2d(kernel='cosine', \\
-                                                  adjust=adjust, bins=bins) + \\
-                                   ggtitle('adjust={0}, bins={1}'.format(adjust, bins)),
-                               j * 400, i * 400, 400, 400)
-        bunch.show()
+        adjustments = [1.5, 2.5]
+        bin_counts = [5, 15]
+        plots = [
+            p + geom_density2d(kernel='cosine',
+                adjust=adjust,
+                bins=bins
+            ) + ggtitle(f'adjust={adjust}, bins={bins}')
+            for adjust in adjustments
+            for bins in bin_counts
+        ]
+
+        gggrid(plots, ncol=2) + ggsize(600, 650)
 
     |
 
@@ -5073,7 +5255,7 @@ def geom_density2df(mapping=None, *, data=None, stat=None, position=None, show_l
 
     .. jupyter-execute::
         :linenos:
-        :emphasize-lines: 12-13
+        :emphasize-lines: 13
 
         import numpy as np
         from lets_plot import *
@@ -5083,14 +5265,16 @@ def geom_density2df(mapping=None, *, data=None, stat=None, position=None, show_l
         x = np.random.normal(size=n)
         y = np.random.normal(size=n)
         p = ggplot({'x': x, 'y': y}, aes(x='x', y='y'))
-        bunch = GGBunch()
-        for i, bw in enumerate([.2, .4]):
-            for j, n in enumerate([16, 256]):
-                bunch.add_plot(p + geom_density2df(kernel='epanechikov', bw=bw, n=n, \\
-                                                   size=.5, color='white') + \\
-                                   ggtitle('bw={0}, n={1}'.format(bw, n)),
-                               j * 400, i * 400, 400, 400)
-        bunch.show()
+        bandwidths = [0.2, 0.4]
+        sample_sizes = [16, 256]
+        plots = [
+            p + geom_density2df(kernel='epanechikov', size=.5, color='white',
+                               bw=bw, n=n
+            ) + ggtitle(f'bw={bw}, n={n}')
+            for bw in bandwidths
+            for n in sample_sizes
+        ]
+        gggrid(plots, ncol=2) + ggsize(600, 650)
 
     |
 
@@ -5106,15 +5290,18 @@ def geom_density2df(mapping=None, *, data=None, stat=None, position=None, show_l
         x = np.random.normal(size=n)
         y = np.random.normal(size=n)
         p = ggplot({'x': x, 'y': y}, aes(x='x', y='y'))
-        bunch = GGBunch()
-        for i, adjust in enumerate([1.5, 2.5]):
-            for j, bins in enumerate([5, 15]):
-                bunch.add_plot(p + geom_density2df(kernel='cosine', \\
-                                                   size=.5, color='white', \\
-                                                   adjust=adjust, bins=bins) + \\
-                                   ggtitle('adjust={0}, bins={1}'.format(adjust, bins)),
-                               j * 400, i * 400, 400, 400)
-        bunch.show()
+        adjustments = [1.5, 2.5]
+        bin_counts = [5, 15]
+        plots = [
+            p + geom_density2df(kernel='cosine', size=.5, color='white',
+                adjust=adjust,
+                bins=bins
+            ) + ggtitle(f'adjust={adjust}, bins={bins}')
+            for adjust in adjustments
+            for bins in bin_counts
+        ]
+
+        gggrid(plots, ncol=2) + ggsize(600, 650)
 
     |
 
@@ -6154,7 +6341,7 @@ def geom_rect(mapping=None, *, data=None, stat=None, position=None, show_legend=
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `MultiPoint`, `Line`, `MultiLine`, `Polygon` and `MultiPolygon`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `limits()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_limits()` function.
 
     ----
 
@@ -6820,7 +7007,7 @@ def geom_text(mapping=None, *, data=None, stat=None, position=None, show_legend=
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Point` and `MultiPoint`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `centroids()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_centroids()` function.
 
     ----
 
@@ -7067,7 +7254,7 @@ def geom_label(mapping=None, *, data=None, stat=None, position=None, show_legend
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Point` and `MultiPoint`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `centroids()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_centroids()` function.
 
     ----
 
@@ -7305,7 +7492,7 @@ def geom_pie(mapping=None, *, data=None, stat=None, position=None, show_legend=N
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Point` and `MultiPoint`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `centroids()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_centroids()` function.
 
     ----
 
@@ -7814,7 +8001,7 @@ def geom_blank(mapping=None, *, data=None, stat=None, position=None, show_legend
 
     The `data` and `map` parameters of `GeoDataFrame` type support shapes `Point` and `MultiPoint`.
 
-    The `map` parameter of `Geocoder` type implicitly invokes `centroids()` function.
+    The `map` parameter of `Geocoder` type implicitly invokes `get_centroids()` function.
 
     ----
 

lets_plot/plot/ggbunch_.py

@@ -0,0 +1,96 @@
+#  Copyright (c) 2025. JetBrains s.r.o.
+#  Use of this source code is governed by the MIT license that can be found in the LICENSE file.
+
+from numbers import Number
+from typing import List, Tuple
+
+from ._global_theme import _get_global_theme
+from .subplots import SupPlotsLayoutSpec
+from .subplots import SupPlotsSpec
+from .subplots_util import _strip_theme_if_global
+
+__all__ = ['ggbunch']
+
+
+def ggbunch(plots: List,
+            regions: List[Tuple[float, float, float, float, float, float]]
+            ) -> SupPlotsSpec:
+    """
+    Combine several plots into a single figure with custom layout.
+
+    Parameters
+    ----------
+    plots : List
+        A list where each element is one of:
+
+        - a plot specification
+        - a subplots specification
+        - None
+
+    regions : List[Tuple]
+        Layout parameters for each plot. Each region is specified as
+        (x, y, width, height, dx, dy) where:
+
+        - x, y: Position of the plot's top-left corner in relative coordinates ([0,0] is top-left corner, [1,1] is bottom-right corner of the container).
+        - width, height: Size of the plot relative to container dimensions (1 equal to the full container width/height).
+        - dx, dy: Pixel offsets to move the region (defaults to 0).
+
+    Returns
+    -------
+    `SupPlotsSpec`
+        A specification describing the combined figure with all plots and their layout.
+
+    Examples
+    --------
+    .. jupyter-execute::
+        :linenos:
+        :emphasize-lines: 10-14
+
+        import numpy as np
+        from lets_plot import *
+        LetsPlot.setup_html()
+        np.random.seed(42)
+        data = {'x': np.random.gamma(2.0, size=100)}
+        p1 = ggplot(data, aes(x='x')) + \\
+            geom_histogram(aes(color='x', fill='x'))
+        p2 = ggplot(data, aes(x='x')) + \\
+            geom_density() + theme_bw() + theme(axis='blank', panel_grid='blank')
+        ggbunch(
+            [p1, p2],
+            [(0, 0, 1, 1),
+             (0.5, 0.1, 0.3, 0.3)]
+        ) + ggsize(400, 300)
+
+    """
+
+    if not len(plots):
+        raise ValueError("Supplots list is empty.")
+
+    # Validate provided regions
+    for i, region in enumerate(regions):
+        if len(region) not in (4, 6):
+            raise ValueError(f"Region {i} must have 4 or 6 values, got {len(region)}")
+        if not all(isinstance(x, Number) for x in region):
+            raise ValueError(f"Region {i} contains non-numeric values: {region}")
+
+        # Validate size is positive
+        if any(x <= 0 for x in region[2:4]):
+            raise ValueError(f"Region {i} sizes must be positive: {region}")
+
+    # Convert regions tuples to lists
+    regions_list = [list(r) for r in regions]
+    layout = SupPlotsLayoutSpec(
+        name="free",
+        regions=regions_list
+    )
+
+    figures = [_strip_theme_if_global(fig) for fig in plots]
+
+    figure_spec = SupPlotsSpec(figures=figures, layout=layout)
+
+    # Apply global theme if defined
+    global_theme_options = _get_global_theme()
+    if global_theme_options is not None:
+        figure_spec += global_theme_options
+
+    return figure_spec

lets_plot/plot/gggrid_.py

@@ -3,10 +3,10 @@
 # Use of this source code is governed by the MIT license that can be found in the LICENSE file.
 #
 
-from lets_plot.plot.core import PlotSpec
 from ._global_theme import _get_global_theme
 from .subplots import SupPlotsLayoutSpec
 from .subplots import SupPlotsSpec
+from .subplots_util import _strip_theme_if_global
 
 __all__ = ['gggrid']
 
@@ -114,27 +114,12 @@ def gggrid(plots: list, ncol: int = None, *,
         align=align
     )
 
-    # Global Theme
-    global_theme_options = _get_global_theme()
-
-    def _strip_theme_if_global(fig):
-        # Strip global theme options from plots in grid (see issue: #966).
-        if global_theme_options is not None and fig is not None and 'theme' in fig.props() and fig.props()[
-            'theme'] == global_theme_options.props():
-            if isinstance(fig, PlotSpec):
-                fig = PlotSpec.duplicate(fig)
-                fig.props().pop('theme')
-                return fig
-            elif isinstance(fig, SupPlotsSpec):
-                fig = SupPlotsSpec.duplicate(fig)
-                fig.props().pop('theme')
-                return fig
-        return fig
-
     figures = [_strip_theme_if_global(fig) for fig in plots]
 
     figure_spec = SupPlotsSpec(figures=figures, layout=layout)
 
+    # Apply global theme if defined
+    global_theme_options = _get_global_theme()
     if global_theme_options is not None:
         figure_spec += global_theme_options
 

lets_plot/plot/plot.py

@@ -142,6 +142,9 @@ def ggsize(width, height):
 
 class GGBunch(FeatureSpec):
     """
+    Class `GGBunch` is deprecated and will be removed in future releases.
+    Please, use function `ggbunch()` to combine several plots into a single figure with custom layout.
+
     Collection of plots created by ggplot function.
     Use method `add_plot()` to add plot to 'bunch'.
     Each plot can have arbitrary location and size.
@@ -176,6 +179,10 @@ class GGBunch(FeatureSpec):
         """
         super().__init__('ggbunch', None)
         self.items = []
+        print("\n(!) WARN: class GGBunch is deprecated and will be removed in future releases.\n\n"
+              "  Please, use function ggbunch() to combine several plots into\n"
+              "  a single figure with custom layout.\n")
+
 
     def add_plot(self, plot_spec: PlotSpec, x, y, width=None, height=None):
         """