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

run_aldex enrich_group source code question #118

Open
possumhound opened this issue Jul 24, 2024 · 0 comments
Open

run_aldex enrich_group source code question #118

possumhound opened this issue Jul 24, 2024 · 0 comments

Comments

@possumhound
Copy link

Hi,

first of all thank you for putting all these microbiome functions into a digestible format - it's great!

However, I am a little confused by the source code below in the run_aldex function:

enrich group

cds <- gsub("rab.win.", "", names(ef)[2:3])
ef <- ef$effect
enrich_group <- ifelse(ef > 0, cds[1], cds[2])

should it not be:

enrich_group <- ifelse(ef > 0, cds[2], cds[1]) ?

Maybe I am misunderstanding something but the way it is currently written it contradicts my abundance plots and my ancombc2 results.

Any input would be much appreciated!

Thank you!

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

No branches or pull requests

1 participant