getFDR: Get FDR for given F statistics based on true and null dataset

Description Usage Arguments Value Examples

View source: R/fdr_funcs.R

Description

Get FDR for given F statistics based on true and null dataset

Usage

1
getFDR(df_out, df_null, squeezeDenominator = TRUE)

Arguments

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

Value

data frame annotating each protein with a FDR based on it's F statistic and number of observations

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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)
 

TPP2D documentation built on Nov. 8, 2020, 4:54 p.m.