模拟钙基线的真实噪音

Question

由于真正令人惊叹的社区，我的项目组更接近于模仿真实的钙基线噪音。

我在数学模型中模拟了典型的钙运动：

感谢社区，我可以将随机噪音添加到不切实际的基线：

然而，噪音动态实际上太快了。有没有办法减慢噪音，创造更广泛的噪音峰值，而不是这些峰值。我添加了一个实际测量值来向您显示我的意思：

如果这个问题太具体，我道歉并删除帖子。

祝愿和非常感谢！

Answer 1

请提出问题和示例reproducible，以便其他人可以提供帮助。话虽如此，看起来基线只是一个随机的正常 - 可能是用x <- rnorm(500)之类的东西创建的。使这种不那么跳跃的一种方法是计算移动平均线。您可以使用TTR或zoo之类的包来执行此操作，也可以创建自己的函数。例如：

x <- rnorm(500)
plot(x, type = "l")

ma <- function(x, n = 5){ filter(x, rep(1/n, n), sides = 2) }
plot(ma(x), type = "l")

plot(ma(x, 10), type = "l")

Answer 2

我现在明白你的意思了。我对这个案子有两点建议，也许他们会有所帮助：

尝试仅将噪音添加到基线的一个子集（以下是10％）

baseline.index = which(App[,2] == min(App[,2]))
baseline.index.subset = sample(x = baseline.index, size = 0.1 * length ( baseline.index) , replace = F)
noise = rnorm( length (baseline.index.subset))
App[ baseline.index.subset,2] = App[ baseline.index.subset,2] + noise

尝试用噪音的平均值和标准差来玩。即：

noise = rnorm( length (baseline.index.subset), mean = 0, sd = 0.1)

如果有帮助，请告诉我们