R操作合并数据等

library(data.table)
setwd("/mnt/10t/mzy/48samples/16s/salmon/")
fread("DNAabd138.csv")->DNA
as.data.frame(DNA)->DNA
row.names(DNA)<-DNA[,1]
unique(DNA[,1])->namelist
x<-c("")
for (i in namelist){
t(data.frame(colSums(DNA[which(DNA$Name==i),3:dim(DNA)[2]]),stringsAsFactors = F))->tmp
dim(tmp)
row.names(tmp)[1]<-i

rbind(x,tmp)->x
}
x<-x[-1,]
write.table(x,file="tmp.unique",quote=F,row.names = T, col.names=T,sep=",")
read.table(file="tmp.unique",header = T,sep=",",stringsAsFactors = F)->unique #我也不知道这个怎么搞，上面合成的数据集不能运算，是字符型的数字，只好写出去再读进来
unique[-which(rowSums(unique)<1),]->unique
c<-c("")
for (i in 1:dim(unique)[1]){
  if (sum(unique[i,]<1 & unique[i,]>0)==0 & sum(unique[i,]==0)<= 12){c<-c(c,i)}
  }
c=c[-1]
subset(unique,row(unique)[,1]%in%c)->unique
row.names(unique)->unique$names
melt(unique,id.vars="names",variable.name="metagenome",value.name="TPM")->unique
info48 <- readxl::read_excel("/mnt/10t/mzy/48samples/info48.xlsx", 
          col_types = c("numeric", "numeric", "numeric", 
          "text", "text", "text", "date", "date", "numeric"))
paste("X",info48$metagenome,sep="")->info48$metagenome
paste("X",info48$metatranscriptome,sep="")->info48$metatranscriptome
head(info48)
head(unique)
row.names(unique)=unique$metagenome
merge(unique,info48,by.x="metagenome",all.x=T)->merged
silvaSSU138.ok.tax <- read.delim("/mnt/10t/database/silva/silvaSSU138.ok.tax.txt", header=FALSE)
merge(merged,silvaSSU138.ok.tax,by.x="names",by.y="V1",all.x=T)->merged
head(merged)
write.csv(merged,file="DNAmerged.csv",col.names = T,row.names = F,quote=F)