@@ -93,13 +93,7 @@ tryCatch( {
9393 ggtitle("Number of sequences after each of the sequence assembly steps") +
9494 theme(axis.text.x= element_text(angle = 45))
9595
96- p <- ggplotly(seqs_plot) %>%
97- layout(shapes = list(list(type = "rect",
98- xref = "paper", yref = "paper",
99- x0 = -0.02, x1 = 1.02, y0 = -0.02, y1 = 1.02,
100- line = list(color = "black", width = 1.5),
101- fillcolor = 'rgba(0,0,0,0)',
102- layer = 'above')))
96+ p <- ggplotly(seqs_plot)
10397 p
10498
10599
@@ -130,6 +124,8 @@ if (is.null(tab_seqs_assembled$sample_id)) {
130124}
131125
132126dat <- tab_seqs_assembled %>%
127+ filter( !grepl("-fail.tsv", output) ) %>%
128+ filter( !grepl("productive-F.tsv", output) ) %>%
133129 tidyr::pivot_wider(id_cols=sample_id,
134130 names_from=task,
135131 values_from=output_size)
@@ -171,11 +167,7 @@ seqs_plot_assembled <- ggplot(data=tab_seqs_assembled,
171167 theme(axis.text.x = element_text(angle = 45)
172168 )
173169
174- p <- ggplotly(seqs_plot_assembled) |>
175- layout(
176- xaxis = list(showline = TRUE, linecolor = "black", linewidth = 0.66, mirror = TRUE),
177- yaxis = list(showline = TRUE, linecolor = "black", linewidth = 0.66, mirror = TRUE)
178- )
170+ p <- ggplotly(seqs_plot_assembled)
179171p
180172```
181173
@@ -218,11 +210,7 @@ g2 <- ggplot(family, aes(x=gene, y=seq_freq, fill=sample_id, group=sample_id)) +
218210 ylab("Frequency") +
219211 xlab("") +
220212 theme(legend.position = "right")
221- p <- ggplotly(g2) |>
222- layout(
223- xaxis = list(showline = TRUE, linecolor = "black", linewidth = 0.66, mirror = TRUE),
224- yaxis = list(showline = TRUE, linecolor = "black", linewidth = 0.66, mirror = TRUE)
225- )
213+ p <- ggplotly(g2)
226214p
227215
228216ggsave(filename = paste0(vfamily_dir, "/V_Family_distribution_patient.pdf"), plot = g2, width = 18, height = 15, units = "cm")
@@ -247,11 +235,7 @@ g2 <- ggplot(family, aes(x=gene, y=seq_freq, fill=sample_id, group=sample_id)) +
247235 theme(axis.text.x=element_text(angle=45, hjust=1, vjust=1)) +
248236 ylab("Frequency") +
249237 xlab("")
250- p <- ggplotly(g2) |>
251- layout(
252- xaxis = list(showline = TRUE, linecolor = "black", linewidth = 0.66, mirror = TRUE),
253- yaxis = list(showline = TRUE, linecolor = "black", linewidth = 0.66, mirror = TRUE)
254- )
238+ p <- ggplotly(g2)
255239p
256240ggsave(filename = paste0(vfamily_dir, "/V_gene_distribution_by_sequence_patient.pdf"), plot = g2, width = 20, height = 40, units = "cm")
257241ggsave(filename = paste0(vfamily_dir, "/V_gene_distribution_by_sequence_patient.png"), plot = g2, width = 20, height = 40, units = "cm")
0 commit comments