if (FALSE) {
# Load example data
data(ko_abundance)
data(metadata)
# Run integrated analysis
integrated_results <- ggpicrust2_extended(
data = ko_abundance,
metadata = metadata,
group = "Environment",
pathway = "KO",
daa_method = "LinDA",
ko_to_kegg = TRUE,
run_gsea = TRUE,
gsea_params = list(
method = "fgsea",
rank_method = "signal2noise",
nperm = 1000
)
)
# Access DAA results
daa_results <- integrated_results$daa_results
# Access GSEA results
gsea_results <- integrated_results$gsea_results
# Access plots
daa_plot <- integrated_results$daa_plot
gsea_plot <- integrated_results$gsea_plot
}
Run the code above in your browser using DataLab