BFAST包的警告错误和季节性问题

Question

关于x86_64-apple-darwin13.4.0的R版本3.3.2（2016-10-31）

我正在尝试将每日分辨率的11年时间序列输入BFAST模型，以查看它是否能够检测出符合我们对输入数据更改日期的知识的中断。该系列本身实际上是针对北极的算法衍生的海冰范围产品与可操作的（人类绘制的）海冰范围产品之间的差异。具体地，差异系列是（人类产品） - （算法产品），使得所有正值表示人类产品具有更多冰，并且相反地表示算法产品。有效数据范围为+/- 1 * 10 ^ 6平方公里，如下所示：

以下是差异系列第一年的前90天数据：

Data <- c(-175276.05, -232276.05, -136336.76, 95715.53, 151779.12, 
164566.44, 218769.37, 200736.84, 105887.64, 14426.35, -41969.07, 
69030.93, 73349.8, -27075.22, -142136.15, -136136.15, 1863.85, 
-16761.15, 100394.9, 113368.56, 61704.79, -16002.55, 122554.58, 
-72809.77, -184809.77, -247407.09, -173722.69, -91140.5, -34049.69, 
-158304.83, -175901.74, -211442.38, -180052.53, -174066.45, -70030.18, 
-82030.18, -33284.22, -4265.03, -41862.19, -30639.13, -62305.53, 
-125305.53, -176305.53, -118305.53, 41406.63, 608.19, -100725.02, 
-76505.16, 7763.92, -34151.19, -131002.89, -110096.16, -72502.6, 
80861.21, 58577.95, -28883.35, -203696.82, -273928.14, -285928.14, 
-322928.14, -107873.1, 19767.96, -84003.32, -89605.11, -68633.45, 
-160704.55, -170663.39, -98237.86, -113962.63, 41085.45, 98463.3, 
87798.23, 148593.44, 90224.71, 84678.51, -3096.93, -512.43, -56539.26, 
-51659.94, -5381.75, 155551.74, -35037.1, -118243.32, -155905.59, 
-154578.81, -343272.84, -363430.58, -316900.35, -345670.7, -306597.62)

如果我使用以下代码将这些数据提供给bfast包（v1.5.7）：

bfast.fit <- bfast(data, h = 0.30, season = "none", max.iter = 2, breaks = NULL, level = 0.10, type = "OLS-MOSUM")

并且情节，我确实得到了结果。但是，有一些“警告”消息被抛出而没有进一步描述：

There were 50 or more warnings (use warnings() to see the first 50)
warnings()
Warning messages:
1: In sqrt(fr) : NaNs produced
2: In sqrt(fr) : NaNs produced
3: In sqrt(fr) : NaNs produced

对这些警告有进一步的解释吗？

其次，我不确定是否有任何季节性模型适用于此数据集。通过傅里叶分析，我知道这个时间序列在335,185和117天有多个季节性。

虚拟或谐波季节对于这种数据集是否有意义？