carb168<-c(6.9,7,4.7,4.6,5.5,5.5,5.6,5.7,6,6.1,6.7,6.8) # A. carbonarius data flavus168<-c(6.1,6.5,5.1,5,5.5,5.5,5.6,5.7,5.5,5.75,6.1,6.2) # A. flavus data results<-c(carb168,flavus168) one_run_concentration<-c(0,0,450,450,300,300,250,250,100,100,50,50) garlic_conc_ppm<-c(one_run_concentration,one_run_concentration) species<-c(rep("carbonarius",12),rep("flavus",12)) data<-as.data.frame(cbind(as.numeric(results), as.numeric(garlic_conc_ppm), as.factor(species))) plot(data$V2,data$V1,col=as.numeric(as.factor(data$V3)),pch=as.numeric(as.factor(species)))