Aggiunta di rumore distribuito di estremo valore (con μ = 0, σ = 10) a un vettore di numeri in R - r, casuale, rumore

Ho la seguente matrice

Measurement Treatment
38 A
14 A
54 A
69 A
20 B
36 B
35 B
10 B
11 C
98 C
88 C
14 C

Voglio aggiungere un valore estremo di rumore distribuito (con media = 0 e sd = 10) ai valori di misurazione. Come posso raggiungerlo in R?

ho trovato revd in extRemes pacchetto, ma non funziona come previsto. fa devd dallo stesso pacchetto fai quello che voglio fare? (ma non consente di definire mean e sd)

risposte:

Se vuoi usare la tua misura come media per il rumore, puoi farlo:

measure = round(runif(10,0,30),0)
data = data.frame(measure)
for(i in 1:nrow(data)){
data$measure1[i] = rnorm(1,data$measure[i],10)
}
data
measure  measure1
1        6  6.281557
2       12 -5.780177
3       18 13.529773
4       26 33.665584
5       14 12.666614
6       24 41.146132
7        5 -1.850390
8       14 16.728703
9       13 26.082601
10      13 14.066475

EDIT: puoi evitare il ciclo for con questo invece:

data$measure1 = data$measure + rnorm(1,0,10)