0

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

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

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

对于一个循环，我想我会写：

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

...不知道如何将我的总数放入我的最终向量中？同样从网上环顾四周，大多数地方都说我应该使用应用代码。我应该使用哪个以及如何将总数放入最终结果？

4

1 回答 1

2

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

如果产生数据的函数比较复杂，那么你可以使用replication来调用一个函数，或者直接使用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
于 2013-01-14T14:35:30.327 回答