很难说出这个问题的标题,但有什么,是一个xts向量:
> test
E1 E2 E3 E4 FP5
2012-05-17 03:34:37 4045.75 4045.75 2835.5 1292 171.9
现在,我想将正态分布的结果添加到每个测试值,但我希望sd参数(标准偏差)是值本身的sqrt
对于isntance,我将采用非自动化方法:
test$E1 = test$E1 + rnorm(1, sd = sqrt(test$E1))
test$E2 = test$E2 + rnorm(1, sd = sqrt(test$E2))
...
以任何方式更简单地执行此操作?
答案 0 :(得分:2)
你可以通过在test:
列上调用apply来执行此操作test = apply(test, 2, function(x) x+rnorm(length(x), sd=sqrt(x)))
更长的形式(没有单行函数语法):
adjust.values = function(x) {
return(x + rnorm(length(x), sd=sqrt(x)))
}
test = apply(test, 2, adjust.values)