gpt4 book ai didi

从 geom_tile 中删除空行

转载 作者:行者123 更新时间:2023-12-04 10:28:10 27 4
gpt4 key购买 nike

我有一个包含 5 列的数据框:

N1 <- an integer between c(125,100,80,70,60,50,40,30,20)

N2 <- an integer between c(1,5,10,15,20,25,30,35,40,50,60,80,100)

Type <- Two different types

Rang <- a number

crit <- a character, only one value
N1总是大于 N2 .
这是我的数据框( Rang_final )开头的示例:
N1  N2     Type      rang     crit
125 1 SST-T_2m 41.86395 BE_proba
125 1 T_2m-SST 163.41217 BE_proba
100 1 SST-T_2m 32.88181 BE_proba
100 1 T_2m-SST 137.44479 BE_proba
80 1 SST-T_2m 22.57176 BE_proba
80 1 T_2m-SST 112.52334 BE_proba
70 1 SST-T_2m 21.30066 BE_proba
70 1 T_2m-SST 99.65523 BE_proba
60 1 SST-T_2m 18.48731 BE_proba
60 1 T_2m-SST 85.36945 BE_proba
50 1 SST-T_2m 18.60074 BE_proba
50 1 T_2m-SST 71.58960 BE_proba
40 1 SST-T_2m 18.58180 BE_proba
40 1 T_2m-SST 62.82670 BE_proba
30 1 SST-T_2m 20.53982 BE_proba
30 1 T_2m-SST 48.58923 BE_proba
20 1 SST-T_2m 27.15737 BE_proba
20 1 T_2m-SST 36.04175 BE_proba
125 5 SST-T_2m 101.82300 BE_proba
125 5 T_2m-SST 168.69954 BE_proba
100 5 SST-T_2m 91.10644 BE_proba
...

我只想在一张图中包含所有这些信息,所以我正在寻找使用 ggplot 和 geom_tile像这样:
p <- ggplot(Rang_final)
p <- (p
+ geom_tile(data=Rang_final[Rang_final$Type=="SST-T_2m",], aes(x=N1, y=N2, fill=rang))
+ geom_tile(data=Rang_final[Rang_final$Type=="T_2m-SST",], aes(x=N2, y=N1, fill=rang))
+ scale_fill_gradient2(name="Rang", low="deepskyblue",mid="yellow",high="red", midpoint=100, na.value = "grey50")
+ theme(axis.title.x = element_text(size=14, face="bold"),
axis.title.y = element_text(size=14, face="bold"),
strip.text.x = element_text(size=14, face="bold"),
strip.text.y = element_text(size=14, face="bold"),
axis.text=element_text(size=14),
axis.title.y=element_text(size=14, face="bold"),
legend.key=element_rect(size=0.5, colour="black"),
legend.text=element_text(size=10),
legend.margin=unit(0,"lines"),
legend.key.size=unit(0.8,"cm"),
legend.text.align=0)
+ theme_bw()

)

print(p)

这就是我得到的(N1 和 N2 的类是数字):
R

但我想删除空行,所以我尝试将 N1 和 N2 转换为因子:
Rang_final$N1 <- factor(Rang_final$N1, levels = c("20", "30", "40", "50", "60", "70", "80", "100", "125") ,ordered = TRUE)
Rang_final$N2 <- factor(Rang_final$N2, levels = c("1", "5", "10", "15", "20", "25", "30", "35", "40", "50", "60", "80", "100", "125") ,ordered = TRUE)

R

但这里我的数字不是按递增顺序排列的。

在这里,我也使用 N1 和 N2 作为因子,但我无法在同一个图表上以这种比例显示两种不同的类型:
p <- ggplot(Rang_final, aes(x=N1, y=N2))
p <- (p
+ geom_tile(aes(fill=rang))
+ facet_grid(Type ~ crit, scales="free")
+ scale_fill_gradient2(name="Rang", low="deepskyblue",mid="yellow",high="red", midpoint=100, na.value = "grey50")
+ theme(axis.title.x = element_text(size=14, face="bold"),
axis.title.y = element_text(size=14, face="bold"),
strip.text.x = element_text(size=14, face="bold"),
strip.text.y = element_text(size=14, face="bold"),
axis.text=element_text(size=14),
axis.title.y=element_text(size=14, face="bold"),
legend.key=element_rect(size=0.5, colour="black"),
legend.text=element_text(size=10),
legend.margin=unit(0,"lines"),
legend.key.size=unit(0.8,"cm"),
legend.text.align=0)
+ theme_bw()

)

R

有人能帮我吗?

最佳答案

这是我的建议。当您将 N1 和 N2 转换为因子时,您的 x 轴和 y 轴编号不按顺序排列的原因是您没有分配所需的顺序。这里我使用了 overflow包,因此您的数据框称为 mydf,与您发布的内容相同。我认为你需要做两件事。

1) 重新排序因子

您需要的是以正确的方式重新排列您的因子。当您在脚本中将字符转换为因子时,会发生类似的事情。

#> mydf$N1 <- as.factor(mydf$N1)
#> levels(mydf$N1)
#[1] "20" "30" "40" "50" "60" "70" "80" "100" "125"
#> unclass(mydf$N1)
#[1] 9 9 8 8 7 7 6 6 5 5 4 4 3 3 2 2 1 1 9 9
#attr(,"levels")
#[1] "20" "30" "40" "50" "60" "70" "80" "100" "125"

如果你看 unclass() ,您会看到“20”有 9。您希望“20”在 unclass() 中有 1 .因此,您希望对数据框进行排序。这导致了奇怪的数字顺序。在将字符转换为因子之前,您希望对 N1 和 N2 进行排序,并将它们转换为因子。
library(dplyr)
# Sort mydf by N1
mydf <- arrange(mydf, N1)

# Convert N1 to factor
mydf$N1 <- as.factor(mydf$N1)

在这里你看到“20”在 unclass 中有 1 个,而“125”有 9 个。所以这是你想要的顺序。
#> levels(mydf$N1)
#[1] "20" "30" "40" "50" "60" "70" "80" "100" "125"
#> unclass(mydf$N1)
#[1] 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 9 9
#attr(,"levels")
#[1] "20" "30" "40" "50" "60" "70" "80" "100" "125"

2) 数据子集

这可能/可能没有必要。但是,到目前为止,我认为这是必要的。以下似乎不太好。
  + geom_tile(data=Rang_final[Rang_final$Type=="SST-T_2m",], aes(x=N1, y=N2, fill=rang))
+ geom_tile(data=Rang_final[Rang_final$Type=="T_2m-SST",], aes(x=N2, y=N1, fill=rang))

所以,我做了一个技巧来安排你的数据。
 # Subset data using dplyr
america <- filter(mydf, Type == "SST-T_2m")
brazil <- filter(mydf, Type == "T_2m-SST")

# Reverse N1 and N2 in brazil (So N1 is actually N2, and N2 is actually N1)
colnames(brazil) <- c("N2","N1", "Type", "rang", "crit")
brazil[,c(2,1,3:5)]

在对原始数据进行子集化并对巴西数据框进行处理后,我将它们组合起来。然后,我按照前面提到的方式安排因子顺序。
 canada <- rbind(america, brazil)
canada <- arrange(canada,N1)
canada$N1 <- as.factor(canada$N1)

canada <- arrange(canada,N2)
canada$N2 <- as.factor(canada$N2)

现在我画一个图。我不知何故对传说相关的事情有问题。所以我在这里删除了它们。
x <- ggplot(canada)
x <- (x
+ geom_tile(aes(x=N1, y=N2, fill =rang))
+ scale_fill_gradient2(name="Rang", low="deepskyblue",mid="yellow",high="red", midpoint=100, na.value = "grey50")
+ theme(axis.title.x = element_text(size=14, face="bold"),
axis.title.y = element_text(size=14, face="bold"),
strip.text.x = element_text(size=14, face="bold"),
strip.text.y = element_text(size=14, face="bold"),
axis.text=element_text(size=14),
axis.title.y=element_text(size=14, face="bold"))
+ theme_bw()

)

让我知道这是否适合您。

enter image description here

关于从 geom_tile 中删除空行,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/25603624/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com