Curvas amplificación RT-qPCR Meta Equipo NADP

if(!require(pacman))
  install.packages("pacman")

## Cargando paquete requerido: pacman

library("pacman")

p_load("vroom","dplyr","ggplot2","tidyr")

Curvas_amplif <- vroom(file="https://raw.githubusercontent.com/ManuelLaraMVZ/resultados_PCR_practica/refs/heads/main/Amplif_grupo2_17022025.csv")

## Rows: 40 Columns: 7
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## dbl (7): Cycle, A1, B1, C1, D1, E1, H1
## 
## ℹ 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_amplif2 <- Curvas_amplif %>% 
  mutate(Ciclos= Cycle, Manuel= A1, Regina=B1, Marian=C1, Ricardo=D1, Neto=E1, Johan=H1) %>% 
  select(-Cycle:-H1) %>% 
  filter(Ciclos>=10)
Curvas_amplif2

Curvas_amplif3 <- Curvas_amplif2 %>% 
  pivot_longer(cols=-Ciclos,
               names_to = "Muestras",
               values_to = "Fluorescencias")
Curvas_amplif3

Grafica_amplif <- ggplot(data = Curvas_amplif3, 
                         mapping = aes(x=Ciclos,
                                       y=Fluorescencias,
                                       color=Muestras)) +
  geom_line(linewidth=1.5)
Grafica_amplif

Umbral_ciclos <- 50
Grafica_amplif2 <- Grafica_amplif+
  geom_hline(yintercept = Umbral_ciclos, linetype="dashed", color="red")+
  theme_classic()+
  labs(title= "Curvas de amplificación RT-qPCR",
       subtitle = "Todas las muestras",
       caption = "Diseñó: NADP",
       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_amplif3$Ciclos),
                                  max(Curvas_amplif3$Ciclos),
                                  by=2))+
  scale_y_continuous(breaks = seq(min(Curvas_amplif3$Fluorescencias),
                                  max(Curvas_amplif3$Fluorescencias),
                                  by=250),
                     labels=scales::number_format(accuracy = 0.05))

## 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 every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

Grafica_amplif2

Grafica_amplif_equipo <- ggplot(Curvas_amplif2, 
                                aes(x = Ciclos)) +
  geom_line(aes(y = Neto, color = "NADP"), linewidth = 1.2) +
  geom_line(aes(y = Manuel, color = "Manuel"), linewidth = 1.5) +
  geom_hline(yintercept = Umbral_ciclos, linetype = "dashed", color = "red") +
  theme_classic() +
  labs(title = "Curvas de amplificación RT-qPCR",
       subtitle = "NADP, Tejido: Riñon",
       caption = "Diseñó: NADP",
       x = "Ciclos",
       y = "Fluorescencia (u.a)",
       color = "Equipo") +  # Etiqueta de la leyenda
  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_amplif3$Ciclos),
                                  max(Curvas_amplif3$Ciclos),
                                  by = 2)) +
  scale_y_continuous(breaks = seq(min(Curvas_amplif3$Fluorescencias),
                                  max(Curvas_amplif3$Fluorescencias),
                                  by = 250),
                     labels = scales::number_format(accuracy = 0.05)) +
  scale_color_manual(values = c("NADP" = "red", "Manuel" = "green"))  # Asignar colores manualmente

Grafica_amplif_equipo