我正在使用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
函数打印用于创建平滑线的值?
答案 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)
告诉我它是否适合你
干杯