首页 > 解决方案 > checkranin(tlim, tt, "tlim") 出错:'tlim[1]' 必须是 <'tlim[2]'

问题描述

所以我目前正在尝试执行一个时空核密度函数,我可以看到核密度分布随时间的变化。这是使用 sparr 包尝试的。我正在运行以下代码:


smell_Cases <- subset(newdata_proj, smell == '1',
                      select=c(x,y, smell))

smell_controls <- subset(newdata_proj, smell == '0',
                       select=c(x,y, smell))
smell_ppp <- list()

smell_ppp$cases<-ppp((smell_Cases$x), smell_Cases$y, marks=as_vector(as.integer(smell_Cases$smell)),
                     window=as.owin(as_Spatial(boundary)))

smell_ppp$controls<-ppp((smell_controls$x), smell_controls$y,
                        window=as.owin(as_Spatial(boundary)))

smell_ppp_Cases <- smell_ppp$cases


hlam <- LIK.spattemp(smell_ppp_Cases)

然后得到以下错误: Checkranin(tlim, tt, "tlim") 中的错误:'tlim[1]' must be < 'tlim[2]'

标签: rvectorspatialtemporalppp

解决方案


错误是您提供的数据的时间窗口无效。根据 help(LIK.spattemp) 上的文档(请参阅“pp”和“tt”参数的条目),如果您不提供每次观察的时间(您在上述函数调用中没有这样做)该函数将尝试使用数据对象的“标记”。您的数据对象的标记是观察时间吗?无论如何,我们需要 MWE 来为您提供全面帮助。


推荐阅读