-3

RRをプロットするRコードは次のとおりです。

# Display posterior means of RR(relative risks)
va$RRmean <- end4$summary.fitted.values[, 1]
library(classInt)
breaks.qt2 <-classIntervals(va$RRmean, n=7, style="quantile", intervalClosure="right")
va.palette<-brewer.pal(9, name="OrRd")
spplot(va, "RRmean", col="transparent", col.regions=va.palette, at=breaks.qt2$brks)

ただし、空白の白として表示されている 1 つの郡があります (ただし、他の郡の中で最も高い値を持っています)。コードに何か問題がありますか?実際の値を確認しましたが、有効な数値です。

4

1 に答える 1