Question

我需要重复以下这样的过程10次：

> random=rnorm(5)  
> total=sum(random)

然后将每次总数放入一个名为final的向量中。

对于循环，我想我会写：

> final <- c(total[[i]])
> for(i in 1:10){  
  random=rnorm(5)  
  total=sum(random)
  ...
}

...不确定如何将我的总数放入我的最终载体中？另外，从在线浏览大多数地方说我应该使用申请代码。我应该使用哪个，如何将总数放入最终？

Answer 1

R为这个用例提供了一个名为replicate的函数，它将按指定的次数重复操作，并将结果存储在向量中：

> replicate(10, sum(rnorm(5)))
[1] -0.2286870  1.5123902 -5.8179539  1.2119908 -2.8749987  5.5817021
[7]  2.5427969  2.1833426  0.4884455 -3.6787912

如果生成数据的函数更复杂，那么您可以使用replicate来调用函数，也可以直接使用sapply。所以

> genRandomData <- function() {
    random <- rnorm(100)
    c(sum(random), var(random))
}
> replicate(10, genRandomData())
         [,1]     [,2]      [,3]      [,4]      [,5]      [,6]       [,7]
[1,] 6.612235 6.895032 -4.004501 -8.250871 2.3098680 2.0388161 23.5202949
[2,] 1.082576 0.769107  1.119480  1.206932 0.9930102 0.9078762  0.9544335
            [,8]      [,9]     [,10]
[1,] -12.2080648 9.0074524 0.2358715
[2,]   0.8256081 0.8669045 0.9017415

或者，使用sapply产生相同的输出（这很好，因为它不会在替换时回复），

> sapply(integer(10), function(x) {
    random <- rnorm(100)
    c(sum(random),var(random))
})

如果你想获得一个数据框，你必须转置;所以

> data.frame(t(sapply(integer(10), function(x) {
    random <- rnorm(100)
    list(s=sum(random),v=var(random))
})))
            s         v
1   -6.343845    1.0397
2    12.94773 0.9480807
3    5.811322 0.9670198
4   -1.941444   1.04547
5  -0.9094589  0.862759
6   -2.852641 0.5504582
7   -8.471266  1.080554
8    -17.2341  1.201679
9    11.60805  1.150254
10  -7.138314  1.080731

使用apply或loop输出具有预定长度的向量

1 个答案: