R语句 蒙特卡罗模拟

chm0508

求教大神！！
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
哪里有问题呢？报错如下

chm0508

已解决！循环语句用{}