Имам следната матрица
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
Искам да добавя екстремна стойност разпределен шум (със средно = 0 и sd = 10) към стойностите на измерването. Как мога да постигна това в R?
намерих revd
в extRemes
пакет, но не работи както се очаква. Прави devd
от същия пакет правя това, което искам да правя? (но не позволява определянето на средно и sd)
Отговори:
0 за отговор № 1Ако искате да използвате мярката си като средно за шума, можете да направите това:
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: Можете да избегнете цикъла за с това вместо това:
data$measure1 = data$measure + rnorm(1,0,10)