R の bfast 関数の結果について質問があります。ブレークが検出されない時系列があるとします。
library(bfast)
library(zoo)
NDVI <- as.ts(zoo(som$NDVI.b,som$Time))
NDVI_w <- window(NDVI, c(2001, 4) , c(2008, 13))
fit <- bfast(NDVI_w, h=1/2, season="dummy", max.iter=1)
plot(fit)
bfast (フィット) オブジェクトからトレンド コンポーネントの勾配値を抽出するにはどうすればよいですか? ANOVA=TRUE のオプション 'plot' は機能しません。
plot(fit, ANOVA=TRUE)$slope
1 つのオプションは、提供されたトレンド コンポーネントに基づいて計算することですが、「フィット」オブジェクトから直接取得する方法はありますか?
out <- fit$output[[1]]
plot(out$Tt)
lm(out$Tt ~ time(out$Tt))$coefficients[2]
どんなヒントでも大歓迎です。