私はspdep の localmoranを使用しています。これは、 p.adjustSPに基づいて複数の比較のために p 値を調整します。ホルム補正を使用したいのですが、ホルム、ボンフェローニ、ホッホベルグのいずれを指定しても、すべて同じ調整済み p 値になることに気付いたので、何か問題があるかどうかはわかりません。補正を指定しない場合、補正されていない p 値は問題ないように見えます。他の誰かがこれに遭遇しましたか?それは何かが間違っていることを意味しますか? それとも、これらすべての補正が同じ値になる可能性はありますか? ありがとう!
例を追加するために編集しました(私の実際のデータではありません):
library(spdep)
dat <- st_read(system.file("etc/shapes/bhicv.shp", package="spdep")[1])
# create neighborhood
neighbourhood <- poly2nb(dat, queen=TRUE)
neighbourhood_weights_list <- nb2listw(neighbourhood,style="W", zero.policy=TRUE)
# add unadjusted p values
dat$lmoran_pNONE <- localmoran(dat$ELCI,neighbourhood_weights_list,p.adjust.method="none",
na.action=na.exclude,zero.policy=TRUE)[,5]
# add adjusted p values
dat$lmoran_pHOLM <- p.adjustSP(dat$lmoran_pNONE, neighbourhood, method = "holm")
dat$lmoran_pBONF <- p.adjustSP(dat$lmoran_pNONE, neighbourhood, method = "bonferroni")
# show p values
head(dat[10:12])