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

Question about small sample size #284

Open
UniAlberta opened this issue Sep 11, 2024 · 2 comments
Open

Question about small sample size #284

UniAlberta opened this issue Sep 11, 2024 · 2 comments
Labels
invalid This doesn't seem right

Comments

@UniAlberta
Copy link

I was wondering to ask about the sample size in each level of a categorical group for analyzing differentially abundant species. I only have 3 samples for each level. I can run ANCOMBC function as below. At the end, it says "You may have unstable results because of sample size <5". I wanted to know if my results are still reliable or if I need to look for another package.
output = ancombc2(data = tse, assay_name = "counts", tax_level = "Family",
fix_formula = "age + region + bmi", rand_formula = NULL,
p_adj_method = "holm", pseudo_sens = TRUE,
prv_cut = 0.10, lib_cut = 1000, s0_perc = 0.05,
group = "bmi", struc_zero = TRUE, neg_lb = TRUE,
alpha = 0.05, n_cl = 2, verbose = TRUE,
global = TRUE, pairwise = TRUE, dunnet = TRUE, trend = TRUE,
iter_control = list(tol = 1e-2, max_iter = 20,
verbose = TRUE),
em_control = list(tol = 1e-5, max_iter = 100),
lme_control = lme4::lmerControl(),
mdfdr_control = list(fwer_ctrl_method = "holm", B = 100),
trend_control = list(contrast = list(matrix(c(1, 0, -1, 1),
nrow = 2,
byrow = TRUE),
matrix(c(-1, 0, 1, -1),
nrow = 2,
byrow = TRUE),
matrix(c(1, 0, 1, -1),
nrow = 2,
byrow = TRUE)),
node = list(2, 2, 1),
solver = "ECOS",
B = 10))

@Maggie8888
Copy link
Collaborator

Since your sample size is too small most of the statistical test won't be stable. If you can somehow increase your sample size that will be a better way to solve the issue.

@Maggie8888 Maggie8888 added the invalid This doesn't seem right label Sep 27, 2024
@UniAlberta
Copy link
Author

Thanks for your reply. There`s no way to increase the sample size because they were taken from agricultural soils. So if I repeat ancombc2 with a low sample size multiple times and get the same result, can I trust them?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
invalid This doesn't seem right
Projects
None yet
Development

No branches or pull requests

2 participants