Как показать значения матрицы на Levelplot
У меня есть матричные данные здесь, и я визуализировал их с levelplot
, Участок находится ниже. Но я просто не мог поместить значения в сюжет, я имею в виду, что прочитал этот вопрос, но все еще не мог понять это.
Как я могу это сделать? Благодарю.
3 ответа
Проблема с кодом в ответе, с которым вы связаны, состоит в том, что он работает только тогда, когда объекты в формуле уровня графика названы x
, y
, а также z
,
Вот пример, который использует более стандартную идиому для обработки аргументов, передаваемых пользовательской функции панели, и поэтому становится более общеприменимым:
library("lattice")
## Example data
x <- seq(pi/4, 5*pi, length.out=10)
y <- seq(pi/4, 5*pi, length.out=10)
grid <- expand.grid(X=x, Y=y)
grid$Z <- runif(100, -1, 1)
## Write a panel function (after examining 'args(panel.levelplot) to see what
## will be being passed on to the panel function by levelplot())
myPanel <- function(x, y, z, ...) {
panel.levelplot(x,y,z,...)
panel.text(x, y, round(z,1))
}
## Try it out
levelplot(Z ~ X*Y, grid, panel = myPanel)
mat <- read.csv("J_H2S1T6_PassTraffic.csv", header=F)
y <- as.numeric(mat[1,-1])
mat <- mat[-1,-1]
n <- dim(mat)[1]
Вот модификация, я генерирую новый масштаб
x <- seq(min(y), max(y), length.out=n)
grid <- expand.grid(x=x, y=x)
mat <- as.matrix(mat)
dim(mat) <- c(n*n,1)
grid$z <- mat
Здесь модификация. Я изменяю размер матрицы на вектор, чтобы поместить его в сетку.
mat <- as.matrix(mat)
dim(mat) <- c(n*n,1)
grid$z <- mat
p <- levelplot(z~x*y, grid,
panel=function(...) {
arg <- list(...)
panel.levelplot(...)
panel.text(arg$x, arg$y,arg$z)},
scales = list(y = list(at=y,labels=y),
x = list(at=y,labels=y)))
print(p)
Другой вариант заключается в использовании layer()
от latticeExtra. Это позволяет вам наложить один сюжет поверх другого, используя +
Оператор знакомый энтузиастам ggplot2:
library(latticeExtra)
## Applied to the example data in my other answer, this will produce
## an identical plot
levelplot(Z ~ X*Y, data = grid) +
layer(panel.text(X, Y, round(Z, 1)), data = grid)