如何提取ggplot平滑函数并保存为小时级数据框？

我想要从ggplot的平滑曲线中提取对应的函数，并将其保存为包含小时级数据点的数据框。

已尝试方案 #

我试过多种插值方法，但效果都不理想：

• 线性插值会产生锯齿状图案；
• na_spline会生成怪异的曲线图案。

实际数据的趋势更接近ggplot中geom_smooth的效果，于是尝试用以下代码复现：

loess.data <- stats::loess(Hallwil2018_2019$Avgstemp~as.numeric(Hallwil2018_2019$datetime), span = 0.5)
loess.predict <- predict(loess.data, se = T)

但这段代码生成的结果列表缺失NA值，且长度远小于预期。

内容的提问来源于stack exchange，提问作者Cla FM