找回密码
 立即注册
查看: 6106|回复: 1

R语句 蒙特卡罗模拟

[复制链接]
发表于 2015-3-10 19:49:50 | 显示全部楼层 |阅读模式
求教大神!!
X服从(0,b)的均匀分布,n=100
两种方案对b进行估计:
       方案A:  b=2倍的X的平均值
       方案B:  b=max{X1,X2,...,Xn}
假设b=10,用蒙特卡洛模拟衡量两个方案的好坏,衡量标准为MSE=Var(X)+(E(X)-b)^2
程序如下:
>A<-rep(0,10^5)
>B<-rep(0,10^5)
>for(i in 1:10^5)
+(X.sam<-runif(100,0,10)
+A<-2*mean(X.sam)
+B<-max(X.sam))
>var(A)+(mean(A)-10)^2
[1]0.335057
>var(B)+(mean(B)-10)^2
[1]0.01949572
哪里有问题呢?报错如下
QQ截图20150310193133.png

回复

使用道具 举报

 楼主| 发表于 2015-3-10 20:46:52 | 显示全部楼层
已解决!循环语句用{}
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|R语言中文网

GMT+8, 2024-11-25 20:53 , Processed in 0.030441 second(s), 20 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表