24小时热门版块排行榜

返回列表

baobao007

铁杆木虫 (正式写手)

应助: 2 (幼儿园)
金币: 9676.8
散金: 100
帖子: 344
在线: 575.9小时
虫号: 2211769
注册: 2012-12-28
专业: 美学

[求助] 蒙特卡洛模型的-R语言代码！！！！！已有1人参与

最近想用蒙特卡洛模型来验证健康风险模型的有效性（如下）
v=rnorm(10000,2,0.2)
b=rnorm(10000,70,7)
c=rlnorm(10000,4.6,2.3)
f=rnorm(10000,0.03,0.003)
for(i in 1:10000){
m=c*v*f*0.001/b}
（假设m就是计算健康风险的公式，下面想输出置信区间在5%到95%之间的计算结果，怎么写R语言代码！！！！）

回复此楼

» 猜你喜欢

1楼 2016-03-26 21:56:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hukaimmiww

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 451.6
帖子: 112
在线: 29.8小时
虫号: 181352
注册: 2006-02-10
专业: 环境生物物理

N <- 9999
m_sim <- rep(NA,N)
for(i in 1:N){
  v=rnorm(100,2,0.2)
  b=rnorm(100,70,7)
  c=rlnorm(100,4.6,2.3)
  f=rnorm(100,0.03,0.003)
  m_sim =c*v*f*0.001/b
}
head(m_sim,10)

这样？

赞一下

回复此楼

2楼2019-08-22 20:39:58

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hukaimmiww

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 451.6
帖子: 112
在线: 29.8小时
虫号: 181352
注册: 2006-02-10
专业: 环境生物物理

【答案】应助回帖

N <- 9999
m_sim <- rep(NA,N)
for(i in 1:N){
  v=rnorm(100,2,0.2)
  b=rnorm(100,70,7)
  c=rlnorm(100,4.6,2.3)
  f=rnorm(100,0.03,0.003)
  m_sim =c*v*f*0.001/b
}
head(m_sim,10)

赞一下

回复此楼

3楼2019-08-22 20:48:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖