Skip to contents

Computes summaries of the data matrix grouped by combinations of design variables.

Usage

# S3 method for class 'multidesign'
summarize_by(x, ..., sfun = colMeans, extract_data = FALSE)

Arguments

x

A multidesign object

...

Unquoted names of variables to group by

sfun

Summary function to apply (default is colMeans)

extract_data

Logical; whether to extract raw data instead of computing summary

Value

A new multidesign object containing:

x

Matrix of summary statistics

design

Design information for each summary

column_design

Original column metadata

Examples

X <- matrix(rnorm(100*20), 100, 20)
Y <- tibble::tibble(
  condition = rep(c("A", "B"), each=50),
  block = rep(1:2, times=50)
)
mds <- multidesign(X, Y)

# Get means by condition
means_by_cond <- summarize_by(mds, condition)

# Get means by condition and block
means_by_both <- summarize_by(mds, condition, block)