Spaghetti plot

Description

Spaghetti plot for longitudinal data

Usage

spaghetti(formula, data, id = "id", group = NULL, type = "l", lty = 1,
  col = 1:10, alpha = 0.3, lwd = 1, trend.formula = formula,
  tau = NULL, trend.lty = 1, trend.join = TRUE, trend.delta = 0.2,
  trend = !is.null(tau), trend.col = col, trend.alpha = 0.2,
  trend.lwd = 3, legend = NULL, xlab = "Time", ylab = "", add = FALSE,
  ...)

Arguments

Author(s)

Klaus K. Holst

Examples

if (interactive() & requireNamespace("mets")) {
K <- 5
y <- "y"%++%seq(K)
m <- lvm()
regression(m,y=y,x=~u) <- 1
regression(m,y=y,x=~s) <- seq(K)-1
regression(m,y=y,x=~x) <- "b"
d <- sim(m,500)
dd <- mets::fast.reshape(d);
dd$num <- dd$num+rnorm(nrow(dd),sd=0.5) ## Unbalance
spaghetti(y~num,dd,id="id",lty=1,col=Col(1,.4),trend=TRUE,trend.col="darkblue")
}

Results