我对R很新,我试图对数据框的每个元素进行操作,如下所示:
> head(data)
V1 V2 V3 V4 V5 V6 V7 V8 V9
1 3.55271e-15 10.7349 21.9238 26.5136 42.38290 37.0355 28.6789 36.6033 26.46900
2 0.00000e+00 14.5497 19.4960 23.9728 37.66170 32.5874 32.9722 46.7856 41.23690
3 0.00000e+00 19.8794 12.6273 15.7562 26.07470 16.3857 33.3016 40.0800 33.45120
4 0.00000e+00 13.9638 33.8372 50.1383 52.74920 55.9552 79.8752 78.2309 76.51520
5 0.00000e+00 25.1614 23.7001 21.1502 7.11349 13.3939 18.4262 1.1183 8.78106
6 3.55271e-15 11.1300 10.5257 21.4103 8.95956 14.8252 42.4594 40.8457 29.57740
我想要做的是将函数应用于数据框中的每个元素,然后平均每个k的结果:
for (k in c(1,100){
sample <- cos(k*data)
average <- mean(sample)}
但是我意识到这将覆盖样本变量,我只留下一个样本(k = 100)。
我也知道R更喜欢矢量化函数而不是循环,所以我该如何处理这个问题呢?谢谢!
答案 0 :(得分:1)
您可以使用sapply()将函数应用于数据并返回平均值。
sapply(1:100, (function(i){
sample <- cos(i*data)
mean(sample)
}))
或者,如果您想存储&#34;样品&#34;在列表中,您可以先在lapply()上使用,然后再在样本列表中使用sapply()。
samples <- lapply(1:100, (function(i){
cos(i*data)
}))
averages <- sapply(samples, mean)
head(averages)
[1] 0.046893344 0.032709183 -0.009479597 -0.185085895 -0.395590101 0.554012071