R 轴图上的实际数据

Question

我会用 R 打印轴上的实际数据。如您所见，$\beta$ 在 x 轴上从 -0.5 到 0.5，但 R 绘制了点数。

library(lattice)
T <- 40
r <- 0.02
sigma <- 0.1
w <- 100
w0 <- 20
P <- 1.0
gamm <- 0.05

t <- 30

sb <- 0.10
prop <- function(beta, mu){
  (mu - r) * (w - w0 - T) * exp(-r * (T - t)) * 
    (mu - (mu - r) *  exp(-2 * beta * r * (T - t))) / 
    (gamm * r * sigma^2 * sb * (w - w0 - t))
}

beta <- seq(-0.5, 0.5, length= 10)
mu <- seq(0.04, 0.09, length= 10)

z <- outer(beta, mu, prop)

#this removes a border from the fig
trellis.par.set("axis.line", list(col="transparent"))

wireframe(z, drape=TRUE, xlab = expression(paste(beta)), ylab=expression(mu),
          zlab=expression(pi(t)), scale=list(arrows=FALSE))

Answer 1

您可以为矩阵 z 指定行名和列名 - 这将更改沿轴显示的内容：

z <- outer(beta, mu, prop)
rownames(z) <- round(beta, digits=2)
colnames(z) <- round(mu, digits=2)

#this removes a border from the fig
trellis.par.set("axis.line", list(col="transparent"))

wireframe(z, drape=TRUE, xlab = expression(paste(beta)), ylab=expression(mu),
          zlab=expression(pi(t)), scale=list(arrows=FALSE))

R 轴图上的实际数据

Actual data on R axis plot

r

lattice