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