我有一个NDVI观测的栅格堆栈，我想计算趋势的Theil斜率。我以前使用过bfast和greenbrown程序包，但是这是我第一次使用spatialEco程序包，并且出现以下错误：
my_stack <-raster.kendall（my_stack，tau = TRUE，截距= TRUE，p.value = TRUE，z.value = TRUE，置信度= TRUE）
错误（函数（y，x = seq（along = y），备选=“ two.side”，： 当ci.slope = TRUE时，必须至少有3个非缺失的有限观测值
我知道该错误与NA值的数量有关，但我没有完全理解。即使在某些情况下，资产净值很多，我也拥有多年来的数据。
有什么想法可以理解/解决此问题吗？