r - stat_summary2dプロットの各セル内の要素数をggplotに入れる方法は?

Question

z 変数の平均 (z_mean) が変数 x および y によってどのように変化するかを示すために、次の 2D プロットを生成しました。

上の図を生成するコードは次のとおりです。

test = data.frame(group = factor(c(rep(1, 100), rep(2, 100), rep(3, 100))),
              x = sample(20, 300, replace =T),
              y = sample(50, 300, replace =T),
              z_mean = sample(5, 300, replace =T),
              z_sd = runif(300))

ggplot(data=test, aes(x=x, y=y, z=z_mean) ) + 
stat_summary2d(fun = function(x) {mean(x)}, binwidth = c(2, 5))+
# stat_bin2d(binwidth = c(2, 5))+
scale_fill_gradientn(colours=c("green", "yellow", "red"), name = "")+
facet_wrap(~group) +
ggtitle("z_mean")

ここで、各セル内のデータ ポイントの数をテキストに入れたいと思います。それを行う方法はありますか、stat_bin2d から値を抽出している可能性があります。stat_summary2d または stat_bin2d によって抽出された値を抽出する方法はありますか?

また、変数 z の標準偏差 (z_sd) を示す別の図があります。下記参照：

これら 2 つの図を 1 つに組み合わせる創造的な方法はありますか。たぶん、長方形を分割するか、長方形内の円のサイズを調整します。ggplot でそれを行う方法が見つかりませんでした。