Как изменить расстояние или поворот заголовка оси графика относительно самой оси - PullRequest
Новая
       51

Как изменить расстояние или поворот заголовка оси графика относительно самой оси

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

Я рисую график ниже, используя пакет plot3D из R. Я хотел бы немного переместить метку оси z ($ f (y1, y2)) со значений оси или оставить ее в горизонтальном положении.Кто-нибудь знает, как я могу действовать?Далее код: 

library("plot3D")
#Function density probability
library(pbivnorm)
bsb <- function(t1,t2){
  a1 <- sqrt(phi1/2)*(sqrt(((phi1+1)*t1)/(phi1*mu1))-sqrt(((phi1*mu1)/((phi1+1)*t1))))
  a2 <- sqrt(phi2/2)*(sqrt(((phi2+1)*t2)/(phi2*mu2))-sqrt(((phi2*mu2)/((phi2+1)*t2))))
  Phi2 <- pbivnorm(a1, a2, rho, recycle = TRUE)
  b1 <- ((phi1+1)/(2*phi1*mu1))*sqrt(phi1/2)*(((phi1*mu1)/((phi1+1)*t1))^(1/2)+((phi1*mu1)/((phi1+1)*t1))^(3/2))
  b2 <- ((phi2+1)/(2*phi2*mu2))*sqrt(phi2/2)*(((phi2*mu2)/((phi2+1)*t2))^(1/2)+((phi2*mu2)/((phi2+1)*t2))^(3/2))
  fdp <- Phi2*b1*b2
  return(fdp)
}

t1 <- seq(0.001,100,length=30)
t2 <- seq(0.001,20,length=40)
#Parameters
mu1=7
phi1=2
mu2=1
phi2=9
rho=0

z<-outer(t1,t2,bsb) # calculate density values
pmat=persp3D(t1, t2, z,
        main="",xlab="$y_{1}$",ylab="$y_{2}$",zlab="$f(y_{1},y_{2})$",cex.axis=1,cex.lab=1,
        col = "gray10",border = "gray40",
        theta=50, phi=15,
        expand=0.9,
        d=2,
        shade=0.3,
        ticktype="detailed",
        nticks=5,
        facets=FALSE,contour = list(nlevels=10,col="gray35"), zlim= c(-0.1, 0.23), bty = "b2")
text(trans3d(0,7.3,0.173,pmat), "(b)",cex=1,col="black")

Картинка выглядит так:

enter image description here

Я хотел бы следующие фотографии:

enter image description here

1 Ответ

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

Проблема с размером окна вывода.Если вы используете больший, текст будет выглядеть лучше.Например, в полноэкранном режиме я получаю это из вашего исходного кода:

enter image description here

Вы также можете указать меньшее значение dpi при открытииокно, если ваше графическое устройство поддерживает это.Например, если я использую dev.new(dpi = 50), я получаю

enter image description here

Я не думаю, что есть способ повернуть метки в persp3D, но выможете нарисовать график без метки, затем используйте text, чтобы добавить метку.Вам также нужно будет увеличить размер поля на этой стороне.Например, 

par(mar = c(5.1, 9.1, 4.1, 2.1))
pmat <- persp3D(t1, t2, z, main="", xlab="$y_{1}$", ylab="$y_{2}$", 
    zlab="", cex.axis=1, cex.lab=1,
    col = "gray10", border = "gray40",
    theta=50, phi=15,
    expand=0.9, d=2, shade=0.3,
    ticktype="detailed", nticks=5,
    facets=FALSE, contour = list(nlevels=10,col="gray35"), 
    zlim= c(-0.1, 0.23), bty = "b2")
text(trans3d(0,7.3,0.173,pmat), "(b)", cex=1, col="black")
text(trans3d(0,-3,0.05,pmat), label= "$f(y_{1},y_{2})$", 
     cex=1, col="black", xpd=NA, pos=2)

Это дает мне:

enter image description here

...