running anova with multigroup analysis

[rbutleriii]

Hello, I am trying a four group analysis, and wanted to pass anova to the plot_context_boxplot function. I get the error:

Error in `mutate()`:
ℹ In argument: `test = map(...)`.
Caused by error in `map()`:
ℹ In index: 1.
Caused by error in `UseMethod()`:
! no applicable method for 'anova' applied to an object of class "formula"
Backtrace:
     ▆
  1. ├─liana::plot_context_boxplot(...)
  2. │ └─... %>% ...
  3. ├─dplyr::mutate(., across(where(is.numeric), ~round(.x, digits = 4)))
  4. ├─dplyr::mutate(., p.adj = p.adjust(p.value, method = "fdr"))
  5. ├─tidyr::unnest(., test)
  6. ├─dplyr::mutate(...)
  7. ├─dplyr:::mutate.data.frame(...)
  8. │ └─dplyr:::mutate_cols(.data, dplyr_quosures(...), by)
  9. │   ├─base::withCallingHandlers(...)
 10. │   └─dplyr:::mutate_col(dots[[i]], data, mask, new_columns)
 11. │     └─mask$eval_all_mutate(quo)
 12. │       └─dplyr (local) eval()
 13. ├─purrr::map(...)
 14. │ └─purrr:::map_("list", .x, .f, ..., .progress = .progress)
 15. │   ├─purrr:::with_indexed_errors(...)
 16. │   │ └─base::withCallingHandlers(...)
 17. │   ├─purrr:::call_with_cleanup(...)
 18. │   └─liana (local) .f(.x[[i]], ...)
 19. │     ├─broom::tidy(...)
 20. │     ├─rlang::exec(...)
 21. │     └─stats::anova(...)
 22. └─base::.handleSimpleError(...)
 23.   └─purrr (local) h(simpleError(msg, call))
 24.     └─cli::cli_abort(...)
 25.       └─rlang::abort(...)
Execution halted

Is this because multiple-group testing hasn't been implemented?

[rbutleriii]

Update, when running with aov instead it works. The warning messages seem like it is because I have paired = FALSE included in the call which is fine.

Warning message:
There were 10 warnings in `mutate()`.
The first warning was:
ℹ In argument: `test = map(...)`.
Caused by warning:
! In lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
 extra argument ‘paired’ will be disregarded
ℹ Run `dplyr::last_dplyr_warnings()` to see the 9 remaining warnings. 

One notable thing though is that each factor is now plotted twice (second padj is NA). is that just trying to plot padj for the residuals? Thoughts on correcting?