simulate_linear: Simulate linear trends

In ssokolen/metcourse: Tools for Metabolic Time-course Data

Description

Simulates linear trends with p fraction decreasing and (1-p) fraction increasing.

Usage

simulate_linear(n, n.samples, p)

Arguments

n	Number of trends to generate.
n.samples	Number of timepoints within each trend.
p	Proportion of trends decreasing (vs. increasing).

Value

A dataframe with each row representing a metabolite concentration time-course scaled between 0 and 1. The corresponding x variables are assumed to be equally spaced between 0 and 1 i.e. x <- seq(0, 1, length.out=n.samples).

Examples

trends <- simulate_linear(1000, 100, 0.75)

# Conversion for plotting
y_mat <- t(as.matrix(trends))
x <- seq(0, 1, length.out=100)
matplot(x, y_mat, type = 'l', lty = 1, lwd = 4, col = grey(0, 0.05),
        xlab = 'Relative culturing time', ylab = 'Relative concentration')

ssokolen/metcourse documentation built on May 30, 2019, 8:43 a.m.