Question

我需要在R中进行一些咨询。我试图找到一种方法来保存＆＃34; 租赁＆＃34;在这个模拟中，在一个单独的矩阵中。

Part_of_my_simulation_is_described_below

for(i in 1:n){

for(j in 1:antall.simuleringer){

**rental**<-rbinom(1,1,ikke.reforhan)

L1KS[j,i]<-if(L1kontrakt>=i){L1kontr.tot} else if.

有人有什么建议吗？

Answer 1

我会像下面那样进行总结，在那里我将简单可重复的求和结果放在一个列表中。

# seed for reproducibility
set.seed(1)

# define a list
my_list <- list()

# my sequence
my_seq <- 1:10

# sumulation block of code
for(i in my_seq){

  # matrices of random numbers
  # store matrices in a list
  my_list[[i]] <- matrix(round(rnorm(20), 2), nrow = 2)
  my_list
}

# a look at results
my_list[1:2]

[[1]]
      [,1]  [,2]  [,3] [,4]  [,5] [,6]  [,7]  [,8]  [,9] [,10]
[1,] -0.63 -0.84  0.33 0.49  0.58 1.51 -0.62  1.12 -0.02  0.82
[2,]  0.18  1.60 -0.82 0.74 -0.31 0.39 -2.21 -0.04  0.94  0.59

[[2]]
     [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10]
[1,] 0.92  0.07  0.62 -0.16 -0.48  1.36  0.39 -1.38 -0.39  1.10
[2,] 0.78 -1.99 -0.06 -1.47  0.42 -0.10 -0.05 -0.41 -0.06  0.76

Answer 2

我不完全理解你的例子，但至少，伯努利随机变量不需要递归创建。

n <- 5
m <- 5
prob <- 0.5
df <- matrix(rbinom(n * m, 1, prob), nrow = n)
df
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    1    0    1    1
[2,]    1    1    0    0    0
[3,]    0    0    0    0    1
[4,]    0    0    1    0    0
[5,]    1    1    0    0    1

我不能再继续了，因为我不确定如何调理。

在模拟中创建矩阵（在R中）

2 个答案: