Question

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。

Answer 1

在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 - 本例中为内行）。

R将轮廓线保存到数据帧中

1 个答案: