I have locally implemented additional functionality to the joint modeling code in the rstanarm package in order to appropriately model negative binomial data in the longitudinal component as part of my doctoral thesis work. I am trying to determine whether it would be better to incorporate these changes into the existing rstanarm package or to create my own R package. I have made the following changes:
- Enabled offset values in the longitudinal component
- Introduced a scaling parameter for the association term in the survival component
- Implemented the muslope association structure
Would these changes would be appropriate for the rstanarm package?