蒙特卡洛中的种子

Question

#year 2022

R.100<-runif(n=100, min = 40, max = 65)
R.100
summary(R.100)

R.norm.100<-rnorm(length(R.100), mean = 50, sd=7)
R.norm.100
summary(R.norm.100)

a=0.6
C=100

f<-seq(from=1.10, to=1.18, by=0.005)
B<-NULL
B<-as.data.frame(B)

for(i in 1:length(f)) {
  for(R in 1:length(R.norm.100)) {
    B[i, "degerler"] <- (R.norm.100[R]-(1-a)*C*f[i]) / R.norm.100[R]
  }
}

summary(B)

[运行上述循环函数并计算B时，我希望分子中R的模拟数量在分母中相同。例如，如果在分子中将R模拟为55时计算B，我希望分母中的R是55。

如何为该功能编写R代码？

最后一件事，我想多次重复执行此功能，例如50。

您能帮我吗？最好的祝福Zehra Civan

＃year 2022 R.100

Answer 1

我不确定我了解您要做什么。 R.norm.100是一个向量，在分子和分母中将是相同的值。您的代码使用两个循环，R中很少使用这两个循环，因为您正在使用的功能已向量化。在您的情况下，B[i, "degerler"]将被覆盖100次，因此仅保留最后一个值。