Perform differential protein expression analysis

Usage

get_DEPresults(
  se,
  condition1 = NULL,
  condition2 = NULL,
  ref_condition = NULL,
  tests = NULL,
  alpha = 0.05,
  lfc = 1,
  type = "manual",
  fdr_type = "BH"
)

Arguments

se: SummarizedExperiment object with protein expression data.
condition1: One of the conditions to be tested present in expDesign.
condition2: One of the conditions to be tested present in expDesign.
ref_condition: One of the conditions in expDesign that is used as reference. Must be used in combination with type = 'control'.
tests: Character vector specifying which contrasts to test when using "type = 'manual'"
alpha: the p.adjust significance cutoff.
lfc: The log2 fold change significance cutoff.
type: Type of comparison to be made. Options are: 'manual', 'control', and 'all'.
fdr_type: Type of fdr correction. Options are 'fdrtool' and 'BH' (Benjamini-Hochberg). Default is 'BH'

Value

A data frame with the results of differential protein expression analysis

Examples

se <- prepare_se(report.pg_matrix, expDesign)

#> Imputing along margin 1 (features/rows).
#> Warning: 36 rows with more than 50 % entries missing;
#>  mean imputation used for these rows
#> Cluster size 5548 broken into 2203 3345 
#> Cluster size 2203 broken into 1569 634 
#> Cluster size 1569 broken into 1504 65 
#> Cluster size 1504 broken into 645 859 
#> Done cluster 645 
#> Done cluster 859 
#> Done cluster 1504 
#> Done cluster 65 
#> Done cluster 1569 
#> Done cluster 634 
#> Done cluster 2203 
#> Cluster size 3345 broken into 2136 1209 
#> Cluster size 2136 broken into 1087 1049 
#> Done cluster 1087 
#> Done cluster 1049 
#> Done cluster 2136 
#> Done cluster 1209 
#> Done cluster 3345 
res <- get_DEPresults(se, 'motif1', 'neg_ctrl', type = 'manual')
#> Tested contrasts: motif1_vs_neg_ctrl