计算x值

Question

data <- dput(data): structure(list(x = 1:16, y = c(-79.62962963, -84.72222222, -88.42592593, -74.07407407, -29.62962963, 51.38888889, 79.62962963, 96.2962963, 87.96296296, 88.42592593, 73.14814815, 12.96296296, -63.42592593, -87.03703704, -87.5, -87.96296296)), .Names = c("x", "y"), row.names = c(NA, 16L), class = "data.frame") 

我中的R计算与GGPLOT2平滑线为我的数据集：

p1 <- ggplot(data, aes(x=x(°), y=(%))) 

library(splines) 
library(MASS) 
(p2 <- p1 + stat_smooth(method = "lm", formula = y ~ ns(x,3)) + 
    geom_point() 
) 

如何可以计算的x值平滑线的曲线最大值？

Answer 1

你需要做一些数学上的计算stat_smooth产生的data.frame：

library(ggplot2) 
library(splines) 
library(MASS) 

data <- structure(list(x = 1:16, 
         y = c(-79.62962963, -84.72222222, -88.42592593, 
          -74.07407407, -29.62962963, 51.38888889, 
          79.62962963, 96.2962963, 87.96296296, 
          88.42592593, 73.14814815, 12.96296296, 
          -63.42592593, -87.03703704, -87.5, 
          -87.96296296)), .Names = c("x", "y"), 
        row.names = c(NA, 16L), class = "data.frame") 

p1 <- ggplot(data, aes(x=x, y=y)) 
p1 <- p1 + stat_smooth(method = "lm", formula = y ~ ns(x,3)) 
p1 <- p1 + geom_point() 
p1 

gb <- ggplot_build(p1) 

exact_x_value_of_the_curve_maximum <- gb$data[[1]]$x[which(diff(sign(diff(gb$data[[1]]$y)))==-2)+1] 

p1 + geom_vline(xintercept=exact_x_value_of_the_curve_maximum) 

exact_x_value_of_the_curve_maximum 
[1] 9.164557 

有more robust ways为好，但你仍然需要ggplot_build部分来获取数据。