24小时热门版块排行榜

>论坛更新日志 (4418)
>基金申请 (960)
>文献求助 (320)
>论文投稿 (212)
>虫友互识 (171)
>论文道贺祈福 (136)
>硕博家园 (122)
>考博 (95)
>休闲灌水 (92)
>导师招生 (75)
>博后之家 (55)
>考研 (48)
>公派出国 (45)
>教师之家 (44)
>绿色求助(高悬赏) (43)
>找工作 (38)

返回列表

baobao007

铁杆木虫 (正式写手)

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

回复此楼

» 猜你喜欢

生命口C13面上会评时间已经有11人回复
内心匮乏已经有16人回复
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？已经有9人回复
有机物中金属离子降到1ppm以下，有什么好的方法吗？用过有机溶剂重结晶还是降不下来已经有4人回复
投英文期刊如何查重已经有6人回复
不要再数国自然申请书的 filecode 的分隔符个数了已经有21人回复
伯胺上甲基求助已经有7人回复
跨出版社商投稿已经有5人回复
5个%2F，非固定段22个字母已经有7人回复
生生命学部会评已经有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的回帖