Adjust local association measures' p-values
Make an adjustment to local association measures' p-values based on the number of neighbours (+1) of each region, rather than the total number of regions.
p.adjustSP(p, nb, method = "none")
p |
vector of p-values |
nb |
a list of neighbours of class |
method |
correction method as defined in |
A vector of corrected p-values using only the number of neighbours + 1.
Danlin Yu and Roger Bivand Roger.Bivand@nhh.no
data(afcon, package="spData") oid <- order(afcon$id) resG <- as.vector(localG(afcon$totcon, nb2listw(include.self(paper.nb)))) non <- format.pval(pnorm(2*(abs(resG)), lower.tail=FALSE), 2) bon <- format.pval(p.adjustSP(pnorm(2*(abs(resG)), lower.tail=FALSE), paper.nb, "bonferroni"), 2) tot <- format.pval(p.adjust(pnorm(2*(abs(resG)), lower.tail=FALSE), "bonferroni", n=length(resG)), 2) data.frame(resG, non, bon, tot, row.names=afcon$name)[oid,]
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.