Change levels for LM modeling in tax_tree

[scilexenko]

Hi David,

I hope you’re doing well!

I’m encountering an issue with my data, which consists of four levels. I’m trying to change the comparison level to generate different heat trees. I attempted to relevel the factor within my phyloseq object using the following code:

sample_data(phyloseq_g_zz)$Code <- factor(sample_data(phyloseq_g_zz)$Code, levels = c("SNA", "CNA", "SMA", "SVA"))

lm_models2 <- phyloseq_g_zz  %>% 
  taxatree_models(
    ranks = c("Phylum", "Class", "Order", "Family", "Genus", "Species"),
    trans = "log2", trans_args = list(add = "halfmin"),
    formula = ~ Code,
    #variables = "Code", type = lm # modelling function
  ) %>%
  taxatree_models2stats(.keep_models = F) 

However, this doesn’t seem to work as expected when using taxatree_models.
Would you be able to help me with this? I would greatly appreciate any guidance you can provide.

Thank you so much in advance!
Best regards,
Ivan