Error in ancombc about ancombc HOT 3 CLOSED

Hi @jakob-wirbel,

Thank you very much for catching this! The package will be updated shortly to fix this BUG!

Best,
Huang

from ancombc.

Hi @lakerwsl,

Thanks for posting your question! In order for me to better dig into the problem, may you please share with your codes and a mock dataset (with all confidential information marked)?

Thank you,
Huang

from ancombc.

Hey hey 👋
i ran into the same problem and i think i found the reason for the error 😄

If you have a single metadata variable (basically, only the group label) AND you would remove samples with the minimum library size cutoff, this part of the code will return the metadata as a vector, not as a data.frame:

# Discard samples with library size < lib_cut
lib_size = colSums(feature_table, na.rm = TRUE)
if(any(lib_size < lib_cut)){
   subj_del = which(lib_size < lib_cut)
   feature_table = feature_table[, - subj_del]
   meta_data = meta_data[- subj_del, ]
}

I think it could be fixed very quickly by changing the line to:

meta_data = meta_data[- subj_del, , drop=FALSE]

Cheers

from ancombc.