Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Python scripts to BYM2 case study #228

Draft
wants to merge 5 commits into
base: master
Choose a base branch
from

Conversation

farhanreynaldo
Copy link

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.

@farhanreynaldo farhanreynaldo marked this pull request as draft August 3, 2024 17:54
@bob-carpenter
Copy link
Contributor

Pinging @mitzimorris, who wrote the original ICAR case study.

If you want a dependency on a low-profile package, there's https://github.com/al-jshen/pyinla, but it looks like it just wraps the R version, so probably not a good idea.

@farhanreynaldo farhanreynaldo changed the title Python scripts to BYM12 case study Python scripts to BYM2 case study Aug 6, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants