Актуальные данные на графике оси R - PullRequest
0 голосов
/ 24 мая 2018

Я бы R, чтобы напечатать фактические данные на осях.Как вы можете видеть, $ \ beta $ составляет от -0,5 до 0,5 по оси x, но 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))

1 Ответ

0 голосов
/ 24 мая 2018

Вы можете указать имена строк и столбцов для вашей матрицы 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))

enter image description here

...