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

Paired heavy-light modeling #92

Open
wants to merge 13 commits into
base: main
Choose a base branch
from
Open

Paired heavy-light modeling #92

wants to merge 13 commits into from

Conversation

matsen
Copy link
Contributor

@matsen matsen commented Dec 10, 2024

No description provided.

@matsen matsen changed the title a sketch Paired heavy-light modeling Dec 10, 2024
@matsen matsen linked an issue Dec 10, 2024 that may be closed by this pull request
@willdumm
Copy link
Contributor

WIP

WIP

next: fix nonzero subs prob for ^

WIP

WIP

remove joined_mode, make everything be joined_mode always

EOD commit: working on recovering OE plotting with changes

the final solution?

partial cleanup

some cleanup

more cleanup
Copy link
Contributor Author

@matsen matsen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looking nice. Just a few suggested clarifications. I wonder if the way you are using "token" is more specialized than it appears.

netam/sequences.py Outdated Show resolved Hide resolved
]
STOP_CODONS = ["TAA", "TAG", "TGA"]
# Each token in RESERVED_TOKENS will appear once in aa strings, and three times
# in nt strings.
TOKEN_TRANSLATIONS = {token * 3: token for token in RESERVED_TOKENS}
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

So these are actually only reserved token translations?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Prepended RESERVED_

netam/sequences.py Outdated Show resolved Hide resolved
netam/framework.py Show resolved Hide resolved
Copy link
Contributor Author

@matsen matsen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Seems like as elegant as one could ask for given the setting 👍

netam/dasm.py Show resolved Hide resolved
netam/dxsm.py Show resolved Hide resolved
netam/dxsm.py Outdated
# # The following can be used when one wants a better traceback.
# burrito = self.__class__(None, dataset, copy.deepcopy(self.model))
# return burrito.serial_find_optimal_branch_lengths(dataset, **optimization_kwargs)
# TODO
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Noting todo here

if chosen_v_families is not None:
chosen_v_families = set(chosen_v_families)
pcp_df = pcp_df[pcp_df["v_family"].isin(chosen_v_families)]
# TODO Does this seem like the right thing to do?
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think we want AND here, not or. We want the H and L to be restricted to the specified set.

@matsen
Copy link
Contributor Author

matsen commented Dec 28, 2024

(BTW, because I initiated this PR, I can't approve it 😁 )

@willdumm willdumm marked this pull request as ready for review December 30, 2024 17:09
@matsen
Copy link
Contributor Author

matsen commented Dec 30, 2024

Great! @willdumm if you approve then you should be able to squash n merge.

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.

Do paired heavy-light modeling
2 participants