[英]Error Switching geom_tile plot to Stat_Density_2D Plot
我一直在嘗試用以下數據制作熱圖。
數據:
structure(list(platelocheight = c(2.594, 3.803, 3.254, 3.599,
3.617, 3.297, 2.093, 3.611, 2.842, 3.316, 2.872, 3.228, 3.633,
4.28, 3.309, 2.8, 2.632, 3.754, 2.207, 3.604, 3.443, 2.188, 3.452,
2.553, 3.382, 3.067, 2.986, 2.785, 2.567, 3.804), platelocside = c(0.059,
-1.596, -0.65, -0.782, -0.301, -0.104, 0.057, -0.807, 0.003,
1.661, 0.088, -0.32, -1.115, -0.146, -0.364, -0.952, 0.254, 0.109,
-0.671, -0.803, -0.212, -0.069, -0.09, -0.472, 0.434, 0.337,
0.723, 0.508, -0.197, -0.635), exitspeed = c(69.891, 73.352,
83.942, 85.67, 79.454, 85.277, 81.078, 73.573, 77.272, 59.263,
97.343, 91.436, 76.264, 83.479, 47.576, 84.13, 60.475, 61.093,
84.54, 69.959, 88.729, 88.019, 82.18, 83.684, 86.296, 90.605,
79.945, 59.899, 62.522, 77.75)), .Names = c("platelocheight",
"platelocside", "exitspeed"), row.names = c(NA, 30L), class = "data.frame")
我首先使用此代碼成功創建了geom_tile熱圖。
library(RODBC)
library(ggplot2)
con=odbcConnect('username',uid='userid', pwd = 'password')
df=sqlQuery(con,"select platelocheight, platelocside, exitspeed from pitches_sample where pitchcall='InPlay'
and exitspeed is not null")
topKzone <- 3.5
botKzone <- 1.6
inKzone <- -0.95
outKzone <- 0.95
kZone <- data.frame(
x=c(inKzone, inKzone, outKzone, outKzone, inKzone),
y=c(botKzone, topKzone, topKzone, botKzone, botKzone)
)
df$h <- round(df$platelocheight)
df$s <- round(df$platelocside)
ggplot(kZone, aes(x,y)) +
geom_tile(data=df, aes(x=s, y=h, fill=exitspeed)) +
scale_fill_distiller(palette = "Spectral") +
geom_path(lwd=1.5, col="black") +
coord_fixed()
我喜歡情節的結構,但我認為樣式/外觀太像瓷磚了。 在理想的世界中,我想創建一個與此相似的情節。
為此,我接受了人們對stackoverflow的建議,並嘗試使用以下代碼創建stat_density_2D圖。
library(RODBC)
library(ggplot2)
con=odbcConnect('username',uid='userid', pwd = 'password')
df=sqlQuery(con,"select platelocheight, platelocside, exitspeed from pitches_sample where pitchcall='InPlay'
and exitspeed is not null")
topKzone <- 3.5
botKzone <- 1.6
inKzone <- -0.95
outKzone <- 0.95
kZone <- data.frame(
x=c(inKzone, inKzone, outKzone, outKzone, inKzone),
y=c(botKzone, topKzone, topKzone, botKzone, botKzone)
)
df$h <- round(df$platelocheight)
df$s <- round(df$platelocside)
df$es<- round(df$exitspeed)
ggplot(kZone, aes(x,y)) +
stat_density_2d(data=df, aes(x=s, y=h, fill=es),geom="polygon") +
scale_fill_distiller(palette = "Spectral") +
geom_path(lwd=1.5, col="black") +
coord_fixed()
該代碼生成了一個帶有以下輸出的圖,但是您可以看到,它只是一種顏色,並且不像geom_tile圖那樣顯示exitspeed顏色輪廓。
我嘗試了各種不同的繪圖,包括geom_ratser和stat_contour來嘗試使這些繪圖起作用,但是我一直在努力尋找一種方法來使代碼正確地說明我的填充變量(exitspeed)。 該代碼要么運行並生成單色圖,要么出現錯誤:
Error in FUN(X[[i]], ...) : object 'exitspeed' not found
有誰知道如何正確更改geom_tile熱圖,以使輸出看起來更像stat_density_2D圖? 提前致謝!
我們可以建立一個新的網格並將您的數據插入到新的網格中。 這將使其看起來不像矩形y。
library(dplyr)
library(ggplot2)
library(gstat)
library(sp)
new_map <- df %>% rename(x = s, y = h)
coordinates(new_map) <- ~x + y
grd <- expand.grid(x = seq(from = -3, to = 3, by = .1), y = seq(from = 0, to = 5, by = .1))
coordinates(grd) <- ~x + y
gridded(grd) <- TRUE
idw <- idw(formula = exitspeed ~ 1, locations = new_map, newdata = grd)
idw.output <- as.data.frame(idw)
ggplot(kZone, aes(x,y)) +
geom_tile(data=idw.output, aes(x=x, y=y, fill=var1.pred)) +
scale_fill_gradientn(colours = rev(RColorBrewer::brewer.pal(10, "Spectral")), breaks = c(60, 70, 80, 90, 100), labels = c(60, 70, 80, 90, 100), limits = c(60,100))+
geom_path(lwd=1.5, col="black") +
labs(fill = "ExitSpeed")+
coord_fixed()
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.