setwd("data/")
CG1 = read.delim('Genome_butyrate_both.txt', header = F)
CG2 = read.delim('Genome_butyrate_either.txt', header = F)
CG3 = read.delim('Genome_butyrate_Ato.txt', header = F)
CG4 = read.delim('Genome_butyrate_buk_ptb.txt', header = F)
CG5 = read.delim('Genome_butyrate_but.txt', header = F)
CG = merge(CG1,CG2,by='V1',all=T)
colnames(CG)=c('Genome','both','either')
CG3 = merge(CG3,CG4,by='V1',all=T)
colnames(CG3)=c('V1','1','2')
CG3 = merge(CG3,CG5,by='V1',all=T)
colnames(CG3)=c('Genome','Ato','buk_ptb','but')
#CG = merge(CG,CG3,by='Genome',all=T)
CG$difference = CG$V2.x - CG$V2.y
table(CG$difference)
table(CG$both)
table(CG$either)
table(CG3$Ato)
table(CG3$buk_ptb)
table(CG3$but)
View(CG)
View(CG3)
nrow(CG)
nrow(C3)
nrow(CG3)
View(CG3)
View(CG3)
View(CG1)
View(CG2)
CG = merge(CG1,CG2,by='V1',all=T)
length(unique(CG$V1))
View(CG)
View(CG1)
View(CG2)
CG1=CG1[!duplicated(CG1$V1),]
write.table(CG1,'Genome_butyrate_both.txt',row.names = F,
col.names = F,quote=F,sep='\t')
CG2=CG2[!duplicated(CG2$V1),]
write.table(CG2,'Genome_butyrate_either.txt',row.names = F,
col.names = F,quote=F,sep='\t')
CG2 = read.delim('Genome_butyrate_Ato.txt', header = F)
CG4 = read.delim('Genome_butyrate_buk_ptb.txt', header = F)
CG5 = read.delim('Genome_butyrate_but.txt', header = F)
CG2=CG2[!duplicated(CG2$V1),]
write.table(CG2,'Genome_butyrate_Ato.txt',row.names = F,
col.names = F,quote=F,sep='\t')
CG2 = read.delim('Genome_butyrate_buk_ptb.txt', header = F)
CG5 = read.delim('Genome_butyrate_but.txt', header = F)
CG2=CG2[!duplicated(CG2$V1),]
write.table(CG2,'Genome_butyrate_buk_ptb.txt',row.names = F,
col.names = F,quote=F,sep='\t')
CG2 = read.delim('Genome_butyrate_but.txt', header = F)
CG2=CG2[!duplicated(CG2$V1),]
write.table(CG2,'Genome_butyrate_but.txt',row.names = F,
col.names = F,quote=F,sep='\t')
CG1 = read.delim('Genome_butyrate_both.txt', header = F)
CG2 = read.delim('Genome_butyrate_either.txt', header = F)
CG3 = read.delim('Genome_butyrate_Ato.txt', header = F)
CG4 = read.delim('Genome_butyrate_buk_ptb.txt', header = F)
CG5 = read.delim('Genome_butyrate_but.txt', header = F)
CG = merge(CG1,CG2,by='V1',all=T)
colnames(CG)=c('Genome','both','either')
CG3 = merge(CG3,CG4,by='V1',all=T)
colnames(CG3)=c('V1','1','2')
CG3 = merge(CG3,CG5,by='V1',all=T)
colnames(CG3)=c('Genome','Ato','buk_ptb','but')
CG = merge(CG,CG3,by='Genome',all=T)
table(CG$difference)
table(CG$both)
table(CG$either)
table(CG$Ato)
table(CG$buk_ptb)
table(CG$but)
CG$difference = CG$both - CG$either
table(CG$difference)
table(CG$both)
table(CG$either)
table(CG$Ato)
table(CG$buk_ptb)
table(CG$but)
View(CG)
write.table(CG,'Butyrate_summary.txt',row.names = F,
col.names = F,quote=F,sep='\t')
