でポイント パターン モデルのレバレッジ プロットを作成しようとして、問題が発生していspatstatます。ただし、これは奇妙です。昨日はまったく同じモデルで完全に正常に機能していたのに、現在は機能していないからです。
3 つの共変量と Geyer 相互作用を持つモデルがあり、警告メッセージを生成せずに正しく実行されます。ただし、レバレッジ プロットを生成しようとすると、最初に次の警告が表示されます。
Large quadrature scheme split into blocks to avoid memory size
limits;6404 dummy points split into 8 blocks, the first 7 blocks
containing 801 dummy points and the last block containing 797 dummy
points
そして、このエラー：
Error in ppmInfluenceEngine(fit, what = what, ..., iScore = iScore, iHessian = iHessian, :
number of items to replace is not a multiple of replacement length
この問題は、R が直交スキームを分割しようとするという事実と関係があると思います。なぜなら、以前はクラッシュせず、まったく同じモデルでレバレッジ プロットを作成できたので、そうする必要がなかったからです。それを分割します（または、少なくともそれについて私に知らせません！）
library(spatstat)
#generate a pattern
pattern <- rpoispp(1454)
summary(pattern)
#generate a covariate
cov1 <- rnoise(rgen = rnorm, dimyx=32, mean=2, sd=1, w = pattern$window)
cov2 <- rnoise(rgen = rnorm, dimyx=32, mean=1, sd=7, w = pattern$window)
cov3 <- rnoise(rgen = rnorm, dimyx=32, mean=5, sd=10, w = pattern$window)
#fit the model
fit <- ppm(pattern ~ cov1 + cov2 + cov3, Geyer(r = 0.01, s = 0.02))
leverage(fit) #this produces the warning and then the error
R バージョン 3.3.1 と spatstat バージョン 1.46-1 を実行しています