getFDR {TPP2D} | R Documentation |
Get FDR for given F statistics based on true and null dataset
getFDR(df_out, df_null, squeezeDenominator = TRUE)
df_out |
data frame containing results from analysis by fitAndEvalDataset |
df_null |
data frame containing results from analysis by bootstrapNull |
squeezeDenominator |
logical indicating whether F statistic denominator should be shrinked using limma::squeezeVar |
data frame annotating each protein with a FDR based on it's F statistic and number of observations
data("simulated_cell_extract_df") temp_df <- simulated_cell_extract_df %>% filter(clustername %in% paste0("protein", 1:5)) %>% group_by(representative) %>% mutate(nObs = n()) %>% ungroup example_out <- fitAndEvalDataset(temp_df) example_null <- bootstrapNull(temp_df, B = 1) getFDR(example_out, example_null)