Tengo la siguiente matriz
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
Quiero agregar ruido distribuido de valor extremo (con media = 0 y sd = 10) a los valores de medición. ¿Cómo puedo lograr eso en R?
encontré revd
en extRemes
paquete, pero no funciona como se esperaba. Lo hace devd
¿Del mismo paquete hago lo que quiero hacer? (pero no permite definir media y sd)
Respuestas
0 para la respuesta № 1Si desea utilizar su medida como la media del ruido, puede hacer esto:
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
EDITAR: Puede evitar el bucle for con esto en su lugar:
data$measure1 = data$measure + rnorm(1,0,10)