modified dturtle.R in R

){

	biocpar <- BiocParallel::MulticoreParam(getDTthreads())

	tx2gene <- tx2gene(organism)

	names(ctl_count_files) <-	gsub("/quant.sf","",ctl_count_files) |>

								basename()

					posthoc = 0.1

				)

	dturtle$dge_analysis <- c(ctl_count_files,obs_count_files) |> 

							import_dge_counts(

								type = "salmon",

								tx2gene = tx2gene[,c("transcript_id", "gene_name")]

							) |> 

	dturtle$dge_analysis <- cts |> 

							run_deseq2(

								pd = pd,

								pd = data.frame(pd),

								id_col = "id",

								cond_col = "group",

								cond_levels = c(obs,ctl),