24小时热门版块排行榜

返回列表

baobao007

铁杆木虫 (正式写手)

应助: 2 (幼儿园)
金币: 9720.8
散金: 100
帖子: 344
在线: 578.5小时
虫号: 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语言代码！！！！）

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有5人回复
面上没中，邀请各位路过的虫友分析一下分数已经有7人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
麻烦专家们看看评委们的意见（F口面上）已经有13人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复
小白求助投论文要求的highlights应该如何写已经有4人回复

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的回帖