[R-sig-ME] compare a neighbourhood value with the city mean - test of significance in presence of spatial autocorrelation

Thierry Onkelinx th|erry@onke||nx @end|ng |rom |nbo@be
Mon Jun 3 09:23:31 CEST 2019

Dear Dani,

What I would do is create a graph indicating which neighbourhoods are
directly connected to each other (share a border). And use this information
to fit a Besag-York-Mollier model which takes care of the auto-correlation
among the neighbourhoods. For an example see

You can find the INLA package at http://www.r-inla.org

Best regards,

ir. Thierry Onkelinx
Statisticus / Statistician

Vlaamse Overheid / Government of Flanders
Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance
thierry.onkelinx using inbo.be
Havenlaan 88 bus 73, 1000 Brussel

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey


Op ma 3 jun. 2019 om 01:50 schreef dani jones <orchidn using live.com>:

> Hi everyone,
> I have a dataset with all neighbourhoods in a city and a value for each
> neighbourhood (a rate). I would like to be able to create a vizualization
> that enables the user to select individual neighbourhoods and find out
> whether the value of that neighbourhood is significantly higher or lower
> than the city average across all neighbourhoods (I have a value that was
> calculated for the city). Given that the neighbourhoods are contiguous
> areas composing the city, there is spatial dependency, so I cannot use z
> scores or apply other classical statistics test since the observations are
> not independent and spatial autocorrelation might exist.
> I am not sure what test to apply to achieve this. Any help and links to
> potential readings would be very much appreciated.
> Thank you very much!
> Dani
> <http://aka.ms/weboutlook>
>         [[alternative HTML version deleted]]
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list