[求助]用R写的这段代码对吗？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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 : 被替换的项目不是替换值长度的倍数

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：warnings warning Result Sample matrix 代码

你这里的t是矩阵，所以应该写成t[i,1]＝result或是把t改成向量试试！