|
发表于 2022-9-10 10:36:21
|
显示全部楼层
你说的用R语言修复已损坏的FCS文件,就是之前我在论坛上帮Smile14站友解决的那件事:
收样结束后电脑死机,重启导致文件读不到数据
https://www.flowcyto.cn/bbs/thread-7600-1-1.html
(出处: 流式中文网)
R语言调用FlowCore是可以将FCS数据以表格形式展示出来的,但是将表格式的数据,再写回FCS文件,相当于就是把表格数据转换为FCS。
请将原来的mIHC中的数据,复制粘贴到Excel中,另存为csv文件。
接着在R语言中,运行下面这段代码:
library('flowCore')
library('Biobase')
library('data.table')
tempdata <- fread("这里填写你的CSV文件路径,windows路径\要用/替代", check.names = FALSE)
data_subset <- rbindlist(as.list(tempdata))
metadata <- data.frame(name=dimnames(data_subset)[[2]],desc=paste('column',dimnames(data_subset)[[2]],'from dataset'))
## 设置FCS文件的metadata
metadata$minRange <- apply(data_subset,2,min)
metadata$maxRange <- apply(data_subset,2,max)
data_subset.ff <- new("flowFrame",exprs=as.matrix(data_subset), parameters=AnnotatedDataFrame(metadata))
write.FCS(data_subset.ff, paste0(a, ".fcs")) |
|