There are 2 ways to get the scripts. If you are familiar with command lines, you can clone the repository where you want:
git clone https://github.com/lldelisle/rnaseq_rscripts.git
Alternatively, you can simply download the repository:
Read the documentation about this file here. Do not forget the 2 columns hts_count_file
and cufflinks_file
.
Modify the file example/configFileRNAseq_step1.R
to fit what you want and what you need. You can have more details on each parameter here
Open the file step1-generateTables.R
in RStudio and source and follow instructions in the console of RStudio.
I know it is counter intuitive but I suggest to start with step3 to see how your samples behave (which one is close to which one).
Again, first modify the config file and use your FPKM table subset, and then source the step3-grapfClusteringPCAGenes.R
.
If you have only 2 conditions, run the step2. If you have more than 1 comparison, I propose to use step2-multi.
For both, modify the config file, and then source step2-DESeq2.R
or step2-multi_DESeq2.R
.