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

时间:2016-11-08 21:53:32

标签: r

我正在使用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)

如何从chronplot函数打印用于创建平滑线的值?

1 个答案:

答案 0 :(得分:1)

这并不复杂,

我不知道为什么chron()函数没有给出平滑值。 但是,您可以使用同一包中的ffcsaps()函数来获取此平滑值。

示例:

使用ffcsaps()函数在数据框中添加一个包含平滑值的新列。

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)
告诉我它是否适合你 干杯