我正在使用函数绘制曲面图persp3D()
。但是,我希望 y 轴 (k) 刻度上的零与我的 z 轴 (R0) 从相同的零开始。目前,我的 y 轴(k)和 x 轴(A)上的零从同一点开始。我需要 y 轴和 z 轴上的零从同一点开始,以便获得更好的观察方向。我附上了两张图片,一张是我的结果,另一张是我想要的情节的一个例子。下面也是我的代码。
library(plot3D)
lambda <- 1/12
sigma <- 0.4
gamma <- 0.4
beta <- 0.00006
c <-0.0034
m <- 0.056/12
mu <- 1
f <- function(x,y){(beta*(x*(m+lambda)/((m+c)*(m+lambda)+m*y))*sigma*gamma)/((m+sigma+c)*(m+mu+c))}
x <- seq(0,1000,10)
y <- seq(0,1,0.01)
z <- outer(x,y,f)
persp3D(x=x, y=y, z=z,
col.palette = heat.colors,
phi = 15, theta = 295, nticks=5,
zlab="", ylab = "k", xlab="A",
xlim = c(0, 1000),
ylim = c(0, 1),
zlim = c(0, 3.5),
box = TRUE, border = NA, shade = .4,ticktype="detailed",scale= TRUE, expand = 0.5)
par(xpd = TRUE)
text3D(0, 1.3, 1, expression(R[0]), add = T)