diff --git a/_reader/RNAseqWorkshop.pdf b/_reader/RNAseqWorkshop.pdf index 2d6b4854ac4f3f593ab6ac8fbfdd38f122c47956..4226a0ca38dd8de51b23e8a3b8de6d06207dc19d 100644 Binary files a/_reader/RNAseqWorkshop.pdf and b/_reader/RNAseqWorkshop.pdf differ diff --git a/workflows/sleuth_differential_expression.R b/workflows/sleuth_differential_expression.R index 02c7e1ad21cf067088703967ce033b41c2b36ca6..a5d5f058d582d61fd3c64067d9638cc5f4d48081 100644 --- a/workflows/sleuth_differential_expression.R +++ b/workflows/sleuth_differential_expression.R @@ -4,10 +4,6 @@ library(sleuth) library(ggplot2) # First we need to specify where the kallisto results are stored. -# If you didn't specify this in your kallisto script, move all kallisto results -# folders (one for each sample) by GUI or the command line into a new folder called -# "kallisto_results". - # Begin by storing the base directory of the kallisto results in a variable base_dir <- "runs/kallisto_results/" @@ -77,8 +73,8 @@ table(treatment.vs.mock$qval <= 0.01) head(treatment.vs.mock) # <<< challenge excercises >>> # -# 1. compare the logFC edgeR calculated to that which we did -# 2. where does the difference comes from? (it's in the edgeR manual) +# 1. compare the logFC sleuth calculated to that which we did +# 2. where does the difference comes from? # now we transfer the result to our compilation data.frame 'dfr' # actually, all we really want is the 'false discovery rate' AKA 'q_value'