ggplot 允许我为数据系列 x 和 y 绘制等高线:
library(ggplot2)
df <- data.frame(x=c(1:100),y=rnorm(100))
ggplot(df,aes(x=x,y=y))+geom_density2d()
我希望能够将主要四个等高线 (95,75,50,25)max 和 (95,75,50,25)min 的数据输出保存在数据框中。我会很感激你的帮助。也许有一种方法可以直接执行此操作,而无需使用 ggplot。
在对象周围使用 function ggplot_build(
),ggplot()
您可以访问用于绘图的所有数据。数据存储在列表元素中data
。
p<-ggplot_build(ggplot(df,aes(x=x,y=y))+geom_density2d())
str(p$data)
List of 1
$ :'data.frame': 1895 obs. of 6 variables:
..$ level: num [1:1895] 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 5e-04 ...
..$ x : num [1:1895] 1 2 3 4 5 ...
..$ y : num [1:1895] 1.42 1.44 1.45 1.45 1.46 ...
..$ piece: int [1:1895] 1 1 1 1 1 1 1 1 1 1 ...
..$ group: Factor w/ 13 levels "1-001","1-002",..: 1 1 1 1 1 1 1 1 1 1 ...
..$ PANEL: int [1:1895] 1 1 1 1 1 1 1 1 1 1 ...
您可以将这些数据存储为单独的数据框。
gg<-p$data[[1]]
head(gg)
level x y piece group PANEL
1 5e-04 1.000000 1.423926 1 1-001 1
2 5e-04 2.000000 1.435286 1 1-001 1
3 5e-04 3.000000 1.445293 1 1-001 1
4 5e-04 4.000000 1.454166 1 1-001 1
5 5e-04 5.000000 1.462106 1 1-001 1
6 5e-04 5.255343 1.463967 1 1-001 1
group
变量显示每个值属于哪一行(在这种情况下,piece
1 - 外线,13 - 内线)。