R 中带有 plot3D 或 plotly 的 3D 曲面图

Question

我想在 R 中生成一个 3D 表面，其中 y 是通量，x 是年龄，z 是降水量。我的年龄范围在 0 到 152 之间，降水量范围在 0 到 2600 之间。

我还有两个函数，其中通量是年龄或降水的函数：

Flux= 0.387217066*(Age^0.328086337)*(exp(-0.004177033*Age)

和

Flux= 1.117997*(1-exp(-exp(-5.426564)*(Preci-(-220.745499))

我想要实现的有点像this:

我尝试用 R 中的包 plot3D 来做，但没有成功（见下文）

Age<- as.matrix(seq(0:152)) 
Preci<-as.matrix(seq(from=10, to=2600, by=17))
Flux= as.matrix(0.387217066*(Age^0.328086337)*(exp(-0.004177033*Age)) - 1.117997*(1-exp(-exp(-5.426564)*(Preci-(-220.745499)))))
surf3D(Age, Preci, Flux, colvar = Flux, colkey = FALSE, facets = FALSE)

我收到这条错误消息

Error in if (is.na(var)) ispresent <- FALSE else if (length(var) == 1) if (is.logical(var)) if (!var) ispresent <- FALSE : 
  argument is of length zero

Answer 1

这是一个开始，使用 emdbook::curve3d() 作为 lattice::wireframe 的包装器（参见 ?curve3d 的 sys3d 参数）。对我来说，为什么将你的通量函数作为年龄函数与降水函数相减是有意义的，但正如你在上面所做的那样......

## for readability, put numeric values in a separate vector
params <- c(a0=0.387217066,a1=0.328086337,a2=0.004177033,
            p0=1.117997,p1=5.426564,p2=-220.745499)

library("emdbook")

curve3d(with(as.list(params),
          a0*(Age^a1)*exp(-a2*Age)-
            p0*(1-exp(-exp(-p1)*(Preci-p2)))),
        varnames=c("Age","Preci"),
        xlim=c(0,200),ylim=c(10,2600),
        sys3d="wireframe",
        col="gray",
        zlab="Flux")

curve3d 也是 returns 一个包含组件 $x、$y、$z 的列表，您可以将其用作其他 3D 图形框架的输入。

R 中带有 plot3D 或 plotly 的 3D 曲面图

3D surface plot in R with plot3D or plotly

r

surface

plot3d