Amplificación RT-qPCR
Lourdes Gallardo
2026-03-17
ctrl+alt+i Instalación de paquetes
if(!require(pacman))
install.packages("pacman")## Loading required package: pacmanLLmar a pacman
library("pacman")Lamar a paquteria
p_load("vroom",
"dplyr",
"ggplot2",
"tidyr")Llamar a base de datos
curvas_ampli <- vroom(file="https://raw.githubusercontent.com/ManuelLaraMVZ/resultados_PCR_practica/refs/heads/main/Ejemplo_amplificacion_RT-PCR.csv")## Rows: 35 Columns: 6
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## dbl (6): Cycle, A1, B1, C1, D1, E1
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.curvas_ampli## # A tibble: 35 × 6
## Cycle A1 B1 C1 D1 E1
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 1 0.0000213 -0.00531 0.0000279 0.000519 0.000468
## 2 2 -0.0000305 -0.00286 0.00179 0.000301 0.000462
## 3 3 -0.0000823 -0.00136 -0.000335 0.000842 0.0000266
## 4 4 -0.0000536 -0.00101 0.00166 0.000522 -0.000943
## 5 5 -0.0000565 -0.000787 0.000199 0.000586 0.0000877
## 6 6 0.000111 -0.000233 0.00231 0.000451 0.00076
## 7 7 0.000096 0.0000674 -0.000511 -0.00112 -0.000618
## 8 8 0.0000698 0.000602 -0.000595 -0.000944 0.00108
## 9 9 0.000202 0.00143 0.000166 0.000342 -0.000105
## 10 10 -0.00001 0.00101 -0.00136 -0.000254 0.00102
## # ℹ 25 more rowsModificar base de datos
curvas_ampli2<- curvas_ampli %>%
mutate(Ciclos = Cycle, H2O = A1, JRWF = B1, PPOF = C1, MLL= D1, DDR = E1) %>% select(-Cycle:-E1) %>%
filter(Ciclos>=17)
curvas_ampli2## # A tibble: 19 × 6
## Ciclos H2O JRWF PPOF MLL DDR
## <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 17 0.0000132 0.00126 -0.00121 0.000599 0.000786
## 2 18 0.0000394 0.000431 -0.000969 -0.00155 -0.00119
## 3 19 -0.000106 0.000227 -0.00145 -0.000874 -0.00113
## 4 20 -0.000104 0.000133 -0.00197 -0.000632 -0.00134
## 5 21 -0.000119 0.00104 0.000464 0.000245 -0.000417
## 6 22 -0.0000153 0.00028 0.0032 -0.00081 0.00158
## 7 23 0.0000704 -0.000086 0.00494 0.00162 0.000497
## 8 24 0.00012 -0.000238 0.0144 0.00139 0.00121
## 9 25 -0.0000899 -0.000484 0.0263 0.00634 0.00425
## 10 26 0.0000114 -0.00133 0.0491 0.0101 0.00894
## 11 27 -0.000104 -0.000824 0.0763 0.0202 0.0202
## 12 28 -0.000101 -0.000922 0.105 0.0365 0.0353
## 13 29 -0.000222 -0.000938 0.134 0.0624 0.0594
## 14 30 -0.0000675 -0.000187 0.162 0.0895 0.0823
## 15 31 -0.000089 -0.0000371 0.191 0.114 0.107
## 16 32 -0.000188 0.0018 0.220 0.141 0.128
## 17 33 -0.0000273 0.00345 0.243 0.170 0.152
## 18 34 0.000358 0.00615 0.268 0.197 0.173
## 19 35 0.000345 0.00913 0.293 0.223 0.192#reordenamiento
#1 seleccionar los datos a oredenar que son las variables dependientes
#2 agrupar por nombre
#3Agrupar las variables por ciclo
curvas_ampli3 <- curvas_ampli2 %>%
pivot_longer(cols=-Ciclos,
names_to = "muestras",
values_to = "fluorescencias")
curvas_ampli3## # A tibble: 95 × 3
## Ciclos muestras fluorescencias
## <dbl> <chr> <dbl>
## 1 17 H2O 0.0000132
## 2 17 JRWF 0.00126
## 3 17 PPOF -0.00121
## 4 17 MLL 0.000599
## 5 17 DDR 0.000786
## 6 18 H2O 0.0000394
## 7 18 JRWF 0.000431
## 8 18 PPOF -0.000969
## 9 18 MLL -0.00155
## 10 18 DDR -0.00119
## # ℹ 85 more rowsGrafica
grafica_ampli <- ggplot(data=curvas_ampli3,
mapping = aes(x=Ciclos, y= fluorescencias,
color=muestras))+
geom_line(linewidth =1.5)
grafica_ampli
Mejorar la grafica Quitar fondo mejorar ejes
umbral_ciclos <- 0.008
grafica_ampli2 <- grafica_ampli +
geom_hline(yintercept = umbral_ciclos, linetype = "dashed", color = "orange") +
theme_classic() +
labs(
title = "curvas de amplificacion RT-qPCR",
subtitle = "todas las muestras",
caption = "diseñó: Lourdes",
x = "Ciclos",
y = "Fluorescencia(u.a.)"
) +
theme(
axis.line = element_line(size = 1.2, color = "black"),
axis.title = element_text(face = "bold"),
axis.text = element_text(face = "bold"),
legend.title = element_text(face = "bold"),
legend.text = element_text(face = "bold")
) +
scale_x_continuous(
breaks = seq(min(curvas_ampli3$Ciclos), max(curvas_ampli3$Ciclos), by = 2)
) +
scale_y_continuous(
breaks = seq(min(curvas_ampli3$fluorescencias),
max(curvas_ampli3$fluorescencias),
by = 0.02),
labels = scales::number_format(accuracy = 0.02)
)## Warning: The `size` argument of `element_line()` is deprecated as of ggplot2 3.4.0.
## ℹ Please use the `linewidth` argument instead.
## This warning is displayed once per session.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.grafica_ampli2
Grafica de los datos del equipo
umbral_ciclos <- 0.008
grafica_ampli_equ <- ggplot(curvas_ampli2, aes(x = Ciclos)) +
# Línea de MLL
geom_line(aes(y = MLL, color = "MLL"), linewidth = 1.2) +
# Línea de H2O
geom_line(aes(y = H2O, color = "H2O"), linewidth = 1.5) +
# Línea horizontal del umbral
geom_hline(aes(yintercept = umbral_ciclos, linetype = "Umbral"), color = "orange") +
theme_classic() +
labs(
title = "curvas de amplificacion RT-qPCR",
subtitle = "todas las muestras",
caption = "diseñó: Lourdes",
x = "Ciclos",
y = "Fluorescencia(u.a.)",
color = "Muestra", # título de la leyenda de colores
linetype = "Referencia" # título de la leyenda de tipos de línea
) +
theme(
axis.line = element_line(size = 1.2, color = "black"),
axis.title = element_text(face = "bold"),
axis.text = element_text(face = "bold"),
legend.title = element_text(face = "bold"),
legend.text = element_text(face = "bold")
) +
scale_x_continuous(
breaks = seq(min(curvas_ampli3$Ciclos), max(curvas_ampli3$Ciclos), by = 2)
) +
scale_y_continuous(
breaks = seq(min(curvas_ampli3$fluorescencias),
max(curvas_ampli3$fluorescencias),
by = 0.02),
labels = scales::number_format(accuracy = 0.02)
) +
# Colores personalizados para cada línea
scale_color_manual(values = c("MLL" = "cyan3", "H2O" = "darkseagreen2")) +
# Tipo de línea para el umbral
scale_linetype_manual(values = c("Umbral" = "dashed"))
grafica_ampli_equ