choose_tx: Select a subset of transcripts per junction that are more...

View source: R/choose_tx.R

choose_txR Documentation

Select a subset of transcripts per junction that are more likely to be affected by a junction.

Description

All possible affected transcripts are required and can be annotated with add_tx()

Usage

choose_tx(df)

Arguments

df

A data.frame with splice junctions in rows and at least the columns:

  • junc_id junction id consisting of genomic coordinates

  • tx_id transcript id of possibly affected transcripts

  • tx_lst a list of GRanges with the transcript

Value

A data.frame as with relevant transcript and junction combinations. If tx_id is NA in the input data.frame such rows are removed from the output data.frame.

This function selects transcripts that are more likely to be affected to reduce the amount of junction and transcript combinations. The function excludes transcripts for which both junction positions are located in an intron. Junctions in a given transcript must either represent an exon skipping, intron retention, exitron, or alternative splice site event or have both junction positions in an exon. Other junction-transcript combinations are also excluded. This function may loose relevant or keep irrelevant junction-transcripts in particular in regions with multiple isoforms with distinct splicing pattern.

See Also

add_tx

Examples

junc_df <- tibble::tibble(
  junc_id = c("chr2:152389996-152392205:-", "chr2:152389996-152390729:-",
              "chr2:152389955-152389956:-")
)

junc_df <- add_tx(junc_df, toy_transcripts)
choose_tx(junc_df)


TRON-Bioinformatics/splice2neo documentation built on July 27, 2024, 3:25 a.m.