-
Notifications
You must be signed in to change notification settings - Fork 74
Represent statespace metadata with dataclasses #607
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
base: main
Are you sure you want to change the base?
Represent statespace metadata with dataclasses #607
Conversation
|
Check out this pull request on See visual diffs & provide feedback on Jupyter Notebooks. Powered by ReviewNB |
jessegrabowski
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This is a great first pass, much cleaner than what we have now.
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
|
We can also keep all of the existing properties like |
|
Reflecting on it, I am convinced this is the way to go. It's 1000x more ergonomic. I made some changes to your initial code to make the API more "dictionary like", and to reduce code duplication. I moved everything to |
|
@jessegrabowski, this is looking really cool! What can I do to help push this forward? |
|
Delete the new We should keep your notebook with the plan to add it as a new example for the docs. Or it can be merged into the custom statespace notebook. So that should also be updated to import from the new |
|
Perfect! I'll work on that today!! It is really looking cool! |
pymc_extras/statespace/models/structural/components/regression.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression_dataclass.py
Outdated
Show resolved
Hide resolved
tests/statespace/models/structural/components/test_regression.py
Outdated
Show resolved
Hide resolved
|
@jessegrabowski, I agree with all of your comments above. I am going to start making those changes. |
…uplicate with warning 2. removed unnecessary imports from __init__ after deleting regression_dataclass 3. updated components and structural classes to only utilize dataclasses and pull other objects from <foo>_info dataclasses 4. updated tests to conform to dataclass api
2. created tests for add and merge methods 3. added utility to convert from snake to pascal and integrated it in error messaging
… and placed default shoch and state setters
jessegrabowski
left a comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Incomplete review, I'll continue tomorrow AM
| # if key in index: | ||
| # raise ValueError(f"Duplicate {self.key_field} '{key}' detected.") # This needs to be possible for shared states |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
That shouldn't happen here though, it should come up in merge or add right? And we handle it there with the allow_duplicates flag
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think what happens is because our data classes are immutable the __post_init__ runs right after our merge/add because we always return new objects of the same dataclass and it see that there are duplicate keys even though the merge/add method had allowed them via allow_duplicates.
| raise AttributeError(f"Items missing attribute '{self.key_field}': {missing_attr}") | ||
| object.__setattr__(self, "_index", index) | ||
|
|
||
| def _key(self, item: T) -> str: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is this used?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Doesn't look like we are using this at the moment. I think it could come in handle later on to check for specific attributes.
pymc_extras/statespace/models/structural/components/regression.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression.py
Outdated
Show resolved
Hide resolved
pymc_extras/statespace/models/structural/components/regression.py
Outdated
Show resolved
Hide resolved
| self.state_names = list(state_names) if state_names is not None else [] | ||
| self.observed_state_names = ( | ||
| list(observed_state_names) if observed_state_names is not None else [] | ||
| self.param_info = ParameterInfo( |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Should we change the component signature to just take the Info objects directly?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think that would be nice. Right now, there is this intermediate conversion step that needs to take place.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Ok let's leave it for the next PR then and make the priority here just getting each model to be represented in the new way
…_duplicates is False 2. converted component attributes into properties 3. removed _combine_property method 4. removed redundant observed_states property 5. fixed indentation bug
|
Hey @jessegrabowski, by switching a lot of the component attributes to properties I was able to simplify a good amount of downstream methods. If you don't mind taking a look at the current state of this before I go ahead and do the same with the rest of the SSM components. |
|
|
||
| self.coords_info = CoordInfo(coords=[regression_state_coord, endogenous_state_coord]) | ||
|
|
||
| def populate_component_properties(self) -> None: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This method won't be unique to regression right? We will want to move it up to the base class.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
@jessegrabowski, in the base class there is a populate_component_properties method that raises a NotImplemented. Did you want to replace that with a generic method that sets _set_<foo> for the 2 defaults (shocks and states) that we provide?
| self.make_symbolic_graph() | ||
|
|
||
| self._component_info = { | ||
| self._component_info = { # Should this be a dataclass?? |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Probably :)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Let's leave it for now though. The mission creep on this PR is already pretty terrible.
| raise NotImplementedError | ||
|
|
||
| def _set_shocks(self) -> None: | ||
| def _set_shocks(self) -> ShockInfo: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Make a comment in each of these methods that these are generic defaults that can/should be replaced with component specific logic. I got momentarily confused because they are called _set but they actually return info (i.e. they do a get, not a set). The name makes sense in the context of populate_component_properties`, but not in a vacuum; hence the comment.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Oh, that is right. These shouldn't be doing a get they should actually be doing a set. I will fix it and put in a docstring with the explanation that these are defaults that should be replaced.
|
Yeah it looks really great! Go ahead and do the others. Excited to get this over the finish line |
This is a draft proposal for #598
The idea is to handle each component separately using
_set_{component}methods and all information are stored using data classes for easy mapping.I believe this will simplify our tests of these components and will reduce redundancies where we have the same information spread across multiple sub-components like
data_namesanddata_info.@jessegrabowski let me know what you think I put a little notebook together to showcase the changes.