Hello,
I am interested in using Morans I for different time intervals to detect
disease clusters.
Ultimately I would like to use CUSUM - or similar monitoring statistic to
monitor the results of Morans I - similar to the work by
Rogerson (2005) Spatial Surveillance and Cummulative Sum Methods in Spatial and
Syndromic Surveillance for Public Health.
Thus far - thanks to the list I have Morans I running in a loop for each day
(however I found that on some days no data is recorded, this caused an error to
get around this error I included an elseif statement to skip the calculation
for days where no disease notifications are recorded.)
I am just wondering if anyone has code/advice for the best way to apprach the
next stage of monitoring Morans I over different time intervals for the
detection of abnormal clusters. My code thus far is below - I would also like
to assemble the results in a neat table if anyone has tips for this I would alo
appreciate it as I am rather new to R!
Thanks
Serryn
T1<-read.dbf("ob.dbf")
N <- colnames(T1)
T2nb <- read.gal("PC.gal",override=TRUE)
col.W <- nb2listw(T2nb, style="W")
for(i in seq(N)) {
Vec1 <- spNamedVec(N[i+7], T1)
Svec1<-sum(Vec1)
ifelse ( Svec1> 0, res[i]<-moran.test(spNamedVec(N[i+7], T1), col.W,
zero.policy=TRUE), res[i]<-"NULL")
}
print(res)
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
