Resolución del problema de Cinética de la Lactasa # 11.4 . Doran

Script en R 

# CINETICA ENZIMATICA DE LA LACTASA 
# Problema 11.4 - Principios de la ingeniería de los bioprocesos - Pauline M. Doran. Ed Acribia
# Dr. Cid Ramón González       22/05/2015      IBQ - ITSA


# Datos de experimentación (generación de vectores):
# Concentración de lactosa= S (mol*L^-1*10^2)
S=c(2.5e-2,2.27e-2,1.84e-2,1.35e-2,1.25e-2,0.72e-2,0.46e-2,0.204e-2)    #(mol*L^-1)
#Velocidad inicial de reacción = v (mol*L^-1*min^-1*10^3)
v=c(1.94,1.91,1.85,1.8,1.78,1.46,1.17,0.779)
v1=v*10^-3                                                              # mol*L^-1*min^-1

# Ecuación linear de Langmuir para resolver Km y Vmax
# s/v= (1/Vmax)*(s)+(Km/Vmax)

# Resolución:
sv=S/v1                                  # 1) Calcular s/v
reg=lm(sv~S)                          # 2) Genera los coeficientes de la regresión linear
summary(reg)                        # Visualiza los valores del modelo linear
c=coef(reg)                           # Extrae los coeficientes de la regresión linear
b=signif(c[1])                        # Asigna valor a "b"(Intercept) a partir de "c"
m=signif(c[2])                       # Asigna valor a "m"(pendiente) a partir de "c"
Vmax=1/m                           # mol*L^-1*min^-1
Km=Vmax*b                          # mol*L^-1

# Graficar usando "plot"- "s/v" contra "s"
plot(S,sv, col="black",  
     main = "Cinética de Lactasa (s/v= (1/Vmax)*(s)+(Km/Vmax))",       # título del gráfico
     abline(lm(sv~S), col="red", lty="dashed"),                                       # linea de regresión
     xlab="[Lactosa]",                                                                          # título eje de las x
     ylab="Velocidad/[Lactosa]")                                                          # título eje de las y 
text(x=0.0075,y=10, labels="R2=0.9985\ny=445.1047*x+1.6959")         # agregar R^2 y ecuación
# Fin del Script