R 循环覆盖

我有R版本3.1.0，正在使用Windows7。我正在尝试为包含474行的数据集编写循环。
对于每个唯一的ID值，我可以为tow列设置多个值。对于ID的所有唯一值，层和编号列都相同

我想用这个等式计算方差：

(sum((sall$tow-sall$strata)^2))/unique(sall$number)* (unique(sall$number)-1) 

对于所有唯一的ID值。从我在网上读到的内容来看，我的循环似乎正在覆盖自身。为了解决这个问题，我尝试添加最后四行代码，但仍然收到一条错误消息：replacement有10行，data有1行

以下是我的数据：

dput(sall) 

structure(list(ID = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 3L, 4L, 
4L, 5L, 5L, 5L, 6L, 6L, 7L, 7L, 7L, 8L, 8L, 9L, 9L, 9L, 10L, 
10L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 13L, 
13L, 14L, 14L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 
16L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 18L, 18L, 19L, 19L, 
19L, 19L, 19L, 19L, 19L, 19L, 20L, 20L, 20L, 20L, 20L, 21L, 21L, 
21L, 21L, 21L, 21L, 21L, 21L, 21L, 22L, 22L, 22L, 22L, 22L, 23L, 
23L, 23L, 23L, 23L, 23L, 24L, 24L, 24L, 24L, 24L, 24L, 24L), .Label = c("2007 Fall 0102", 
"2007 Fall 0201", "2007 Fall 0202", "2007 Fall 0301", "2007 Fall 0302", 
"2007 Fall 0401", "2007 Fall 0402", "2007 Fall 0501", "2007 Fall 0502", 
"2007 Fall 0601", "2007 Fall 0602", "2007 Fall 0701", "2007 Fall 0702", 
"2007 Fall 0801", "2007 Fall 0802", "2007 Fall 0901", "2007 Fall 0902", 
"2007 Fall 1001", "2007 Fall 1002", "2007 Fall 1101", "2007 Fall 1102", 
"2007 Fall 1201", "2007 Fall 1202", "2007 Fall 1301", "2007 Fall 1302", 
"2007 Fall 1401", "2007 Fall 1402", "2007 Fall 1501", "2007 Fall 1502", 
"2007 Fall BI03", "2007 Fall BI04", "2007 Fall RI03", "2007 Fall RI04", 
"2008 Fall 0102", "2008 Fall 0201", "2008 Fall 0202", "2008 Fall 0301", 
"2008 Fall 0302", "2008 Fall 0401", "2008 Fall 0402", "2008 Fall 0501", 
"2008 Fall 0502", "2008 Fall 0601", "2008 Fall 0602", "2008 Fall 0701", 
"2008 Fall 0702", "2008 Fall 0801", "2008 Fall 0802", "2008 Fall 0901", 
"2008 Fall 0902", "2008 Fall 0903", "2008 Fall 1001", "2008 Fall 1002", 
"2008 Fall 1101", "2008 Fall 1102", "2008 Fall 1201", "2008 Fall 1202", 
"2008 Fall 1301", "2008 Fall 1302", "2008 Fall 1401", "2008 Fall 1402", 
"2008 Fall 1501", "2008 Fall 1502", "2008 Fall BI03", "2008 Fall BI04", 
"2008 Fall RI03", "2008 Fall RI04", "2008 Spring 0102", "2008 Spring 0201", 
"2008 Spring 0202", "2008 Spring 0301", "2008 Spring 0302", "2008 Spring 0401", 
"2008 Spring 0402", "2008 Spring 0501", "2008 Spring 0502", "2008 Spring 0601", 
"2008 Spring 0602", "2008 Spring 0701", "2008 Spring 0702", "2008 Spring 0801", 
"2008 Spring 0802", "2008 Spring 0901", "2008 Spring 0902", "2008 Spring 1001", 
"2008 Spring 1002", "2008 Spring 1101", "2008 Spring 1102", "2008 Spring 1201", 
"2008 Spring 1202", "2008 Spring 1301", "2008 Spring 1302", "2008 Spring 1401", 
"2008 Spring 1402", "2008 Spring 1501", "2008 Spring 1502", "2008 Spring BI03", 
"2008 Spring BI04", "2008 Spring RI03", "2008 Spring RI04", "2009 Fall 0102", 
"2009 Fall 0201", "2009 Fall 0202", "2009 Fall 0301", "2009 Fall 0302", 
"2009 Fall 0401", "2009 Fall 0402", "2009 Fall 0501", "2009 Fall 0502", 
"2009 Fall 0601", "2009 Fall 0602", "2009 Fall 0701", "2009 Fall 0702", 
"2009 Fall 0801", "2009 Fall 0802", "2009 Fall 0901", "2009 Fall 0902", 
"2009 Fall 0903", "2009 Fall 1001", "2009 Fall 1002", "2009 Fall 1101", 
"2009 Fall 1102", "2009 Fall 1201", "2009 Fall 1202", "2009 Fall 1301", 
"2009 Fall 1302", "2009 Fall 1401", "2009 Fall 1402", "2009 Fall 1501", 
"2009 Fall 1502", "2009 Fall BI03", "2009 Fall BI04", "2009 Fall RI03", 
"2009 Fall RI04", "2009 Spring 0102", "2009 Spring 0201", "2009 Spring 0202", 
"2009 Spring 0301", "2009 Spring 0302", "2009 Spring 0401", "2009 Spring 0402", 
"2009 Spring 0501", "2009 Spring 0502", "2009 Spring 0601", "2009 Spring 0602", 
"2009 Spring 0701", "2009 Spring 0702", "2009 Spring 0801", "2009 Spring 0802", 
"2009 Spring 0901", "2009 Spring 0902", "2009 Spring 0903", "2009 Spring 1001", 
"2009 Spring 1002", "2009 Spring 1101", "2009 Spring 1102", "2009 Spring 1201", 
"2009 Spring 1202", "2009 Spring 1301", "2009 Spring 1302", "2009 Spring 1401", 
"2009 Spring 1402", "2009 Spring 1501", "2009 Spring 1502", "2009 Spring BI03", 
"2009 Spring BI04", "2009 Spring RI03", "2009 Spring RI04", "2010 Fall 0102", 
"2010 Fall 0201", "2010 Fall 0202", "2010 Fall 0301", "2010 Fall 0302", 
"2010 Fall 0401", "2010 Fall 0402", "2010 Fall 0501", "2010 Fall 0502", 
"2010 Fall 0601", "2010 Fall 0602", "2010 Fall 0701", "2010 Fall 0702", 
"2010 Fall 0801", "2010 Fall 0802", "2010 Fall 0901", "2010 Fall 0902", 
"2010 Fall 0903", "2010 Fall 1001", "2010 Fall 1002", "2010 Fall 1101", 
"2010 Fall 1102", "2010 Fall 1201", "2010 Fall 1202", "2010 Fall 1301", 
"2010 Fall 1302", "2010 Fall 1401", "2010 Fall 1402", "2010 Fall 1501", 
"2010 Fall 1502", "2010 Fall BI03", "2010 Fall BI04", "2010 Fall RI03", 
"2010 Fall RI04", "2010 Spring 0102", "2010 Spring 0201", "2010 Spring 0202", 
"2010 Spring 0301", "2010 Spring 0302", "2010 Spring 0401", "2010 Spring 0402", 
"2010 Spring 0501", "2010 Spring 0502", "2010 Spring 0601", "2010 Spring 0602", 
"2010 Spring 0701", "2010 Spring 0702", "2010 Spring 0801", "2010 Spring 0802", 
"2010 Spring 0901", "2010 Spring 0902", "2010 Spring 0903", "2010 Spring 1001", 
"2010 Spring 1002", "2010 Spring 1101", "2010 Spring 1102", "2010 Spring 1201", 
"2010 Spring 1202", "2010 Spring 1301", "2010 Spring 1302", "2010 Spring 1401", 
"2010 Spring 1402", "2010 Spring 1501", "2010 Spring 1502", "2010 Spring BI03", 
"2010 Spring BI04", "2010 Spring RI03", "2010 Spring RI04", "2011 Fall 0102", 
"2011 Fall 0201", "2011 Fall 0202", "2011 Fall 0301", "2011 Fall 0302", 
"2011 Fall 0401", "2011 Fall 0402", "2011 Fall 0501", "2011 Fall 0502", 
"2011 Fall 0601", "2011 Fall 0602", "2011 Fall 0701", "2011 Fall 0702", 
"2011 Fall 0801", "2011 Fall 0802", "2011 Fall 0901", "2011 Fall 0902", 
"2011 Fall 0903", "2011 Fall 1001", "2011 Fall 1002", "2011 Fall 1101", 
"2011 Fall 1102", "2011 Fall 1201", "2011 Fall 1202", "2011 Fall 1301", 
"2011 Fall 1302", "2011 Fall 1401", "2011 Fall 1402", "2011 Fall 1501", 
"2011 Fall 1502", "2011 Fall BI03", "2011 Fall BI04", "2011 Fall RI03", 
"2011 Fall RI04", "2011 Spring 0102", "2011 Spring 0201", "2011 Spring 0202", 
"2011 Spring 0301", "2011 Spring 0302", "2011 Spring 0401", "2011 Spring 0402", 
"2011 Spring 0501", "2011 Spring 0502", "2011 Spring 0601", "2011 Spring 0602", 
"2011 Spring 0701", "2011 Spring 0702", "2011 Spring 0801", "2011 Spring 0802", 
"2011 Spring 0901", "2011 Spring 0902", "2011 Spring 0903", "2011 Spring 1001", 
"2011 Spring 1002", "2011 Spring 1101", "2011 Spring 1102", "2011 Spring 1201", 
"2011 Spring 1202", "2011 Spring 1301", "2011 Spring 1302", "2011 Spring 1401", 
"2011 Spring 1402", "2011 Spring 1501", "2011 Spring 1502", "2011 Spring BI03", 
"2011 Spring BI04", "2011 Spring RI03", "2011 Spring RI04", "2012 Fall 0102", 
"2012 Fall 0201", "2012 Fall 0202", "2012 Fall 0301", "2012 Fall 0302", 
"2012 Fall 0401", "2012 Fall 0402", "2012 Fall 0501", "2012 Fall 0502", 
"2012 Fall 0601", "2012 Fall 0602", "2012 Fall 0701", "2012 Fall 0702", 
"2012 Fall 0801", "2012 Fall 0802", "2012 Fall 0901", "2012 Fall 0902", 
"2012 Fall 0903", "2012 Fall 1001", "2012 Fall 1002", "2012 Fall 1101", 
"2012 Fall 1102", "2012 Fall 1201", "2012 Fall 1202", "2012 Fall 1301", 
"2012 Fall 1302", "2012 Fall 1401", "2012 Fall 1402", "2012 Fall 1501", 
"2012 Fall 1502", "2012 Fall BI03", "2012 Fall BI04", "2012 Fall RI03", 
"2012 Fall RI04", "2012 Spring 0102", "2012 Spring 0201", "2012 Spring 0202", 
"2012 Spring 0301", "2012 Spring 0302", "2012 Spring 0401", "2012 Spring 0402", 
"2012 Spring 0501", "2012 Spring 0502", "2012 Spring 0601", "2012 Spring 0602", 
"2012 Spring 0701", "2012 Spring 0702", "2012 Spring 0801", "2012 Spring 0802", 
"2012 Spring 0901", "2012 Spring 0902", "2012 Spring 0903", "2012 Spring 1001", 
"2012 Spring 1002", "2012 Spring 1101", "2012 Spring 1102", "2012 Spring 1201", 
"2012 Spring 1202", "2012 Spring 1301", "2012 Spring 1302", "2012 Spring 1401", 
"2012 Spring 1402", "2012 Spring 1501", "2012 Spring 1502", "2012 Spring BI03", 
"2012 Spring BI04", "2012 Spring RI03", "2012 Spring RI04", "2013 Fall 0102", 
"2013 Fall 0201", "2013 Fall 0202", "2013 Fall 0301", "2013 Fall 0302", 
"2013 Fall 0401", "2013 Fall 0402", "2013 Fall 0501", "2013 Fall 0502", 
"2013 Fall 0601", "2013 Fall 0602", "2013 Fall 0701", "2013 Fall 0702", 
"2013 Fall 0801", "2013 Fall 0802", "2013 Fall 0901", "2013 Fall 0902", 
"2013 Fall 0903", "2013 Fall 1001", "2013 Fall 1002", "2013 Fall 1101", 
"2013 Fall 1102", "2013 Fall 1201", "2013 Fall 1202", "2013 Fall 1301", 
"2013 Fall 1302", "2013 Fall 1401", "2013 Fall 1402", "2013 Fall 1501", 
"2013 Fall 1502", "2013 Fall BI03", "2013 Fall BI04", "2013 Fall RI03", 
"2013 Fall RI04", "2013 Spring 0102", "2013 Spring 0201", "2013 Spring 0202", 
"2013 Spring 0301", "2013 Spring 0302", "2013 Spring 0401", "2013 Spring 0402", 
"2013 Spring 0501", "2013 Spring 0502", "2013 Spring 0601", "2013 Spring 0602", 
"2013 Spring 0701", "2013 Spring 0702", "2013 Spring 0801", "2013 Spring 0802", 
"2013 Spring 0901", "2013 Spring 0902", "2013 Spring 0903", "2013 Spring 1001", 
"2013 Spring 1002", "2013 Spring 1101", "2013 Spring 1102", "2013 Spring 1201", 
"2013 Spring 1202", "2013 Spring 1301", "2013 Spring 1302", "2013 Spring 1401", 
"2013 Spring 1402", "2013 Spring 1501", "2013 Spring 1502", "2013 Spring BI03", 
"2013 Spring BI04", "2013 Spring RI03", "2013 Spring RI04", "2014 Spring 0102", 
"2014 Spring 0201", "2014 Spring 0202", "2014 Spring 0301", "2014 Spring 0302", 
"2014 Spring 0401", "2014 Spring 0402", "2014 Spring 0501", "2014 Spring 0502", 
"2014 Spring 0601", "2014 Spring 0602", "2014 Spring 0701", "2014 Spring 0702", 
"2014 Spring 0801", "2014 Spring 0802", "2014 Spring 0901", "2014 Spring 0902", 
"2014 Spring 0903", "2014 Spring 1001", "2014 Spring 1002", "2014 Spring 1101", 
"2014 Spring 1102", "2014 Spring 1201", "2014 Spring 1202", "2014 Spring 1301", 
"2014 Spring 1302", "2014 Spring 1401", "2014 Spring 1402", "2014 Spring 1501", 
"2014 Spring 1502", "2014 Spring BI03", "2014 Spring BI04", "2014 Spring RI03", 
"2014 Spring RI04"), class = "factor"), tow = c(0, 0, 0, 0, 0, 
0, 2.23605155512019, 0, 0, 0, 0, 0, 9.82973042531156, 0, 13.2468201825844, 
0.891911996342454, 5.21675816713178, 0, 0, 15.2033670694891, 
1.89572410870025, 48.7412108012422, 0, 0, 1.95759259768715, 0, 
3.88468129963282, 0.929188282053128, 1.91044228649375, 0, 3.00339331972541, 
15.6004700096039, 1.89642204148978, 3.81625917124503, 4.93924505216363, 
0.928742305487474, 32.559107185651, 0, 0.967770730644331, 0, 
1.88722551009165, 1.06020636196734, 0.958346832196337, 37.9680437072546, 
0, 2.73697536715327, 3.95587835585049, 0, 2.45847705019412, 0, 
0, 0, 28.8384678112108, 0, 0, 0, 6.45504724311346, 4.23466581629433, 
0, 0, 8.87505991044715, 0, 0.977534888388763, 0, 0, 13.1789290003378, 
0, 0, 0, 0, 0, 1.0244278082888, 0, 0, 0, 0, 0, 0, 0, 2.7408696865954, 
0, 0, 0, 2.13057544602497, 0, 0, 0.974322578911629, 1.01665761052166, 
0, 0, 0, 0, 1.92641091919418, 1.8606348312538, 4.52036131549582, 
0, 1.02462932795132, 0, 0, 0), strata = c(0, 0, 0, 0, 0.745350518373397, 
0.745350518373397, 0.745350518373397, 0, 0, 0, 0, 0, 4.91486521265578, 
4.91486521265578, 6.45183011535287, 6.45183011535287, 6.45183011535287, 
0, 0, 21.9467673264772, 21.9467673264772, 21.9467673264772, 0, 
0, 1.94742463243999, 1.94742463243999, 1.94742463243999, 1.46075597206807, 
1.46075597206807, 1.46075597206807, 1.46075597206807, 9.9567076276068, 
9.9567076276068, 9.9567076276068, 9.9567076276068, 9.9567076276068, 
9.9567076276068, 0.483885365322166, 0.483885365322166, 6.3729711112376, 
6.3729711112376, 6.3729711112376, 6.3729711112376, 6.3729711112376, 
6.3729711112376, 6.3729711112376, 1.60358885151115, 1.60358885151115, 
1.60358885151115, 1.60358885151115, 4.94102260882732, 4.94102260882732, 
4.94102260882732, 4.94102260882732, 4.94102260882732, 4.94102260882732, 
4.94102260882732, 4.94102260882732, 0, 0, 2.87894047489672, 2.87894047489672, 
2.87894047489672, 2.87894047489672, 2.87894047489672, 2.87894047489672, 
2.87894047489672, 2.87894047489672, 0.204885561657759, 0.204885561657759, 
0.204885561657759, 0.204885561657759, 0.204885561657759, 0.304541076288377, 
0.304541076288377, 0.304541076288377, 0.304541076288377, 0.304541076288377, 
0.304541076288377, 0.304541076288377, 0.304541076288377, 0.304541076288377, 
0.62097960498732, 0.62097960498732, 0.62097960498732, 0.62097960498732, 
0.62097960498732, 0.490511421619306, 0.490511421619306, 0.490511421619306, 
0.490511421619306, 0.490511421619306, 0.490511421619306, 1.05794649638585, 
1.05794649638585, 1.05794649638585, 1.05794649638585, 1.05794649638585, 
1.05794649638585, 1.05794649638585), number = c(2, 2, 2, 2, 3, 
3, 3, 2, 2, 3, 3, 3, 2, 2, 3, 3, 3, 2, 2, 3, 3, 3, 2, 2, 3, 3, 
3, 4, 4, 4, 4, 6, 6, 6, 6, 6, 6, 2, 2, 7, 7, 7, 7, 7, 7, 7, 4, 
4, 4, 4, 8, 8, 8, 8, 8, 8, 8, 8, 2, 2, 8, 8, 8, 8, 8, 8, 8, 8, 
5, 5, 5, 5, 5, 9, 9, 9, 9, 9, 9, 9, 9, 9, 5, 5, 5, 5, 5, 6, 6, 
6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7)), .Names = c("ID", "tow", "strata", 
"number"), row.names = c(NA, 100L), class = "data.frame")

我的循环尝试

#empty dataframe to hold loop output 
stratavar<-as.data.frame(matrix( 0, nrow=length(unique(one[,1])), ncol= 3 )) 

for(i in 1:length(unique(sall$ID))){ 

    new<-sall[sall$ID==i,] 
    stratavar[i,3]<-(sum((sall$tow-sall$strata)^2))/unique(sall$number)*
    (unique(sall$number)-1)   
    stratavar[i,1]<-as.character(unique(sall$ID)) 
    stratavar[i,2]<-unique(sall$strata) 

        if(i== "2007 Fall 0102" ) 
                newline<-data.frame(cbind(stratavar[i,1],stratavar[i,2],stratavar[i,3]))else{
                results<-data.frame(cbind(stratavar[i,1],stratavar[i,2],stratavar[i,3]))} 
                total<-rbind(newline,results) 
} 

#空数据帧以保持循环输出
stratavar您需要的是一个表，其中包含每个元素的名称和标识符（列名：id和name，按此顺序排列；它可以与其他数据一起位于输入表中）。在此之后，您可以编写以下内容：

#load the table with the name of the image & id 
TABLE<-read.table("Table.txt",header=T)

这一部分非常重要，您需要一个具有所需结构的文件，我指的是至少列名称

#saves XXX in c:/CSV.csv
for(n in 1:length(TABLE[,1])) {if (TABLE[n,1]==i){write.table(XXX,"c:/CSV.csv",sep=";", append=TRUE, col.names=FALSE)}}
 }

因此，您实际上必须粘贴
dput（sall）
的结果才能共享数据。。。另外，OmespacesinYourCode如何提高可读性？4个空格然后16个空格作为缩进宽度也是一个非常有趣的选择。@Dason缩进可能是我的错（或者是前一个引入任何代码格式的编辑器），但我鼓励OP修复它，我不确定我是否在跟踪你。我的数据是R中的一个合并数据帧，所以我需要导出它，然后再将其作为表导入。看起来我应该能够在R中写一个循环来处理数据。
#saves XXX in c:/CSV.csv
for(n in 1:length(TABLE[,1])) {if (TABLE[n,1]==i){write.table(XXX,"c:/CSV.csv",sep=";", append=TRUE, col.names=FALSE)}}
 }