nair_pc1 = read.table('data/hNAD_PC1')
b_compartment = read.table('data/Other_PC1')
df = data.frame(pc1 = c(nair_pc1$V6,b_compartment$V6),group = c(rep('hNAD',dim(nair_pc1)[1]),
                                                                rep('non-hNAD',dim(b_compartment)[1])),
                dose = 'B compartment')
library(ggplot2)
bp <- ggplot(df, aes(x=dose, y=pc1, fill=group)) + 
  geom_boxplot(width=0.4)+scale_fill_manual(values=c("#2BA809", "#00BFC4"))+
  labs(title="",x="", y = "PC1")+ylim(-0.05,0)
bp + theme_classic()