我正在使用下面的循环获取一个串联的表,并将其平均分割,并将该部分附加到flowSet对象内的flowFrames上。到目前为止,这对我的数据非常有用,因为我一直需要查看> 50000个单元格。现在,我有了一个没有统一数量的单元格的数据集,而且我负担不起从分析中删除/删除那些样本。
cellnum<-50000 #change this to number of cells per file
start_cell<-1
end_cell<-cellnum
for (i in 1:length(d_flowSet)) {
print(start_cell)
print(end_cell)
app.dataSmall<-app.data[start_cell:end_cell,]
new_frame <- fr_append_cols(d_flowSet[[i]], app.dataSmall)
filename_fcs <- filenames[i]
Sys.sleep(1)
write.FCS(new_frame, filename = filename_fcs,what="numeric", delimiter = "|", endian="big")
end_cell<-end_cell+cellnum
start_cell<-start_cell+cellnum
}
我想做的是更改循环,以便从称为eventCount
的列表对象中获取单元格的数量。
我该怎么做?如果您能解释为什么那也很好:)