...
The resultant plot is here.
Exercise 7: Pull out a subset of the genes using a log_fold_change greater than 1.5. Generate expression bar plots for just those genes.
...
| Code Block | ||
|---|---|---|
| ||
gene_diff_data <- diffData(genes(cuff_data))
sig_gene_data <- subset(gene_diff_data, (log2_fold_change > 1.5))
head(sig_gene_data)
sig_geneids <- c(sig_gene_data$gene_id)
pdf("expressionboxplot2.pdf")
myGenes <- getGenes(cuff_data, sig_geneids)
expressionBarplot(myGenes)
dev.off() |
The resultant graph is here.