example-models
example-models copied to clipboard
Python scripts to BYM2 case study
The following PR is a draft on Fitting the Model to the Scotland Lip Cancer Dataset
and Bigger data: from 56 counties in Scotland to 1921 census tracts in New York City
sections.
There are several issues need to be addressed.
- The BYM2 model cannot be directly translated from R to Python because as far as I know, there is no equivalent INLA library in Python.
- The raw number of 2001 events plot has different value scale compared to the R implementation. The Python plot has a range of value from 0 to 300, while the R plot has 0 to 15.
- There is an exception for BYM model in Python, however the summary still produced.