anjutalq 发表于 2013-10-16 21:18 
同学,你好。首先谢谢你的热心回复,对我帮助很大。我觉得我还有几个疑惑想请教你,我是工科的学生,学习 ...
同学,不好意思,我好像误导你了,回复完你之后我越想越觉得奇怪,可能不是这样的。
分两步解释:
关于mle,可以参见一面这篇文章:
http://wenku.baidu.com/view/a9fac2946bec0975f465e2bf.html
文章里第一页讲了MLE的两个性质:无偏性和渐进正态性,尤其是渐进正态性是解决你的问题的关键。
渐进正态性是指,当样本量足够大的时候,极大似然估计值服从一个正态分布。由于你的样本有10000个,非常大的一个样本,所以可以很有信心地说,你的估计值的分布一定是正态的。那么由此就可以直接说出估计量的置信区间了。因此,验证这个区间实际上是一个伪命题,因为已经可以由理论指导估计量的置信区间了。
换句话说,matlab给你的置信区间在中心极限极限定理下是一定对的,不需要检验。你老师的要求怎么满足我也不知道了。。。
另一部分是,我的解释和检验方法是错的,我有点想当然了。希望没有造成太大的困扰。
关于你的疑问的回答:
1. 你的理解不太对,也是我犯错的地方。90%置信区间是指,这个分布的真实的参数落在这个区间的概率是90%。更具体的说,如果你从一个标准分布产生一组随机数(N(0,1)),并根据这组随机数估计均值的置信区间,那么这个区间包含0的概率是90%(0是真值)。。理解了吗?
matlab实验代码:
k=0;
for i=1:10000
x=randn(1000,1); %standard normal random numbers
mu_est=mean(x); sigma_est=std(x);
lower=mu_est-1.96*sigma_est; upper=mu_est+1.96*sigma_est;
%lower是95%置信区间下限,upper是上限,理论上这个区间包含0的概率是95%,不过你再check一下1.96这个数对不对
if lower<0 & upper>0
k=k+1;
end
end
计算完成后,k的值大约为9500。希望能帮助你理解。
2. 参数置信区间是不能检验的,这是更正。
剩下的应该没什么问题了。欢迎讨论。