code: def test_points():
      p = (
          p0
          + geom_point(
              aes(
                  fill=after_stat("density"), size=after_stat("density")
              ),
              stat="pointdensity",
          )
          + scale_size_radius(range=(10, 20))
      )
  
      assert p == "points"
  
imports: import numpy
  import pandas
  from plotnine import aes, after_stat, geom_point, geom_pointdensity, ggplot, scale_size_radius
label_desc: Use image node_id 2 for calling a global function with example usage: test_points() without return types