如何解决错误因素在 R 中的错误级别

Question

我在R非均匀 G 函数应用于我的点模式时遇到了一些困难。 为了使用GmultiInhom ，我首先尝试将我的点模式bci.tree8pppa转换为多类型模式：

bci.tree8multi = ppp(bci.tree8pppa$x, bci.tree8pppa$y, window=owin(c(0,1000), c(0,500)), marks = factor(bci.tree8pppa$marks[,3]))

然后应用G函数如下：

G = GmultiInhom(bci.tree8multi, marks(bci.tree8multi) == species1, marks(bci.tree8multi) == species2, lambdaI = lambda1points, lambdaJ = lambda2points, lambdamin = min(lambda2points), r = c(0,r1,r2,r3))

但这会产生错误： "Error in split.default(X, group): factor has bad level"

我该如何解决这个问题？ 先感谢您！

Answer 1

为了任何R程序员的利益：我在基础R系统中的.Internal(split.default(x,f))的 C 源代码中追踪了错误消息"factor has bad level" 。 仅当f是list而不是factor时，才会出现此错误消息。 该代码使用执行字符串操作的 function interaction将f转换为因子。 这种转换可能 go 是错误的，因为结果因子具有“不良水平”：因子的 integer 表示包括小于 1 或大于因子水平数的值。 然后发生错误。

原始帖子没有提供一个工作示例，因此很难弄清楚spatstat::GmultiInhom中的代码是如何导致将错误类型的数据提供给split.default 。 但是，它必须与参数r的滥用有关。 spatstat中的代码将被收紧，以对r的格式执行更严格的要求。