Simulation de bruit réaliste pour une ligne de base de calcium

Question

Grâce à la communauté vraiment incroyable, mon groupe de projet est un pas de plus pour imiter le bruit de base réaliste de calcium.

J'ai simulé un mouvement de calcium typique dans un modèle mathématique:

Merci à la communauté que je pourrais ajouter un bruit aléatoire à la ligne de base irréaliste:

Cependant, la dynamique du bruit est en fait trop vite. Y at-il un moyen de ralentir le bruit et de créer des pics de bruit plus larges au lieu de ces pics. J'ajoute une mesure réelle pour vous montrer ce que je veux dire:

Si cette question est trop spécifique, je présente mes excuses et supprimerai le poste.

Meilleurs voeux et merci!

Answer 1

Veuillez faire votre question et les exemples reproducible afin que d'autres puissent vous aider. Cela étant dit, il semble que la ligne de base est juste une normale aléatoire - probablement créée avec quelque chose comme x <- rnorm(500). Une façon de rendre cela moins nerveux est de calculer une moyenne mobile. Vous pouvez utiliser un package comme TTR ou zoo pour ce faire, ou vous pouvez créer votre propre fonction. Par exemple:

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

Je vois votre point maintenant. J'ai deux suggestions pour ce cas, peut-être qu'ils seront utiles:

Essayez d'ajouter du bruit à un sous-ensemble de votre ligne de base (ci-dessous est un 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 

Et essayer de jouer un bit avec la moyenne et l'écart type du bruit. à savoir:

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

Faites nous savoir si cela aide