添加極值分佈的噪聲（帶μ= 0，σ= 10）至數的向量中的R

Question

我有以下矩陣

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中實現這一點？

我在extRemes包中發現了revd包，但它不能按預期工作。來自同一包裝的devd做我想做的事情嗎？ （但它不允許均值和標準差來定義）

Answer 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 

編輯：您可避免for循環代替：

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