如何排列两个热图的色阶（plot3D软件包中的image2D函数）？

Question

编辑：原来是densite_cauchy_norm中一个放错了括号的“歪斜”结果，现在一切正常。
我需要比较三个3D直方图（热图）。
由于某种原因，其中一个（在封闭的图像上是第二个）将值映射到与其他两个颜色不同的颜色（例如，第二个直方图上的密度为0.020的颜色与其他两个上的密度为0.015的颜色相同）。

这些来自蒙地卡罗模拟，在这里我需要将通过具有两个不同密度包络（图1和图2）的拒绝采样模拟的2个样本（X_1, X_2）的直方图与实际的“密度”（{{ 1}}（未归一化）进行采样（图3）。 因此，在它们都具有相同的颜色映射之前，我无法真正表明它们与理论的接近程度。

对于前两个直方图，我对样本X_1和X_2使用了相同的命令。

看下面的代码；您可以在底部找到数据：

f_tilde

第二个直方图应该近似于最佳，但没有其他颜色的相同色标，这似乎是完全错误的。

制作可重现示例的数据和所需功能：

library(plot3D)  
par(mfcol = c(1, 3)); 

n <- 100000;  
X_1 <- simu_f_1(n);  
X_2<- simu_f_2(n);  
x_s <- 50;  y_s <- 50;    

mon_histo(X_1, x_s, y_s, opt_3D = FALSE);  
mon_histo(X_2, x_s, y_s, opt_3D = FALSE);  

z <- outer(seq.int(0, 4, length= x_s), seq.int(0, 2, length = y_s), f_tilde);

image2D(z=z/sum(z), x = seq.int(0, 4, length= x_s),y = seq.int(0, 2, 
                    length = y_s), border="black", contour = FALSE);