domingo, 19 de abril de 2015

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: 

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

No hay comentarios.:

Publicar un comentario