我想知道spatstat包中的marks命令是否可以用于分类变量而不是连续变量。我有一个xlsx表(A2N)有三列:大小(分类),CoordX(数字)和CoordY(数字)。
以下是我使用的脚本:
attach(A2N)
MPN2 <- ppp(CoordX, CoordY, c(621470,623141), c(9412006,9414455))
plot(MPN2)
marks(MPN2) <- A2N[, c("Size")]
plot(Smooth(MPN2))
这就是出现的信息:
> plot(Smooth(MPN2))
Error in r[i1] - r[-length(r):-(length(r) - lag + 1L)] :
non-numeric argument to binary operator
In addition: Warning messages:
1: some mark values are NA in the point pattern X
2: some mark values are NA in the point pattern X
3: some mark values are NA in the point pattern X