Dodanie ekstremalnej wartości szumu rozproszonego (z μ = 0, σ = 10) do wektora liczb w R - r, losowy, szum

Mam następującą matrycę

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

Chcę dodać szum rozproszony o ekstremalnych wartościach (ze średnią = 0 i sd = 10) do wartości pomiaru. Jak mogę to osiągnąć w R?

znalazłem revd w extRemes pakiet, ale nie działa zgodnie z oczekiwaniami. Robi devd z tego samego pakietu zrobić to, co chcę zrobić? (ale nie pozwala na zdefiniowanie średniej i SD)

Odpowiedzi:

Jeśli chcesz użyć swojej miary jako średniej hałasu, możesz to zrobić:

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

EDYCJA: Zamiast tego możesz uniknąć pętli for:

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