放弃！R里面做不出来的。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

放弃！数据太大，R软件无法解决！！！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：r软件 论坛 求助 悬赏

加大悬赏额!!!

1000块大洋啊~~

1# elephann

仅供参考~~~

4# 楚韵荆风
谢谢你花费时间编程，待我试验验证若能有助于解决问题一定兑现承诺！先仅表示感谢！此外，其他朋友若有好的想法也不妨继续帮忙，若有助于问题解决或者提供好的思路都未尝不可继续奖励！！

4# 楚韵荆风
1 关于描述性统计有点疑问，每个时间序列的统计量这容易，但如果我有1000或者2000，我难道都必须去逐个验证？可否有方法，比如类似F统计量的联合检验？
2 我要求对一个总体中N任意取出m个序列的等权新组合、再减去另一对照序列差额的方差，m从1可以到N，任意取m构成的新序列共有N!/m!(N-m)!个，也共有N!/m!(N-m)!个差额方差，求这N!/m!(N-m)!个差额方差的平均值，验证差额方差平均值是否显著异于零。
你的解似乎偏离了！！

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：对于别人热情的参与拟的帖子要给予评分

7# zhangtao


不好意思，因为我没有楼主的数据，只是在逻辑上给出了程序得思路，要有相关数据才能验证，可能会存在程序的语法错误。

8# 楚韵荆风
因为操作显示有问题，所以还在考虑怎样修改！本打算能做出合理结果就给你奖励的！好吧，评分支持一下！！