How many transcripts?
Remember in the video how we discovered the length of the exons in one of the transcripts of our gene of interest? It is your turn to find out how many transcripts the gene ABCD1 has. You can find out by using:
transcriptsBy(x, by = "gene")
Once you get all the transcripts by gene, you can then select any gene by its id. The gene id of ABCD1 is 215. A little tip to select the information on a specific gene is to use back-tick marks around the gene id, for example transcripts$`215`.
Deze oefening maakt deel uit van de cursus
Introduction to Bioconductor in R
Oefeninstructies
- Load the human transcripts DB to
hg. - Pre-filter chromosome X using
seqlevels(). - Get all transcripts in
hgby"gene", store results inhg_chrXt, and print it. - Select gene
`215`fromhg_chrXtusing$.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Load the human transcripts DB to hg
library(TxDb.Hsapiens.UCSC.hg38.knownGene)
hg <- ___
# Prefilter chromosome X "chrX" using seqlevels()
___(___) <- c(___)
# Get all transcripts by gene and print it
hg_chrXt <- transcriptsBy(___, by = ___)
hg_chrXt
# Select gene `215` from the hg_chrXt
___$___