From 1bd5e41c8faac19311238d809d86ec835319ae38 Mon Sep 17 00:00:00 2001 From: Elias Projahn Date: Wed, 22 Jun 2022 19:34:39 +0200 Subject: [PATCH] Add additional information page --- R/server.R | 2 ++ R/sysdata.rda | Bin 0 -> 3515 bytes R/ui.R | 31 +++++++++++++++++++++++++++---- scripts/gsea.R | 3 +++ 4 files changed, 32 insertions(+), 4 deletions(-) create mode 100644 R/sysdata.rda diff --git a/R/server.R b/R/server.R index 00600c3..18cd5d4 100644 --- a/R/server.R +++ b/R/server.R @@ -159,6 +159,8 @@ server <- function(input, output, session) { DT::formatRound("p_value", digits = 4) |> DT::formatPercentage(c("total_ratio", "query_ratio"), digits = 1) }) + + output$gsea_plot_ranking <- plotly::renderPlotly(gsea_plot_ranking) } #' Create a displayable data table from the gene results data. diff --git a/R/sysdata.rda b/R/sysdata.rda new file mode 100644 index 0000000000000000000000000000000000000000..c785d0cf64f0a8541f1f79d591fca193b2d9f600 GIT binary patch literal 3515 zcmV;s4Mg%nT4*^jL0KkKSzb3#qW~-D|NsC0|NsC0|NsC0|NsC0|M~E#5kZcCAfQkH zkVF6#PAlLLeEIOs*!F#Wb0~JP(4GYhW-8jyp|PmbAOHZM1xLg}DD@cw)X;iCrhouo zfEo<|13&`;X`!IR00EF>&@)5^Kn+3~O{t+zVJ0;-Jx?jM8%WwsGG>Gg8W@1c&|)wM z3<7BNGJ2UbdYKv+N2#OIdW;QIOeyIgGyu>v00TonGynhq00000000^Q000Rzl1fw2 zg&RaLCYc&GMAJqS1T--LriMlUhCsq-U<7DnUMn)ivjMURX zrqs<5qfM#goI*>!xC{kEb%Ve^JlFJH9C{kHrVPy&y6qG3` zC|M{L78Zq~#LK<4z{I5`-zrDRBV>z_#}6cRM>KXvXm=C3o#V|sNZle zmTsB==$X|9aUA-tb-aDas#^W`u1(u-MrAiA5=x9TaybnIdQkf0d`Gx3Y30hVO3Kq4M_4FMKGvo|O}Do`#&Bg~G_hv>q; z>zH{L*&Q9v$#1TfNrDtWnx)G-Fx?Q+G80{Kvl^JB5h5xBAcF*h7%>!#OhThKwi2>I z7!!s8`M_d>BN+)I$s|rBgxiugS-#*|Rm9VE^ETG!g|`2Lt_gnTN%7Qn2q%;nwgwBF z)8d&ToMF9Jo}fH{XSz#53o_maM@IuAQt7nUFXf4OsG*K!DUpmZp!6aowt7J#dI);2CaC1rI|S){94Pyr6b1!)L^p;DEJL5Tt#Kwd(u!DWr9iW21M z`5?hA3IaE&h5}(RNT9@9rFyMh2{zK@{kpFSeDkX6t458EZ9UxBnRTN`h|-Pr%~}0-Q^7zO6=t>qB{hb766$BY zq%%k2k$@nDRVh+{=!}ps83Z8-5m?R2;LJl2U!kr_;XnFlkQR8W^~=BnAP$>RxJ(y_ zgi9&r0|m*ZCVd+Bb^U&&v$2TgW_;wQSm|?Gd}lLiGrY~zX^?&g)zg!rW0E#u(cb5} z;oW9g!-^*M_Fv+<$=@zFCYl^)f*^mU;3DDi2ne(}E`UffNkG#<$9(L&6|Mx=a$}c0 zHmd?_Wy$eqon>WtOsMLsuJv6eEO4BwCAqO$QsVW@>ltoe>srdE$7nC2=2aPQWn8#m zZmK)gFR3=sKdysID$U|yUfuJf+Pta809O$P7soP%F;r>15_uqy6=(whbhc&{zg!L?L=a#{6Q?ZLXfV+F@`Au1J$nGg zpragYS|@X{K?%`mB0-AEU+*I%5d`gCk|k$inQ$@4*)ifAG75!4hg1y0GpuZ^xjjj0 zL!pVhRFqPsQmTpL!$M)9o@0gx`(J5k85&sKO3St5X*Do6)OxB@P_w^I21W^j_qAJ+ zn96%BMP%jQZtpV9s#csvEfsz2*RI}dR|}S`+awEcFB9+MpaS0AUZ@mX*i;nS{hybI z!GZ|Ji0=D-mBiz3wb9GFXSrBiSWw=K(?>2_|-i-Ruw=p#ICD6=+hWA86MI`%W8W@Wi4Egn60X=8P? z?O}^6MH>Yuq9hfF1O;NOA|gm4vW*lP5D{97Vib}bazK|95stD4#kAa8 zMp`JMiYaQTqR9eDCPU;z_hY`CP0e3 z@0Siqp%%I-H{0G%Yc{q?UA3_aXhm9bEflMfNs&)E6kV%g<_of->8hh)(WbVuGZrke zl@^_uV{}teiApJ15HZE1K~-?nVG&jkQIG)?kR(BZhyg(&3SUMnMoGqufmI;^Kw}hu z1kf)Jqr|2c2H{*(rc$?OPBP8H1T-KM7=#F;P@%fQ1mBF1l1v?tfkle4C?QA|;wY;~ zj3NvmA|Zf`N+d-DL{JDtgdoI85J?dPBNEI}vjQ=VD2Og7Frbf$MMMx_Qctr&eYymg zL4b@wARD3^$jHQF3@ovVR8THakrfGfM2b>El#vT?8*o~PSW%<~j7cRCf>=T3^FAnZCNRj<&3tMm2|5XwS9Bq zHo7ADkv*3~iJ$GIJGrbkPxj7u6qJ`u^y=HhBygE{h3(SX zZ4agvz>dh*fmX!KR@BOz_%fWbPC=4n6?*97FL^o+&@IJW1k?k>;Kw#U3czud<7( z!W#;Z1%!hcBtTfI2!kTTSg=?SAd!l|g8+iWipC=tEEo((tQJTC7@-y_L?Re*<6y!N zgn+R~oYV)TW!)QCz=69xP2pg}Kf6B!N}^R!N+}fyF-SThXj*Mz<(i{#EnYEPLstgyB{n^iH=ORicTuny=kZXfh z23neI;DZ}|YzFiZH~Ee+H5dWJK3fQ}PXUx9v0zDAFexB?vcqzsKNETh@WHH2)&#Lq zeZ$y(MbaH456$Fx>*Q*r;}ht9&B6EeA($OOVsUjStJZ{2OTCub)WDVkq4gzE(CbH=hGcbfEKj_2FjGC0wCy9wO{e3% zgWhQ6l#T}rUAF7Owj;jiyJu(L9OD)mambs~kI$*qk8P4ov*(V`(n@7Tl5Z1}aV(SZ zd!ALhokxcUH;%I8w*e`Qu~paJ!YQ8gW=2czX>aCbw|N(7(3c|Yp_v{F8Te7_^P%nw-It=_4m&y2Cig|VI|$V5;+*ckTQ?F%Q5yY6<63CdSX9eA zmH5M?@Vo6J)jut*&orkWl(A_R6l(3RUnOxma^l=*vNiUNN3w~QJ`#5fE@FRhI-A73qi5dD$ zs}-G=rIGh6k1Gk9MsA+HDyb9JCWW&!DIK-<6qJ(ul8J_CK5x8iI~%-n!8A6Cb-sJ2 zElpA>nq;_{-m%zuo%$^@B9csIDeH8NhDuXUZ@C6u+LL*VI7wP$zr)v&^-req-kyqg zj+@MRZ#XN*0iapYr!6O zNtH6>jGw@V7LiWR6=Y zf-X#FSs1QlvnNnxP~FBdBBPU}jAUt-+q~Roy~XOJ&)ehUnkM9lwIhOV>|od}7S4@G z!CvJ|h}^=s-N&@Qvpwx^8eHd`5@`>AtKG(HaLY;BQI|@`SsX1DJ4@!rP@IX~y@Rl~ zE>cvUW-oNH4^H-V-)4r9H0 zXN2(GG}9e5(?f8N>3q&FQYxfc+Z!Wh!(?nWM#ExkO@m|%0l$Y)1}n0UIsj@wLZT!> zMSwwwq#y#qA(IIZ%uTcBo~d;D8PTha83XK literal 0 HcmV?d00001 diff --git a/R/ui.R b/R/ui.R index 9040149..b7e071f 100644 --- a/R/ui.R +++ b/R/ui.R @@ -141,10 +141,33 @@ ui <- function() { ) ), tabPanel( - title = "Help" - ), - tabPanel( - title = "Publication" + title = "Additional information", + div( + class = "container", + h2("Number of interesting genes along the ranking"), + p(HTML(paste0( + "The notion of ubiquitous genes or housekeeping genes implies ", + "some kind of enrichment of important genes. Within groups of ", + "genes with certain biological associations, these genes should ", + "be overrepresented. We use GeneOntology terms as well as some ", + "other gene set sources to represent this concept. The following ", + "plot shows the number of associated terms for a non-overlapping ", + "sliding window of 500 genes along the ranking of ubiquity using ", + "our default parameters. The terms have been obtained using a ", + "gene set enrichment analysis with the tool ", + "g:Profiler. We observe that the most ", + "ubiquitous genes have many more known biological implications ", + "than any other bucket of genes. The genes of average ubiquity ", + "have almost no associations with GeneOntology terms. The number ", + "of associations rises again for the least ubiquitous genes." + ))), + p(HTML(paste0( + "Note: Click on the legend items to toggle single sources. A ", + "double-click will isolate a single source of interest." + ))), + plotly::plotlyOutput("gsea_plot_ranking", height = "600px") + ) ) ) ) diff --git a/scripts/gsea.R b/scripts/gsea.R index e0d0edb..6fb92ca 100644 --- a/scripts/gsea.R +++ b/scripts/gsea.R @@ -45,3 +45,6 @@ fig <- plotly::plot_ly() |> ) plotly::save_image(fig, image_path, width = 1200, height = 800) + +gsea_plot_ranking <- fig +usethis::use_data(gsea_plot_ranking, internal = TRUE, overwrite = TRUE)