Performs an ANOVA statistical test

anova_protti(data, grouping, condition, mean_ratio, sd, n)

Arguments

data

a data frame containing at least the input variables.

grouping

a character column in the data data frame that contains precursor or peptide identifiers.

condition

a character or numeric column in the data data frame that contains the conditions.

mean_ratio

a numeric column in the data data frame that contains mean intensities or mean intensity ratios.

sd

a numeric column in the data data frame that contains the standard deviation corresponding to the mean.

n

a numeric column in the data data frame that contains the number of replicates for which the corresponding mean was calculated.

Value

a data frame that contains the within group error (ms_group) and the between group error (ms_error), f statistic and p-values.

Examples

data <- data.frame(
  precursor = c("A", "A", "A", "B", "B", "B"),
  condition = c("C1", "C2", "C3", "C1", "C2", "C3"),
  mean = c(10, 12, 20, 11, 12, 8),
  sd = c(2, 1, 1.5, 1, 2, 4),
  n = c(4, 4, 4, 4, 4, 4)
)

anova_protti(
  data,
  grouping = precursor,
  condition = condition,
  mean = mean,
  sd = sd,
  n = n
)
#> # A tibble: 2 × 5
#>   precursor ms_group ms_error     f      pval
#>   <chr>        <dbl>    <dbl> <dbl>     <dbl>
#> 1 A            112       2.42 46.3  0.0000183
#> 2 B             17.3     7     2.48 0.139