#Aspergillus moulds

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_conc<-c(0,0,450,450,300,300,250,250,100,100,50,50)
garlic_conc_ppm<-c(one_run_conc,one_run_conc)
type<-c(rep("carbonarius",12),rep("flavus",12))
data<-as.data.frame(cbind(as.numeric(results), as.numeric(garlic_conc_ppm), as.factor(type)))
data
model1<-lm(data[,1] ~ data[,2]*data[,3])
par(mfrow=c(2,2))
plot(model1)

model4<-lm(data[,1] ~ log((data[,2]+1)*data[,3])) 
plot(model4)

summary(model4)