drbats.simul: Main simulation function

In DrBats: Data Representation: Bayesian Approach That's Sparse

Main simulation function

Description

Main simulation function

Usage

drbats.simul(
  N = 10,
  P = 150,
  t.range = c(0, 1000),
  b.range = c(0.2, 0.4),
  c.range = c(0.6, 0.8),
  b.sd = 2,
  c.sd = 2,
  a.range = c(-0.4, 0.4),
  y.range = c(0, 10),
  amp = 10,
  per = 12,
  data.type = "sparse",
  breaks = 15,
  sigma2 = 0.2,
  seed = NULL
)

Arguments

N	integer number of functions to simulate (default = 10)
P	a number of observation times (default = 150)
t.range	a range of times in which to place the P observations (default = c(1, 1000))
b.range	a vector giving the range of values for the mean of the first mode (default b.range = c(0.2, 0.4))
c.range	a vector giving the range of values for the mean of the second mode (default c.range = c(0.6, 0.8))
b.sd	the standard deviation for the first mode (default b.sd = 2)
c.sd	the standard deviation for the second mode (default c.sd = 2)
a.range	a vector giving the range of values for the slope (default a.range = c(-0.4, 0.4))
y.range	a vector giving the range of values for the intercept (default y.range = c(0, 10))
amp	the amplitude of the cosine function (default = 10)
per	the periodicity of the cosine function (default = 12)
data.type	string indicating type of functions (options :sparse, sparse.tend, sparse.tend.cos)
breaks	number of breaks in the histogram basis
sigma2	the precision of the error terms (default = 0.2)
seed	integer specification of a seed (default = NULL)

Value

Y.simul a list containing a matrix Y, a matrix beta, and a matrix epsilon

t.simul a matrix of simulated observation times

X the underlying signal to build the data, see DataSimulationandProjection vignette

proj.pca the outputs of the function pca.proj.Xt

wlu the outputs of the function W.QR

Author(s)

Gabrielle Weinrott

Examples

res <- drbats.simul(N = 5, P = 100, t.range = c(5, 100), breaks = 8)
X <- res$X
t <- res$t.simul
# To plot the observations, ie the rows
matplot(t(t), t(X), type = 'l', xlab = "Time", ylab = "X")

DrBats documentation built on March 18, 2022, 5:15 p.m.