0

我正在使用spdeplocalmoran,它根据 p.adjustSP 调整 p 值以进行多重比较。我想使用 holm 校正,但我不确定是否有问题,因为我注意到无论我指定 holm、bonferroni 还是 hochberg,它都会导致相同的调整 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])
4

0 回答 0