Use this function to find for each possible subset of drug combinations out of a given list of synergies, the number of models that predicted it given the models' predictions. So, if for example the set of synergies is this one: {'A-B','C-D','E-F'}, we want to know how many models predicted none of them, just the single subsets (e.g. the {'A-B'}), the two-element subsets (e.g. the {'A-B','C-D'}) and all 3 of them.
get_synergy_subset_stats(model.predictions, synergies)
a data.frame
object with rows the models and
columns the drug combinations. Possible values for each model-drug combination
element are either 0 (no synergy predicted), 1 (synergy was
predicted) or NA (couldn't find stable states in either the drug
combination inhibited model or in any of the two single-drug inhibited models).
a character vector with elements the synergistic drug
combinations. Note that these synergies should be a subset of the column
names of the model.predictions
data.frame.
an integer vector with elements the number of models the predicted each synergy subset. The names attribute has the names of each synergistic drug combination subset, which are the drug combinations comma separated (e.g. 'A-B,C-D').
Note that if the synergies
vector has more than 10-15 elements, then
this function might take long time to execute even with an optimal
implementation of count_models_that_predict_synergies
.