Script Clase introduccion a R Cinética enzimática

Ejercicio: 

La lactasa, también conocida como B-galactosidasa, cataliza la hidrólisis de la lactosa para producir glucosa y galactosa a partir de la leche y del suero. Se realizaron experimentos para calcular los parámetros cinéticos para la enzima. Los datos de velocidad inicial se muestran a continuación.

Concentración de lactosa (mol*L-1 * 102)	Velocidad inicial (mol*L-1*m-1 *103)
2.5	1.94
2.27	1.91
1.84	1.85
1.35	1.8
1.25	1.78
0.73	1.46
0.46	1.17
0.204	0.779

      Calcular Vmax y Km.

Script in "R"

# Asignar Valores de la concentración de Lactosa-sustrato (S) y de la velocidad de reacción (V)
S<-c(2.5,2.27,1.84,1.35,1.25,0.73,0.46,0.204)
V<-c(1.94,1.91,1.85,1.8,1.78,1.46,1.17,0.779)
iS=1/S # Devuelve el valor inverso de la concentracion del sustrato
iV=1/V # Devuelve el valor inverso de la velocidad de reacción
# Graficar los valores de los inversos del sustrato contra velocidad
plot(iS,iV,type="p",col="blue",main="Cinética lactasa",sub="Linearización",xlab="1/[lactasa]",ylab="1/Velocidad")
abline(lm(iV~iS), col="red")
 #Calcular los coefficientes de la regresion
  lm(iV~iS)
 Pendiente<-coef(lm(iV~iS))[2]
Intercept<-coef(lm(iV~iS))[1]
# Aplicando la linearizacion Lineweaver-Burk de la ecuación de Michaelis-Menten "v = Vmax / (1 + (Km/[S]))" entonces 
# "1 / v = 1 / Vmax + Km / Vmax x 1 / [S]". Entonces tenemos que: Intercept = 1/Vmax y pendiente=Km/Vmax por lo tanto:
# Vmax=1/Intercept y Km=Vmax*Pendiente,  asi:
Vmax=1/Intercept
Km=Vmax*Pendiente
print (Vmax)
print (Km)