Skip to content

Simplify MOFA2 pipeline with mofa2 wrapper - #871

Open
RiboRings wants to merge 6 commits into
develfrom
mofa2
Open

RiboRings wants to merge 6 commits into
develfrom
mofa2

Conversation

@RiboRings

Copy link
Copy Markdown
Member

With this PR, the chapter on MOFA2 is updated to reflect the streamlined MOFA2 pipeline using the mofa2 wrapper that was recently implemented with this PR: bioFAM/MOFA2#144.

@RiboRings
RiboRings requested a review from TuomasBorman July 28, 2026 12:37
@RiboRings

Copy link
Copy Markdown
Member Author

@TuomasBorman @antagomir This is ready to merge

@TuomasBorman TuomasBorman left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Looks good! Couple points (see comments)

Comment thread inst/pages/multiassay_ordination.qmd Outdated
Comment on lines +101 to 103
Note that agglomeration by taxonomic rank might produce duplicate feature names,
such as "uncultured". If present, they can be made unique as follows:

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Should this be handled by default, so that there should not be duplicated names?

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

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

Do you mean model mofa should return an error/warning if names are duplicated? Or making them unique inside the mofa function?

I think users would better off make the names unique before running mofa, otherwise the link between original names and unique names is lost and users cannot trace the results back to their original set of features.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Sorry, I checked this quickly.

I meant that

agglomerateByPrevalence() should output unique feature names by default

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

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

Oh ok, so this chunk can be safely removed, correct?

This branch has not been deployed

No deployments
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.

2 participants