R: dplR，chron，提取平滑的年表

Question

我正在使用detrend包中的dplR函数，使用负指数/负线性方法来降低树环数据库的性能。

dt.rwi <- detrend(rwl=d, verbose=TRUE, method="ModNegExp", pos.slope = FALSE)

然后，我根据每个系列中的去趋势值构建一个年表，并使用nyrs=30命令平滑。

neg.crn <- chron(dt.rwi, biweight = TRUE, prewhiten = FALSE, nyrs=30)

问题是，当我打印neg.crn对象时，我只得到了每年的平均去趋势指数，而不是平滑的值。该绘图返回平滑的行，但同样，我无法提取用于构造这一行的数据。

plot(neg.crn, add.spline=T, nyrs=30)

如何打印用于从chron plot 或plot函数创建平滑行的值？

Answer 1

这并不复杂，

我不知道chron()函数为什么不给出平滑的值。但是，您可以使用来自同一个包的ffcsaps()函数来获得这个平滑的值。

示例：

使用ffcsaps()函数在dataframe中添加一个带有平滑值的新列。

neg.crn$smooth32 <- ffcsaps(neg.crn[,1], nyrs = 32) # that is all! Now just check these values in your dataframe.

测试它绘制你的数据

plot(neg.crn[,1], type = "l")
lines(neg.crn$smooth32, col = "red", lwd = 2)

告诉我这对你有用吗，干杯