我有两个数据集的形式$q \times k$矩阵（变量$x$和$y$）。两个变量都不同（假设它们相差一个常数）。矩阵的列是后续时间点，行是观察值（样本分位数）。问题是可视化这两个变量以显示它们之间的差异及其随时间的变化。原始数据可用。情节需要外行人可读。

下面我贴上可以作为说明的合成数据。

set.seed(123)
mu <- 0
n <- 1000
k <- 10
q <- 10
x <- matrix(NA, q, k)
y <- matrix(NA, q, k)

for (i in 1:k) {
  x[, i] <- quantile(rnorm(n, mu[i], 1), c(seq(0, 1, by=.1)-.05)[-1])
  y[, i] <- quantile(rnorm(n, mu[i]+1.5, 1), c(seq(0, 1, by=.1)-.05)[-1])
  mu[i+1] <- mu[i] + runif(1, -1, 1)
}

你有什么想法如何可视化？