Skip to contents

This function is useful for plotting time-courses for individual features.

Usage

mergeMeasMeta(y, metadata, features = NULL, sampleColname = "sample")

Arguments

y

Matrix-like object of measurements, with rows corresponding to features and columns to samples.

metadata

data.frame containing experimental design information for each sample. Rows of metadata must correspond to columns of y. Row names are ignored.

features

Vector of names, row numbers, or logical values for subsetting the features. NULL indicates all features.

sampleColname

String indicating the column in metadata containing the name of each sample, which must correspond to the column names of y.

Value

A data.table with one row for each sample-feature pair.

Examples

library('data.table')

y = GSE34018$y
metadata = GSE34018$metadata

fit = getModelFit(y, metadata)
fit = getPosteriorFit(fit)

measObs = mergeMeasMeta(y, metadata, features = c('13170', '12686'))
measFitMean = getExpectedMeas(
  fit, times = seq(0, 24, 0.5), features = c('13170', '12686'))