4楼朋友，您好！
您编的程序为什么在R2.12.0中运行不了？有以下错误提示，您能帮我看看吗？
如何解决？
非常非常感谢！
已加载了fBasics包

> #画出r=1，2，…，N和相应m个序列平均方差的关系图
> var.mean=rep(0,N)
错误: 找不到对象'N'
> for(r in 1:N)
+ { m=choose(N,r)
+ SA=array(0,dim=(m,T,r));ind=sample(1:N,r)
错误: 意外的',' 于
"{ m=choose(N,r)
SA=array(0,dim=(m,"
> SA.mean=rep(0,m)
错误: 找不到对象'm'
> for(i in 1:m)
+ {  SA[i,,]=Q[,ind]
+    SA.mean=var(apply(SA[i,,],2,mean))
+ }
错误: 找不到对象'm'
> var.mean[r]=mean(SA.mean)
错误于mean(SA.mean) : 在为'mean'函数选择方法时评估'x'参数出了错
> }
错误: 意外的'}'在"}"里
> plot(1:N,var.mean,type=”o”)
错误: 意外的输入在"plot(1:N,var.mean,type=?"}"里
>

6# elephann

第一：如果你的数据量非常大，计算量将非常大，你的题目没有讲清楚，如果这N个样本是来自统一总体，想要验证其正态性，直接把N个样本数据合并成一组数据来整整其正态性！！
第二：对于你的第二个问题   “要求对一个总体中N任意取出m个序列的等权新组合、再减去另一对照序列差额的方差.............”  下面这段程序就是根据你的意思来编的，在逻辑上是行得通的，如果有偏差只能说明你的数据或者所要求的方法有问题，另外我没有你的数据，我也无法验证你的结果。
画出r=1，2，…，N和相应m个序列平均方差的关系图
var.mean=rep(0,N)
for(r in 1:N)
{ m=choose(N,r)
SA=array(0,dim=(m,T,r));ind=sample(1:N,r)
SA.mean=rep(0,m)
for(i in 1:m)
{
SA[i,,]=Q[,ind]

SA.mean=var(apply(SA[i,,],2,mean))
}
var.mean[r]=mean(SA.mean)
}


Ps：对于别人热情的参与拟的帖子要给予评分