| 191.712 | 200.815 | 215.146 | 249.401 | 217.788 | 212.715 |
| 174.938 | 169.162 | 203.946 | 196.35 | 210.391 | 179.661 |
| 228.01 | 202.719 | 245.253 | 280.987 | 265.614 | 223.254 |
| 487.342 | 446.748 | 494.293 | 601.477 | 562.337 | 479.632 |
| 581.103 | 573.794 | 634.982 | 734.85 | 729.229 | 633.034 |
| 1024.6 | 930.128 | 1098.17 | 1185.52 | 1232.88 | 1013.8 |
| 2214.55 | 2089.11 | 2291.14 | 2611.22 | 2700.72 | 2393.33 |
| 6827.95 | 6450.03 | 6961.99 | 7390.9 | 7580.08 | 6813.42 |
| 8718.63 | 8198.6 | 8944.41 | 9430.22 | 9446.58 | 8743.45 |
| 893.031 | 1096.28 | 1135.21 | 1818.05 | 1374.99 | 1252.97 |
| 1812.13 | 1990.45 | 2083.7 | 2992.22 | 2451.62 | 2221.81 |
| 2837.58 | 3058.13 | 3284.25 | 4118.02 | 3402.24 | 3215.58 |
| 119.767 | 132.518 | 153.262 | 212.91 | 158.974 | 147.066 |
| 237.259 | 271.16 | 289.426 | 382.05 | 305.85 | 277.902 |
| 370.406 | 404.413 | 468.395 | 640.628 | 481.54 | 454.303 |
| 215.427 | 247.536 | 269.567 | 409.154 | 300.616 | 273 |
| 216.853 | 248.486 | 285.843 | 390.737 | 312.847 | 268.478 |
| 346.921 | 388.559 | 429.21 | 566.067 | 425.542 | 387.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 : 被替换的项目不是替换值长度的倍数