Problems using facet_wrap about ggtranscript HOT 2 CLOSED

Ah yes! ok, it makes sense. I was actually almost there earlier but I couldn't figure out exactly what the problem was. I was trying to apply to_intron to only subsets of the dataframe but this is much easier. Thank you very much and great package, I was really missing something simple like ggtranscript, that works well with the ggplot2 philosophy.
Thanks for the prompt reply, I am closing this as it's resolved.
Andrea

from ggtranscript.

Thanks for the feedback on ggtranscript!

to_intron() assumes each group_var is represented once in the data. In your data each "feature_id" is repeated 3 times (one repeat per DAI), which results in the weird output.

The easiest solution is to create a unique group_var for each of your transcript + DAIs combinations. For example:

prop_plot <- prop_plot |> dplyr::mutate(unique_id = stringr::str_c(feature_id, "-", DAI)

Then, the rest of your code should work as intended.

# using the created "unique_id" as the group_var
prop_plot_introns <- ggtranscript::to_intron(prop_plot, "unique_id")

prop_plot %>%
  ggplot2::ggplot(ggplot2::aes(
    xstart = start,
    xend = end,
    y = feature_id
  )) +
  ggplot2::facet_grid(cols = vars(DAI)) +
  ggtranscript::geom_range(
    ggplot2::aes(fill = Proportion)
  ) +  
  ggtranscript::geom_intron(
    data = prop_plot_introns,
    ggplot2::aes(strand = strand)
  ) +
  ggplot2::scale_fill_gradient2(
    high = "darkorchid",
    low = "white",
    mid = "darkorange",
    midpoint = 0.5
    # midpoint = min(prop_plot$Proportion) + ((max(prop_plot$Proportion) - min(prop_plot$Proportion)) / 2)
  ) +
  ggplot2::guides(fill = ggplot2::guide_legend("Proportion"))

from ggtranscript.