这篇博客展示了如何使用R语言读取多个数据文件并进行数据处理。通过循环读取'HSC-result.txt'、'CMP-result.txt'等文件，并进行筛选操作，将结果保存为新的文件。接着，博主应用了Levene's test和t-test进行统计分析，生成了不同输出文件，如'mlpmep_levene_p.txt'和'mlpmep_ttest_pvalue.txt'，以存储测试结果。
hsc<-read.table('HSC-result.txt',header = T,row.names = 1)
cmp<-read.table('CMP-result.txt',header = T,row.names = 1)
gmp<-read.table('GMP-result.txt',header = T,row.names = 1)
mep<-read.table('MEP-result.txt',header = T,row.names = 1)
mlp<-read.table('MLP-result.txt',header = T,row.names = 1)
mpp<-read.table('MPP-result.txt',header = T,row.names = 1)
mlpmep<-cbind(mlp,mep)
#write.table(mlpmep,'mlp-mep-use2.txt',quote = F,sep = '\t')
mlpmep<-mlpmep[apply(mlpmep, 1, mean)>30,]
write.table(mlpmep,'mlp-mep-tpm-morethan30.txt',sep = '\t',quote = F)
yinzi1<-factor(c(rep(1,67),rep(2,224)))
get_len <- function(x){
a = leveneTest(as.numeric(x)~yinzi1)
return(a$'Pr(>F)')
}
filter_res<-apply(mlpmep,
最低0.47元/天 解锁文章
2万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
