# Use after the first script of competitive inhibition
# Script for double reciprocal (Lineweaver-Burk) plot
# Created by S. Fushinobu (2018.8.17)

Mtitledr <- "Competitive doble reciprocal" # main title of the plot

kic.dr <- data.frame(rconc=1/kic.df$conc,rrate=1/kic.df$rate,inhibitor=kic.df$inhibitor)
plot(kic.dr$rconc,kic.dr$rrate,pch="",xlim=c(-1/Km/(1+min(kic.dr$inhibitor)/Kic),max(kic.dr$rconc)),ylim=c(0,max(kic.dr$rrate)),main=Mtitledr,xlab=paste("1/[",Substname,"]"),ylab="1/v")

for (i in 1:length(inhibitor.conc)) {
 points(subset(kic.dr$rconc, kic.dr$inhibitor==inhibitor.conc[i]), subset(kic.dr$rrate, kic.dr$inhibitor==inhibitor.conc[i]), col=inhib.color[i], pch=1)
}

for (i in 1:length(inhibitor.conc)) {
  abline(1/Vmax, (1+(inhibitor.conc[i]/Kic))*Km/Vmax,col=inhib.color[i])
}

abline(v = 0)

# write legends
legend("topleft", paste(Inhname,labels,Inhunit), col=inhib.color, pch=1, bg="white")