Resolución de la constante de velocidad en una cinética de reacción de primer orden (R-Script)

# Método para calcular la constante de velocidad de una reacción de primer orden (k)

# Este ejemplo es una reacción de Bromo + ácido fórmico cuyos datos se encuentran en la tabla 13.1 del libro de Química  de Chang & Goldsby, undécima edición, p 567.

# Si la reacción es de primer grado entonces podemos comprobarlo con la ecuación ln[A]t=-kt+ln[A]0

# Paso 1: Importar los datos desde la tabla de excel en formato CSV usando:

R <- read.delim("~/Dropbox/ITSA/MATERIAS/CINETICA QUIMICO BIOLOGICA/R-exercises/Kinetics/reaccion_bromoac formico.txt")
# Cuando "R" corresponde a la tabla de datos:

View (R)

# Paso 2: Definir los valores de las variables a través de vectores: t=Tiempo; Br=concentracion del bromo

t<-R$Tiempo

Br<-R$Br

# Paso 3: Calcular el logaritmo de la concentración del Bromo

lnBr<-log(Br)

# Paso 4: Graficar el tiempo (t) contra logaritmo de la concentration de bromo (lnBr) en M

# para comprobar si la reacción es de primer orden, con la regresión lineal con el comando "abline" y "lm"

plot(t,lnBr, type="p", col="blue", xlab="Tiempo (s)", ylab="ln[Br](M)")

abline(lm(lnBr~t), col="red")
# Al cruzar la linea todos los puntos de concentración indica que la reacción es de primer orden. 

#  Paso 5: Realizar el cálculo de regresión lineal para encontrar los coeficientes de la recta:

Reg<-lm(lnBr~t)

summary(Reg)
# Donde la intersección será igual a la concentración inicial de Br. 
#  y coeficiente de t, es igual a la pendiente y por tanto igual a la constante de velocidad de reacción (k)

Video: 

https://youtu.be/mNHMSqe1oZY

Cualquier duda escribir en el espacio de comentarios en la parte de abajo. 

Saludos!