Skip to content

Commit 8cfb217

Browse files
committed
Finalised clusterProfiler
1 parent 0b90d6b commit 8cfb217

File tree

2 files changed

+24
-18
lines changed

2 files changed

+24
-18
lines changed

11-gprofiler2.Rmd

Lines changed: 20 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -37,33 +37,39 @@ Since we are doing ORA, we will need a filtered gene list, and a background gene
3737

3838
<p>&nbsp;</p> <!-- insert blank line -->
3939

40-
## End of activity summary
41-
42-
- We have extracted a gene list and background gene list from a DE dataset and run ORA with `gprofiler2` `gost` function
43-
- We have plotted the data with `gostplot` Manhattan plots and `ggplot2` dotplots
44-
- We have run a `gost` multi-query separating up and down regulated genes
45-
- We have verified that `gprofiler2` results match the results from `g:Profiler` web
46-
- We have captured all version details relevant to the session within the R notebook
47-
48-
<p>&nbsp;</p> <!-- insert blank line -->
49-
50-
The last task is to `knit` the notebook. Our notebook is editable, and can be changed. Deleting code deletes the output, so we could lose valuable details. If we knit the notebook to HTML, we have a permanent static copy of the work.
51-
40+
**To finish up:** The last task is to `knit` the notebook. Our notebook is editable, and can be changed. Deleting code deletes the output, so we could lose valuable details. If we knit the notebook to HTML, we have a permanent static copy of the work.
5241

5342
&#x27A4; Knit the notebook to HTML
5443

5544
On the `editor` pane toolbar, under `Preview`, select `Knit to HTML`.
5645

5746
If you have already run `Preview`, you will see `Knit` instead of `Preview`.
5847

59-
The HTML file will be saved in the same directory as the notebook, and with the same filename, but the `.Rmd` prefix will be replaced by `.html`. The knit HTML will typically open automatically once complete. If you receive a popup blocker error, click `cancel`, and in the `Files` pane of RStudio, single click the `gprofiler.html` file and select `View in Web Browser`.
48+
The HTML file will be saved in the same directory as the notebook, and with the same filename, but the `.Rmd` prefix will be replaced by `.html`. The knit HTML will typically open automatically once complete.
49+
50+
❗If you receive a popup blocker error, click `cancel` (not `try again`), and in the `Files` pane of RStudio, single click the `gprofiler.html` file and select `View in Web Browser`.
6051

6152
<img src="images/open-html-rstudio.png" style="border: none; box-shadow: none; background: none; ">
6253

6354
<p>&nbsp;</p> <!-- insert blank line -->
6455

6556
Note that the notebook will only successsfully knit if there are <span style="color: red;">no errors</span> in the code. You can 'preview' HTML with code errors.
6657

67-
**If your knit fails, please ask for assistance resolving the errors**
58+
**If your knit fails, please ask for assistance resolving the errors**
59+
60+
## End of activity summary
61+
62+
- We have extracted a gene list and background gene list from a DE dataset and run ORA with `gprofiler2` `gost` function
63+
- We have plotted the data with `gostplot` Manhattan plots and `ggplot2` dotplots
64+
- We have run a `gost` multi-query separating up and down regulated genes
65+
- We have verified that `gprofiler2` results match the results from `g:Profiler` web
66+
- We have saved a rendered HTML document as a static record of our analysis including all code executed, parameters applied, as well as version information for R, RStudio, R packages and gProfiler database to ensure reproducibility
67+
68+
<p>&nbsp;</p> <!-- insert blank line -->
69+
70+
71+
72+
73+
6874

6975

12-clusterprofiler.Rmd

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
# clusterProfiler and enrichplot
22

33

4-
[clusterProfiler](https://bioconductor.org/packages/release/bioc/html/clusterProfiler.html) is a comprehensive suite of enrichment tools. It has functions to run ORA or GSEA over commonly used databases (GO, KEGG, KEGG Modules, DAVID, Pathway Commons, WikiPathways) as well as universal enrichment functions to perform ORA or GSEA with custom gene sets. We will use these universal tools in the final activity of this workshop, focusing on the supported organisms and datbases for the present activity.
4+
[clusterProfiler](https://bioconductor.org/packages/release/bioc/html/clusterProfiler.html) is a comprehensive suite of enrichment tools. It has functions to run ORA or GSEA over commonly used databases (GO, KEGG, KEGG Modules, DAVID, Pathway Commons, WikiPathways) as well as universal enrichment functions to perform ORA or GSEA with novel species or custom gene sets. We will use these universal tools in the final activity of this workshop, focusing on the supported organisms and datbases for the present activity.
55

66
One of the key advantages of using R over web tools is flexibility with visualisations.
77

@@ -46,7 +46,7 @@ Since we have covered ORA with `gprofiler`, we will perform a GSEA with `cluster
4646
3. Extract the gene IDs and sort by log2 fold change to create the GSEA ranked gene list R object
4747
4. Use `bitr` to convert gene IDs from ENSEMBL to ENTREZ for compatability with `gseKEGG`
4848
5. Perform GSEA with `gseKEGG`
49-
6. Visualise results with `enrichplot`
49+
6. Visualise results with many different plot types from `enrichplot`
5050

5151
<p>&nbsp;</p> <!-- insert blank line -->
5252

@@ -69,10 +69,10 @@ Let's head over to RStudio now and try out some functions! 🏃
6969

7070
## End of activity summary
7171

72-
- We have explored the supported organisms and namespaces of the `clusterProfiler` enrichment functions
72+
- We have explored the supported organisms, namespaces and databases of the `clusterProfiler` enrichment functions
7373
- We have extracted a ranked gene list for GSEA and converted the gene IDs for compatability with `gseKEGG`
7474
- We have performed GSEA on the KEGG database with `gseKEGG` and visualised the results with multiple plot types
75-
- We have captured all version details relevant to the session within the R notebook
75+
- We have captured all version details relevant to the session within the R notebook and knit the file to HTML for record keeping
7676

7777
<p>&nbsp;</p> <!-- insert blank line -->
7878

0 commit comments

Comments
 (0)