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

Figure legend in lefserPlotClad output #75

Open
shbrief opened this issue Oct 6, 2024 · 2 comments
Open

Figure legend in lefserPlotClad output #75

shbrief opened this issue Oct 6, 2024 · 2 comments
Assignees

Comments

@shbrief
Copy link
Contributor

shbrief commented Oct 6, 2024

Why are there only two circle sizes in the figure legend?

> library(lefser)
> ## Old cMD data in the lefser package
> data(zeller14)
> zeller14 <- zeller14[, zeller14$study_condition != "adenoma"]
> zeller14$study_condition <- factor(zeller14$study_condition,
+                                    levels = c("CRC", "control")) # re-levels
> tn <- get_terminal_nodes(rownames(zeller14))
> zeller14_tn <- zeller14[tn,]
> zeller14_tn_ra <- relativeAb(zeller14_tn)
> zeller14_input <- rowNames2RowData(zeller14_tn_ra)
> resAll <- lefserClades(zeller14_input, classCol = "study_condition")
18 features don't have species information.
96 features don't have strain information.
Dropped 96 features without full taxonomy information.
lefser will be run at the phylum, class, order, family, genus, species, strain level.

>>>> Running lefser at the phylum level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.

>>>> Running lefser at the class level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.

>>>> Running lefser at the order level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.

>>>> Running lefser at the family level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.

>>>> Running lefser at the genus level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.

>>>> Running lefser at the species level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.

>>>> Running lefser at the strain level. <<<<
The outcome variable is specified as 'study_condition' and the reference category is 'CRC'.
 See `?factor` or `?relevel` to change the reference category.
> lefserPlotClad(resAll) 
Using palette: colorblind

image

@shbrief
Copy link
Contributor Author

shbrief commented Oct 6, 2024

@sdgamboa Just double-checking: the circle size is 'LDA score' but not 'relative abundance'?

@sdgamboa
Copy link
Contributor

sdgamboa commented Oct 7, 2024

Yeah. It's the absolute LDA score. I think the scores are higher and more stable because the features abundance is agglomerated by clade.

@shbrief shbrief changed the title Figure legend in lefserClades output Figure legend in lefserPlotClad output Oct 8, 2024
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

2 participants