Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Exercise: Pull out a subset of the genes using a p_value cutoff of 0.01. Generate a heatmap for just those genes.

Hint
Expand
Solution
Solution
Expand
Code Block
titleOne possible solution

gene_diff_data <- diffData(genes(cuff_data))
sig_gene_data <- subset(gene_diff_data, (p_value < 0.01))
head(sig_gene_data)
sig_geneids <- c(sig_gene_data1$gene_id)

myGenes <- getGenes(cuff_data, sig_geneids)
csHeatmap(myGenes, cluster = "both")

Hint

If you have trouble sourcing cummeRbund, try this:
module swap gcc intel
Reenter R and redo above steps.