我有一个RMA规范化数据(来自CEL文件),并想将其写入可以在excel中打开但有一些问题的文件中。
library(affy)
cel <- ReadAffy()
pre<-rma(cel)
write.table(pre, file="norm.txt", sep="\t")
write.table(pre, file="norma.txt")
输出是在使用上述命令写入的文本文件中按行排列的,因此当导出到excel时,输出格式错误,并且随着最大行数的使用,许多信息将被切断。以下方式:
GSM 133971.CEL 5.85302 3.54678 6.57648 9.45634
GSM 133972.CEL 4.65784 3.64578 3.54213 7.89566
GSM 133973.CEL 6.78543 3.54623 2.54345 7.89767
如何以正确的格式将R中的CEL文件写入记事本或excel?
最佳答案
您需要使用exprs
函数从标准化探针中提取值。就像是:
write.csv(exprs(pre), file="output.csv", row.names=FALSE)
应该可以。