This function is useful for plotting time-courses for individual features.
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
metadatamust correspond to columns ofy. Row names are ignored.- features
Vector of names, row numbers, or logical values for subsetting the features.
NULLindicates all features.- sampleColname
String indicating the column in
metadatacontaining the name of each sample, which must correspond to the column names ofy.
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'))