R 레벨 레벨 축 조정

Question

"격자"패키지의 레벨 플롯 기능을 사용하여 확률 분포 함수 (PDF)를 R의 히트 맵으로 플롯하려고합니다. PDF를 함수로 구현 한 다음 값 범위와 외부 함수에 대한 두 벡터를 사용하여 레벨 플롯에 대한 행렬을 생성했습니다. 축에 행 또는 열 수 대신 두 개의 실제 값 범위를 표시하는 두 축에 적절히 간격을 둔 눈금을 추가 할 수 없다는 내 문제를 표시하려고합니다.

: 열과 행 I는 (다른 R-함수에서 사용 된)하지만 값이 로우 & 열의 숫자 눈금 자연스럽게 이격되어 다음의 플롯을 수신 이름을 지정하지 않고

# PDF to plot heatmap 
P_RCAconst <- function(x,tt,D) 
    { 
    1/sqrt(2*pi*D*tt)*1/x*exp(-(log(x) - 0.5*D*tt)^2/(2*D*tt)) 
    } 

# value ranges & computation of matrix to plot 
tt_log <- seq(-3,3,0.05) 
tt <- exp(tt_log) 
tt <- c(0,tt) 
x <- seq(0,8,0.05) 
z <- outer(x,tt,P_RCAconst, D=1.0) 
z[,1] <- 0 
z[which(x == 1),1] <- 1.5 
z[1,] <- 0.1 

# plot heatmap using levelplot 
require("lattice") 
colnames(z) <- round(tt, 2) 
rownames(z) <- x 
levelplot(z, cex.axis=1.5, cex.lab=1.5, col.regions=colorRampPalette(c("blue", "yellow","red", "black")), at=seq(0,1.9,length=200), xlab="x", ylab="time t", main="PDF P(x,t)") 

행과 열 내가 눈금이 전혀 읽을 수 있지만하지 않는 다음과 같은 플롯을받을 수있는 이름을 할당으로

적어도 실제 값에 해당합니다

나는이 겉보기에 사소한 문제에 너무 많은 시간을 할애하여 여러분의 도움을 많이 주셔서 감사합니다!

Answer 1

어쩌면이 기본 예제가 도움이 될 수 있습니다,

d = data.frame(x=rep(seq(0, 10, length=nrow(volcano)), ncol(volcano)), 
       y=rep(seq(0, 100, length=ncol(volcano)), each=nrow(volcano)), 
       z=c(volcano)) 

library(lattice) 

# let lattice pick pretty tick positions 
levelplot(z~x*y, data=d) 

# specific tick positions and labels 
levelplot(z~x*y, data=d, scales=list(x=list(at=c(2, 5, 7, 9), 
              labels=letters[1:4]))) 

# log scale 
levelplot(z~x*y, data=d, scales=list(y=list(log=TRUE))) 

틀림없이이 문서의 긴 페이지가 있지만 그것은 모두 ?xyplot 설명입니다. 마침내 levelplot 및 data.frame 객체를 사용하여 내 함수 P (X, t)의 히트 맵 플롯을 구현하는 방법

Answer 2

이것을 (뱁에서 제공하는 솔루션에 따라)

# PDF to plot heatmap 
P_RCAconst <- function(xx,tt,D) 
    { 
    1/sqrt(2*pi*D*tt)*1/xx*exp(-(log(xx) - 0.5*D*tt)^2/(2*D*tt)) 
    } 

# value ranges & computation of matrix to plot 
tt_end <- 20        # set end time 
xx_end <- 8        # set spatial boundary 
tt <- seq(0,tt_end,0.01)     # variable for time 
xx <- seq(0,xx_end,0.01)     # spatial variable 
zz <- outer(xx,tt,P_RCAconst, D=1.0)  # meshgrid for P(x,t) 
zz[,1] <- 0        # set initial condition 
zz[which(xx == 1),1] <- 4.0    # set another initial condition 
zz[1,] <- 0.1       # set boundary condition 

# plot heatmap using levelplot 
setwd("/Users/...")      # set working dirctory 
png(filename="heatmapfile.png", width=500, height=500) #or: x11() 
par(oma=c(0,1.0,0,0))     # set outer margins 
require("lattice")      # load package "lattice" 
d = data.frame(x=rep(seq(0, xx_end, length=nrow(zz)), ncol(zz)), 
       y=rep(seq(0, tt_end, length=ncol(zz)), each=nrow(zz)), 
       z=c(zz)) 
levelplot(z~x*y, data=d, cex.axis=1.5, cex.lab=1.5, cex.main=1.5, col.regions=colorRampPalette(c("blue", "yellow","red", "black")), at=c(seq(0,0.5,0.01),seq(0.55,1.4,0.05),seq(1.5,4.0,0.1)), xlab="x", ylab="time t", main="PDF P(x,t) for RC A, constant population") 
dev.off() 

그리고 어떻게 최종 내 함수 P (X, t)의 히트 맵의 줄거리는 다음과 같습니다