使用r绘制特定区域

Question

我是R代码的绿色手。现在我在使用R代码绘制某个特定区域时遇到了一些麻烦。

我的问题如下：我有一个等式取决于两个参数（a，b），我想绘制与（a，b）相关的区域，在那里我可以找到我方程的解。

等式定义为：

其中F（x）是mean = 0和sd = 1的正态分布的累积函数，f（x）是法线和f的密度函数。 id是f。

的第一个导数

我的目标是构建由：

定义的区域（Surface）

我尝试使用轮廓来解决我的问题，但我无法获得S的情节。

我成功地将等式（x）的解作为x = g（a）的函数绘制了b的所有值，这是我的代码：

E<-function(a,u) {  
((a-u)*(dnorm(q(u),m,s, log = FALSE))/(qnorm(u)))-((qnorm(u))^2)/2 +a*qnorm(u)-(a^2)/2 
}
a <- u <- seq(0,1,0.01)
z <- outer(a,u,E)
b<-seq(0,1,0.001)
contour( x=a, y=a, z=z,levels=b, las=1, drawlabels=FALSE, lwd=3,xlab="a", ylab="x")

结果如下图：

这不是我想要的，我希望有一个b的函数，a = h（b）。

有没有办法解决我的问题？

非常感谢任何帮助。干杯

谢谢！