全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
10764 1
2009-02-27

191.712200.815215.146249.401217.788212.715
174.938169.162203.946196.35210.391179.661
228.01202.719245.253280.987265.614223.254
487.342446.748494.293601.477562.337479.632
581.103573.794634.982734.85729.229633.034
1024.6930.1281098.171185.521232.881013.8
2214.552089.112291.142611.222700.722393.33
6827.956450.036961.997390.97580.086813.42
8718.638198.68944.419430.229446.588743.45
893.0311096.281135.211818.051374.991252.97
1812.131990.452083.72992.222451.622221.81
2837.583058.133284.254118.023402.243215.58
119.767132.518153.262212.91158.974147.066
237.259271.16289.426382.05305.85277.902
370.406404.413468.395640.628481.54454.303
215.427247.536269.567409.154300.616273
216.853248.486285.843390.737312.847268.478
346.921388.559429.21566.067425.542387.018

     .下面省略,一共1000行,6列

——————————————————————————————————————————————

这个是代码:出现下面的提示,得不到结果怎么回事,哪位好心人帮着解读下,感激不尽!

> m=3;n=3;
> Data=as.matrix(read.table(paste("data1.txt")))
> t<-matrix(nrow=1000,ncol=1)
> for(i in 1:1000){
+ rowData=Data[i,]
+ sample1=rowData[c(2:4)]
+ sample2=rowData[c(5:7)]
+ sam1m=mean(t(sample1))
+ sam2m=mean(t(sample2))
+ sam1v=var(t(sample1))
+ sam2v=var(t(sample2))
+ result=(sam1m-sam2m)/(sam1v/m+sam2v/n)^0.5
+ t=result
+ }
警告多于50个(用warnings()来显示第一个到第五十个)
> warnings()
警告信息:
1: In t = result : 被替换的项目不是替换值长度的倍数
2: In t = result : 被替换的项目不是替换值长度的倍数
3: In t = result : 被替换的项目不是替换值长度的倍数
4: In t = result : 被替换的项目不是替换值长度的倍数
5: In t = result : 被替换的项目不是替换值长度的倍数
6: In t = result : 被替换的项目不是替换值长度的倍数
7: In t = result : 被替换的项目不是替换值长度的倍数
8: In t = result : 被替换的项目不是替换值长度的倍数
9: In t = result : 被替换的项目不是替换值长度的倍数
10: In t = result : 被替换的项目不是替换值长度的倍数
11: In t = result : 被替换的项目不是替换值长度的倍数
12: In t = result : 被替换的项目不是替换值长度的倍数
13: In t = result : 被替换的项目不是替换值长度的倍数
14: In t = result : 被替换的项目不是替换值长度的倍数
15: In t = result : 被替换的项目不是替换值长度的倍数
16: In t = result : 被替换的项目不是替换值长度的倍数
17: In t = result : 被替换的项目不是替换值长度的倍数
18: In t = result : 被替换的项目不是替换值长度的倍数
19: In t = result : 被替换的项目不是替换值长度的倍数
20: In t = result : 被替换的项目不是替换值长度的倍数
21: In t = result : 被替换的项目不是替换值长度的倍数
22: In t = result : 被替换的项目不是替换值长度的倍数
23: In t = result : 被替换的项目不是替换值长度的倍数
24: In t = result : 被替换的项目不是替换值长度的倍数
25: In t = result : 被替换的项目不是替换值长度的倍数
26: In t = result : 被替换的项目不是替换值长度的倍数
27: In t = result : 被替换的项目不是替换值长度的倍数
28: In t = result : 被替换的项目不是替换值长度的倍数
29: In t = result : 被替换的项目不是替换值长度的倍数
30: In t = result : 被替换的项目不是替换值长度的倍数
31: In t = result : 被替换的项目不是替换值长度的倍数
32: In t = result : 被替换的项目不是替换值长度的倍数
33: In t = result : 被替换的项目不是替换值长度的倍数
34: In t = result : 被替换的项目不是替换值长度的倍数
35: In t = result : 被替换的项目不是替换值长度的倍数
36: In t = result : 被替换的项目不是替换值长度的倍数
37: In t = result : 被替换的项目不是替换值长度的倍数
38: In t = result : 被替换的项目不是替换值长度的倍数
39: In t = result : 被替换的项目不是替换值长度的倍数
40: In t = result : 被替换的项目不是替换值长度的倍数
41: In t = result : 被替换的项目不是替换值长度的倍数
42: In t = result : 被替换的项目不是替换值长度的倍数
43: In t = result : 被替换的项目不是替换值长度的倍数
44: In t = result : 被替换的项目不是替换值长度的倍数
45: In t = result : 被替换的项目不是替换值长度的倍数
46: In t = result : 被替换的项目不是替换值长度的倍数
47: In t = result : 被替换的项目不是替换值长度的倍数
48: In t = result : 被替换的项目不是替换值长度的倍数
49: In t = result : 被替换的项目不是替换值长度的倍数
50: In t = result : 被替换的项目不是替换值长度的倍数

二维码

扫码加我 拉你入群

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

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

全部回复
2009-3-1 19:42:00
你这里的t是矩阵,所以应该写成t[i,1]=result或是把t改成向量试试!
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群