Hi there!
I'm trying to using the MultiGroupPower function to run an 'omic power analysis between microbiome data (MGS) and metabolite data for a group variable with three levels. I've tried to execute the function as follows:
MultiGroupPower(data_list, groups, type = c(1,2), comparisons = NULL, omicPower = 0.6,
averagePower = 0.8, fdr = 0.05, cost = 1, equalSize = TRUE, max.size = 10000,
omicCol = NULL, powerPlots = TRUE, summaryPlot = TRUE)
but keep getting the error:
Error in if (pow.mid > ave.pow) { : missing value where TRUE/FALSE needed
I'm not sure what's gone wrong, as I've input the data types correctly, filtered out all NA values in my grouping vectors, and was able to get the same data to work using a two-group comparison by Sex with the vanilla MultiPower function.
Any advice would be appreciated when interpreting this error!