From 18505989a5903c8707fab5291eb146e093e7ed97 Mon Sep 17 00:00:00 2001 From: shashank-KU Date: Mon, 6 Mar 2023 11:35:32 +0100 Subject: [PATCH] Initial commit --- .Rbuildignore | 5 + .gitignore | 6 + CODE_OF_CONDUCT.md | 126 ++++++++++++++++++ DESCRIPTION | 29 ++++ LICENSE | 2 + LICENSE.md | 21 +++ NAMESPACE | 3 + OmicsIntegrator.Rproj | 22 +++ R/OmicsIntegrator-package.R | 6 + R/calculate_correlations.R | 67 ++++++++++ README.md | 108 +++++++++++++++ data-raw/DATASET.R | 5 + data-raw/wgcna_metagenomics_eigengenes.rds | Bin 0 -> 1857 bytes ...wgcna_transcriptomics_Liver_eigengenes.rds | Bin 0 -> 21655 bytes data/Metagenomics.rda | Bin 0 -> 2191 bytes data/Transcriptomics.rda | Bin 0 -> 22187 bytes man/OmicsIntegrator-package.Rd | 15 +++ man/calculate_correlations.Rd | 39 ++++++ tests/testthat.R | 12 ++ tests/testthat/test-OmicsIntegrator.R | 3 + vignettes/.gitignore | 2 + vignettes/OmicsIntegrator.Rmd | 19 +++ 22 files changed, 490 insertions(+) create mode 100644 .Rbuildignore create mode 100644 .gitignore create mode 100644 CODE_OF_CONDUCT.md create mode 100644 DESCRIPTION create mode 100644 LICENSE create mode 100644 LICENSE.md create mode 100644 NAMESPACE create mode 100644 OmicsIntegrator.Rproj create mode 100644 R/OmicsIntegrator-package.R create mode 100644 R/calculate_correlations.R create mode 100644 README.md create mode 100644 data-raw/DATASET.R create mode 100644 data-raw/wgcna_metagenomics_eigengenes.rds create mode 100644 data-raw/wgcna_transcriptomics_Liver_eigengenes.rds create mode 100644 data/Metagenomics.rda create mode 100644 data/Transcriptomics.rda create mode 100644 man/OmicsIntegrator-package.Rd create mode 100644 man/calculate_correlations.Rd create mode 100644 tests/testthat.R create mode 100644 tests/testthat/test-OmicsIntegrator.R create mode 100644 vignettes/.gitignore create mode 100644 vignettes/OmicsIntegrator.Rmd diff --git a/.Rbuildignore b/.Rbuildignore new file mode 100644 index 0000000..8297bf0 --- /dev/null +++ b/.Rbuildignore @@ -0,0 +1,5 @@ +^OmicsIntegrator\.Rproj$ +^\.Rproj\.user$ +^data-raw$ +^LICENSE\.md$ +^CODE_OF_CONDUCT\.md$ diff --git a/.gitignore b/.gitignore new file mode 100644 index 0000000..0542f56 --- /dev/null +++ b/.gitignore @@ -0,0 +1,6 @@ +.Rproj.user +inst/doc +.Rhistory +.Rdata +.httr-oauth +.DS_Store diff --git a/CODE_OF_CONDUCT.md b/CODE_OF_CONDUCT.md new file mode 100644 index 0000000..63a79b6 --- /dev/null +++ b/CODE_OF_CONDUCT.md @@ -0,0 +1,126 @@ +# Contributor Covenant Code of Conduct + +## Our Pledge + +We as members, contributors, and leaders pledge to make participation in our +community a harassment-free experience for everyone, regardless of age, body +size, visible or invisible disability, ethnicity, sex characteristics, gender +identity and expression, level of experience, education, socio-economic status, +nationality, personal appearance, race, caste, color, religion, or sexual +identity and orientation. + +We pledge to act and interact in ways that contribute to an open, welcoming, +diverse, inclusive, and healthy community. + +## Our Standards + +Examples of behavior that contributes to a positive environment for our +community include: + +* Demonstrating empathy and kindness toward other people +* Being respectful of differing opinions, viewpoints, and experiences +* Giving and gracefully accepting constructive feedback +* Accepting responsibility and apologizing to those affected by our mistakes, + and learning from the experience +* Focusing on what is best not just for us as individuals, but for the overall + community + +Examples of unacceptable behavior include: + +* The use of sexualized language or imagery, and sexual attention or advances of + any kind +* Trolling, insulting or derogatory comments, and personal or political attacks +* Public or private harassment +* Publishing others' private information, such as a physical or email address, + without their explicit permission +* Other conduct which could reasonably be considered inappropriate in a + professional setting + +## Enforcement Responsibilities + +Community leaders are responsible for clarifying and enforcing our standards of +acceptable behavior and will take appropriate and fair corrective action in +response to any behavior that they deem inappropriate, threatening, offensive, +or harmful. + +Community leaders have the right and responsibility to remove, edit, or reject +comments, commits, code, wiki edits, issues, and other contributions that are +not aligned to this Code of Conduct, and will communicate reasons for moderation +decisions when appropriate. + +## Scope + +This Code of Conduct applies within all community spaces, and also applies when +an individual is officially representing the community in public spaces. +Examples of representing our community include using an official e-mail address, +posting via an official social media account, or acting as an appointed +representative at an online or offline event. + +## Enforcement + +Instances of abusive, harassing, or otherwise unacceptable behavior may be +reported to the community leaders responsible for enforcement at bioinformaticssrm2011@gmail.com. +All complaints will be reviewed and investigated promptly and fairly. + +All community leaders are obligated to respect the privacy and security of the +reporter of any incident. + +## Enforcement Guidelines + +Community leaders will follow these Community Impact Guidelines in determining +the consequences for any action they deem in violation of this Code of Conduct: + +### 1. Correction + +**Community Impact**: Use of inappropriate language or other behavior deemed +unprofessional or unwelcome in the community. + +**Consequence**: A private, written warning from community leaders, providing +clarity around the nature of the violation and an explanation of why the +behavior was inappropriate. A public apology may be requested. + +### 2. Warning + +**Community Impact**: A violation through a single incident or series of +actions. + +**Consequence**: A warning with consequences for continued behavior. No +interaction with the people involved, including unsolicited interaction with +those enforcing the Code of Conduct, for a specified period of time. This +includes avoiding interactions in community spaces as well as external channels +like social media. Violating these terms may lead to a temporary or permanent +ban. + +### 3. Temporary Ban + +**Community Impact**: A serious violation of community standards, including +sustained inappropriate behavior. + +**Consequence**: A temporary ban from any sort of interaction or public +communication with the community for a specified period of time. No public or +private interaction with the people involved, including unsolicited interaction +with those enforcing the Code of Conduct, is allowed during this period. +Violating these terms may lead to a permanent ban. + +### 4. Permanent Ban + +**Community Impact**: Demonstrating a pattern of violation of community +standards, including sustained inappropriate behavior, harassment of an +individual, or aggression toward or disparagement of classes of individuals. + +**Consequence**: A permanent ban from any sort of public interaction within the +community. + +## Attribution + +This Code of Conduct is adapted from the [Contributor Covenant][homepage], +version 2.1, available at +. + +Community Impact Guidelines were inspired by +[Mozilla's code of conduct enforcement ladder][https://github.com/mozilla/inclusion]. + +For answers to common questions about this code of conduct, see the FAQ at +. Translations are available at . + +[homepage]: https://www.contributor-covenant.org diff --git a/DESCRIPTION b/DESCRIPTION new file mode 100644 index 0000000..9ce9ee3 --- /dev/null +++ b/DESCRIPTION @@ -0,0 +1,29 @@ +Package: OmicsIntegrator +Title: Calculate correlation between two data frame +Version: 0.1 +Authors@R: + person("Shashank", "Gupta", , "bioinformaticssrm2011@gmail.com", role = c("aut", "cre"), + comment = c(ORCID = "0000-0001-8839-7003")) +Description: This function takes in two data frames df1 and df2, and calculates the correlation matrix between the columns of both data frames. The user can specify the correlation method to be used by passing the method parameter. By default, the function uses Pearson correlation. + The function calculates the p-values for the correlation tests and adjusts them for multiple testing using the specified adjustment method (default method is False Discovery Rate adjustment). The function returns a list of the correlation results including the correlation matrix, the p-value matrix, the adjusted p-value matrix, and the significance matrix indicating the significance level of the correlations. +License: MIT + file LICENSE +Encoding: UTF-8 +Roxygen: list(markdown = TRUE) +RoxygenNote: 7.2.3 +Suggests: + knitr, + rmarkdown, + testthat (>= 3.0.0) +VignetteBuilder: knitr +Imports: + cowplot, + dplyr, + ggplot2, + magrittr, + pheatmap, + RColorBrewer, + WGCNA +Depends: + R (>= 2.10) +LazyData: true +Config/testthat/edition: 3 diff --git a/LICENSE b/LICENSE new file mode 100644 index 0000000..92464a7 --- /dev/null +++ b/LICENSE @@ -0,0 +1,2 @@ +YEAR: 2023 +COPYRIGHT HOLDER: OmicsIntegrator authors diff --git a/LICENSE.md b/LICENSE.md new file mode 100644 index 0000000..db3e795 --- /dev/null +++ b/LICENSE.md @@ -0,0 +1,21 @@ +# MIT License + +Copyright (c) 2023 OmicsIntegrator authors + +Permission is hereby granted, free of charge, to any person obtaining a copy +of this software and associated documentation files (the "Software"), to deal +in the Software without restriction, including without limitation the rights +to use, copy, modify, merge, publish, distribute, sublicense, and/or sell +copies of the Software, and to permit persons to whom the Software is +furnished to do so, subject to the following conditions: + +The above copyright notice and this permission notice shall be included in all +copies or substantial portions of the Software. + +THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR +IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, +FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE +AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER +LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, +OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE +SOFTWARE. diff --git a/NAMESPACE b/NAMESPACE new file mode 100644 index 0000000..c61ced3 --- /dev/null +++ b/NAMESPACE @@ -0,0 +1,3 @@ +# Generated by roxygen2: do not edit by hand + +export(calculate_correlations) diff --git a/OmicsIntegrator.Rproj b/OmicsIntegrator.Rproj new file mode 100644 index 0000000..69fafd4 --- /dev/null +++ b/OmicsIntegrator.Rproj @@ -0,0 +1,22 @@ +Version: 1.0 + +RestoreWorkspace: No +SaveWorkspace: No +AlwaysSaveHistory: Default + +EnableCodeIndexing: Yes +UseSpacesForTab: Yes +NumSpacesForTab: 2 +Encoding: UTF-8 + +RnwWeave: Sweave +LaTeX: pdfLaTeX + +AutoAppendNewline: Yes +StripTrailingWhitespace: Yes +LineEndingConversion: Posix + +BuildType: Package +PackageUseDevtools: Yes +PackageInstallArgs: --no-multiarch --with-keep.source +PackageRoxygenize: rd,collate,namespace diff --git a/R/OmicsIntegrator-package.R b/R/OmicsIntegrator-package.R new file mode 100644 index 0000000..a65cf64 --- /dev/null +++ b/R/OmicsIntegrator-package.R @@ -0,0 +1,6 @@ +#' @keywords internal +"_PACKAGE" + +## usethis namespace: start +## usethis namespace: end +NULL diff --git a/R/calculate_correlations.R b/R/calculate_correlations.R new file mode 100644 index 0000000..c8c3830 --- /dev/null +++ b/R/calculate_correlations.R @@ -0,0 +1,67 @@ +#' Calculate correlation between two data frames +#' +#' @description +#' This function calculates the correlation matrix between the columns of two data frames. +#' +#' @param df1 A data frame containing the first set of variables to be correlated. The samples should be in the rows and the MEs in the columns. +#' @param df2 A data frame containing the second set of variables to be correlated. The samples should be in the rows and the MEs in the columns. +#' @param method (optional) The correlation method to be used. Default is "pearson". +#' @param adjust_method (optional) The method for adjusting the p-values for multiple testing. Default is "fdr". +#' +#' @return A list containing the correlation results +#' @return$correlation Matrix of correlation coefficients. +#' @return$p_value Matrix of p-values for the correlation coefficients. +#' @return$p_value_adj Matrix of adjusted p-values for the correlation coefficients. +#' @return$signif_matrix Matrix indicating the significance level of the correlations. +#' +#' @examples +#' # Calculate correlation between two data frames using default parameters +#' result_list <- calculate_correlations(df1 = transcriptomics_data, df2 = metagenomics_data) +#' +#' # Calculate correlation between two data frames using Kendall correlation and Bonferroni adjustment +#' result_list <- calculate_correlations(df1, df2, method = "kendall", adjust_method = "bonferroni") +#' +#' @export +calculate_correlations <- function(df1, df2, method="pearson", adjust_method="fdr") { + # Calculate correlation matrix + cor_mat <- cor(df1, df2, method=method) + + # Calculate p-values and adjust for multiple testing + p_val_mat <- matrix(ncol=ncol(df2), nrow=ncol(df1), + dimnames=list(colnames(df1), colnames(df2))) + corr_val_mat <- p_val_mat + + for (i in 1:ncol(df1)) { + for (j in 1:ncol(df2)) { + cor_res <- cor.test(df1[,i], df2[,j], method=method) + p_val_mat[i,j] <- cor_res$p.value + corr_val_mat[i,j] <- cor_res$estimate + } + } + + p_val_adj_mat <- p.adjust(p_val_mat, method=adjust_method) + + # Add significance level. + # One star means a p-value of less than 0.05; + # Two stars is less than 0.01, and three, is less than 0.001. + signif_matrix <- rep("", length(p_val_mat)) + three_star <- which(p_val_mat <= 0.001) + signif_matrix[three_star] <- "***" + two_star <- which((p_val_mat <= 0.01) & (p_val_mat > 0.001)) + signif_matrix[two_star] <- "**" + one_star <- which((p_val_mat <= 0.05) & (p_val_mat > 0.01)) + signif_matrix[one_star] <- "*" + dim(signif_matrix) <- dim(p_val_mat) # Give textMatrix the correct dimensions + + # Collect all results into a list + results <- list(p_value=p_val_mat, + p_value_adj=p_val_adj_mat, + signif_matrix=signif_matrix, + correlation=corr_val_mat) + + return(results) +} + + + + diff --git a/README.md b/README.md new file mode 100644 index 0000000..e099e50 --- /dev/null +++ b/README.md @@ -0,0 +1,108 @@ + +# OmicsIntegrator + + + + +The code provided in this repository can be used as a framework for exploring correlations between different omics datasets, including transcriptomics, metatranscriptomics, metagenomics,metaproteomics, and others. By utilizing packages such as WGCNA, pheatmap, and cowplot, it provides a flexible and customizable way to generate integrated heatmaps and visualize correlations between datasets. + +With some modifications to the input datasets, this code can be adapted to work with different omics data types. For example, in Step 3, users can substitute in their own data frames to calculate correlations between different datasets. Additionally, the color schemes used in the heatmaps can be adjusted to best suit the user's data. + +Overall, this code provides a powerful tool for exploring correlations between different omics datasets, and can be adapted to work with a wide range of data types. + +## Installation + +Install R (version >= 3.6.0) and RStudio (optional) +Install the devtools package in R using install.packages("devtools") + +You can install the development version of OmicsIntegrator like so: + +``` r +devtools::install_github("shashank-KU/OmicsIntegrator") +``` +## Steps + +The following steps use OmicsIntegrator to integrate transcriptomics and metagenomics data and visualize the results: + +### Step 1: Perform hierarchical clustering of transcriptomics data using WGCNA + +Load the WGCNA package in R +Calculate the bicorrelation matrix of transcriptomics data with outlier removal using the bicor() function +Convert the matrix to a distance matrix using as.dist() +Perform hierarchical clustering on the distance matrix using hclust() with "ward.D2" method + +``` r +library(OmicsIntegrator) +Metagenomics <- readRDS("data-raw/wgcna_metagenomics_eigengenes.rds") +Transcriptomics <- readRDS("data-raw/wgcna_transcriptomics_Liver_eigengenes.rds") + +dendro <- hclust(as.dist(1 - WGCNA::bicor(Transcriptomics, maxPOutliers = 0.05)), method = "ward.D2") +``` + +### Step 2: Generate a heatmap of transcriptomics data with dendrogram + +Transpose the transcriptomics data using t() +Create a data frame from the transposed data using data.frame() +Generate a heatmap of the data using pheatmap() with the dendrogram from Step 1 and no column tree + +``` r +Transcriptomics<-data.frame(t(Transcriptomics)) +result2 <-pheatmap(Transcriptomics, + cluster_rows = dendro) +``` + +### Step 3: Calculate correlations between transcriptomics and metagenomics data using OmicsIntegrator + +Use the calculate_correlations() function from OmicsIntegrator to calculate Pearson correlations between the transposed transcriptomics data and the metagenomics data + +``` r +result3 <- calculate_correlations(df1 = t(Transcriptomics), + df2 = Metagenomics) +``` + +### Step 4: Generate a heatmap of the correlations between transcriptomics and metagenomics data + +Create a color ramp for the heatmap using colorRampPalette() +Generate a heatmap of the correlations using pheatmap() with the dendrogram from Step 1, the color ramp, and no row tree +Add significant correlations to the heatmap using the display_numbers parameter from pheatmap() + +``` r +heatmap_colors <- colorRampPalette(rev(RColorBrewer::brewer.pal(n = 6, name ="RdBu")))(51) + +result4 <- pheatmap::pheatmap(result2$correlation, + color = heatmap_colors, + treeheight_col = 0, + treeheight_row = 0, + cluster_rows = dendro, + #cutree_rows = row_cut, + display_numbers = result2$signif_matrix, + breaks = seq(from = -1, to = 1, length.out = 51), + show_rownames = T, legend = T, + labels_row = paste0(rownames(result2$correlation)), + labels_col = paste0(colnames(result2$correlation))) +``` + +Step 5: Combine the heatmap of transcriptomics data and the heatmap of correlations + +Use the plot_grid() function from the cowplot package to combine the two heatmaps into a single figure with two columns +Set the relative widths of the two columns using the rel_widths parameter +Adjust the margins of the plot using ggplot2::theme() with the plot.margin parameter + + +``` r +cowplot::plot_grid(result2$gtable, + result4$gtable, + ncol = 2, + rel_widths = c(3.5, 1)) + + ggplot2::theme(plot.margin = ggplot2::unit(c(1,1,1,1), "cm")) +``` + + + + + + + + + + diff --git a/data-raw/DATASET.R b/data-raw/DATASET.R new file mode 100644 index 0000000..03a4333 --- /dev/null +++ b/data-raw/DATASET.R @@ -0,0 +1,5 @@ +## code to prepare `DATASET` dataset goes here +Metagenomics <- readRDS("data-raw/wgcna_metagenomics_eigengenes.rds") +Transcriptomics <- readRDS("data-raw/wgcna_transcriptomics_Liver_eigengenes.rds") +usethis::use_data(Transcriptomics, Metagenomics, overwrite = TRUE) +usethis::use_data(Metagenomics, overwrite = TRUE) diff --git a/data-raw/wgcna_metagenomics_eigengenes.rds b/data-raw/wgcna_metagenomics_eigengenes.rds new file mode 100644 index 0000000000000000000000000000000000000000..bd8350a86c53a5de3d07a1c62268734b836debfa GIT binary patch literal 1857 zcmV-H2fp|piwFP!000001BI4(R8w~r$3u`R$Ve1W5G+`*t0KIo2#Wkt6$%Jpl-0Ba z5k)BhBm@i)6a&Ek2}=N3iXb$CVNqmX#Cw>j1+mID$`T9Au!v}BN2`L4%w%c^oxJ0j zKi>J?`@Q%6e(&=;@167W)*z8cs-$IPGD%fU7LwFF-0inm$vkaYvRwA8l0BxtqkN-l z_Y)w#bTF=b&kD#^D^4{(iuxwy{I!N{_I}PccbIH&5IVt+nsy#Y?4!2x9C&XkMKZ}|$hxCyR z9zUkNfwV^&milRp5Wnta{fdXGAh3zQQ(^UINLeVQH@sf~DQrXP%=9p%qz0yo^Ll_0 z{t?u!@Bbg_Mq~Z!C(w3`nWe{c$40laGKt4 zQD!@+v~L+ix^2onVc8EckBdUjel!fMOGV#byQP8n?wJeoTt12}FBG4MD^>I-l*}<~ zMvp4asgfm0;6@jsI)QaQr@$CheL31(rrC$ux(sy&+8!YA=#F%uu0*YBpEVful%twm z{Ef!`bX4!%B6xa1AGP?r_O6=Vh8kvscXbckL(P#6rC6m9)&Af&lyt`p%2}TY)o+f$ z@#;!JoM#hUqc2FiT!W!%MmxVO?GEIfc$~4?;S`9r`)K~Kdj?4&e+;y&PDY{tiQ|0_ z4k{U8bde^kQL0ZumLrV=IcFx0;FAvM%FaD4ka!tdF7c(ocYUGGmu7GF$t`GPkL=)j z-G#ct9SNRW(x7%eN7(*Y3u^w>_xj)5wNP_2Y^cIT2n|J@3BBj@p{C>=qW4>vbJQ=kr(n=0vVW zX*|ZU&#w=myVVlkUHLXBW{OWU%*+H~-Sd&sXLN`u8Jif_5~Ik}&10faE(*!{xY5qT z95~XAPIqEeAcNPZ(!MNTu31WXJ_sH^XAKm|zN5T8A6K_QMPG`;+zK5NO>}N_CQ|Fo(m5puZmYXq`UIkO2Q>`Uvmt$_ zbL#ylM}_ic-2-dy)uAur4C}WJd=E)m&yZcMP7(T&JIU6$CSCH}0eSk_lS^I*6!nw> z3Z+W-1^6*y@7urE2gyCf3q#uW@;X(n`{in+W^P?Ml9sN&%cSz zEQJW&F(*6n5aOKF$)+C;RrH_D_^_e%>q2yyA)(|N^#k9lOJnnV7kCy5iri9a73+z! zV>v?A$i-Z7P8mc0d{t9XPEg^n`7c=+$L0AILfvQ9rs+^zvyxvb*Q$4{dmfT(P2+R4 z(%(@pl&k0+b?RHK)LQgWUTs7~aBmdKC_7v>N52n=&(aSxo(qDQj0MisUNJ;}wftv2 zi8%@@zWAa!Vh8Z7_djB{Xrhyx@_PPGCaQKz!N*n}Ma6HN2c~nzQSY8Yo73DW6!!Fe za;z{^aZVx!-``AjUHCux9q?U;i|17UVo1}MPpor+L~Hi!hcs*88ECe)&UDKAz6HTG z-J&-y`xWyk8wUe>8s0;WJ17t7o`Zp?}c5ddYSuiu(Efe&0{|ecpd~IQQ{d&*ybs&(}HEb-6X^C@3hXD5xnZ zDX3`HO$wT8>gP^~ultzSKPavLFs}b_gZ;O`1=rPwSaa!Fy!ZBpSiah_DPSZB-MA%} zJ}#DFf#|N;;NyQVcX@kF(!p(5J$B&Q!pv1nV|_eDeo=r`LIRvDO14PSJZD0gV1`XL zBO-JR3z*R;XXj|O32I)H<^FC+f%4d;vipd#hErb--0tm_F#2DCY1)ByabNh8T921@p^Y zNIYJ6aEkL9ILq93+M#KSrEl&UzuGy4ew<=AxT$%dyw0NXB=t@V^4+>}SU(8U$240O zRL`LC`vk>t_yxXfiBijFZ7|~2iB@A-4XF8j+R1=h5K@O32W*4up`KNii`HxzNWKFS zbpk;c6?R{;la3Y=#;Ts3e3J)-7VP7XG^fF4i%ehYrmK*c>rZ>L$p{i} zGyFQN_8y6%PUW?kdoW>T>S0*EKN7#36&s|9K(d&q`tt!lTT)>-?L$g#6B9@-O*Te z^G5&ai~LwU_dANe5U}>k9?|lS_gFb>x~#Ng6>DR|9}4Dd!loU?=hNlokyxhZJ+kX9 zR##YG`dxAklhV1nls_1srFtv$5cUbW4? zTLjGu`8K7K!%&v768rM~DQFfuxu0U+EYx&XB&mJKg9eVIy&(gRP*yz}p5E#O#7$Gr z1sb`b;oQ;Y;Ih}yBBOPdN}&Rp?Wak!2Bg0{gU7*pITq6$ETrnux>l} z{>)#P=(oa(9U2FNNA_U4ucL&(KUQddw>#Nq zjt#{Z`blr6QYwm$RJLcCPx#}FMn;|E z&ndAguzuu*{ck8_Wf^|T9EM@lpN@PY4?<42xMJ*UZ^-QkJ3Pr*hh^+_-`jbYki2Nc zg=e^cL~;Ck>5FD)xxV80s^=0^ykNT8{)Y>)Cn#y6LgJxqbg0fTkO^|OzkcU$*9oP{ z4J?s;;uE8%P|+pyxXbhckj1Q2C?AOfDQExW%`0~Q!xIt`j@upe z2*!r1HoYy^*s#v!7tM;CC6dp@T%fL#MdCHdS|+7QY^M66Y1Z)$D};G{Hr z^<`D8OU<)qv#&=2*ZiXLOJ%Grd&@}p<%M-S=3h`N)kBnNhTCP;PH?4TeOSJ)9^?6$ zuPpQTV@?;IkQ}pzq{~mFc2a*vgNq-!N`IaQUpZwv+1@8u@_GNuYqXV6e0`vp&EpY7 zi!WqPE?vZ|;|U!Gj+U5x;FR!0JPS0oryJfKzX{FFv!CYG%%D*)r_WK58?xu0p3OO^ z2t;|p+V)3hfJ|Klw_;?Vj@hX6BL6HDYRoAAW)%QZ#pB$>hjCCRdSy=4t_tcWStq_f z@qrc(L-)s~^?&2Zo6@#xFaLM!5HD_IY-YtYMTL6xt=?F*Q(5`mR|;$&dfC70ha5(9 zC^DUg zQ2MgfA*F82_+-H+Zz_uUmd)8m+Y}&*`BPqe!V^fDad6j`(89dO^@pckT!th|t+L4q z9!!3Gz+Hqx6LN}09JQPip~@Xdj5gn}`W9* zqcl|a2%LESZ4D}IM@Js`O+uZ?<1~hROCVO(ggty;3q<>;Tn|Kyp{ir|wnCR_C|u3D zyU$Dz8f;{Gf^6-8+;rti-gz+# zI!d)dL#LR{-usMDd-gT!L^>PPBtdB$R|@1klPBnXdkHl4AgK3 zJP&x_06u|KyGZnOs21hpr#s(+*I(XtlKRAkX(ukaeAX7g)PBS0W(s#m_(`?4^KBY< zzjw>Z3n63jzRZg4F1w&$+vhVqU8)ej>v4y)OeSQF^o)SlBk+&T6~ZHde|aJ~UsHa( zzX`*aN5{9aoyQW@uS@UpT(N9r`@J2&{P(_I(iQ0|(VbIJZGFdC=ZeY9};K zrkokIPld*GhF#t}Wub*Tak%{+9Uz}u=56h2z#BsCqTgDe;UKMxx=Sn&4L|h8zFvj8 zmLprWE+|1W|CZDS!)H*pTm1EwrkhyY+h*s(LV@K<+=q=1Eg)G}_F9gyDb^^x>tZgL z!)kY-QJOb}SRo#o;zB!xmARK>vpYkvp{ynB*$=Vx_;fpVV!Q`S@zJg_z5AFODO;4! z$AILE%?Afl|MApJdFB_rS}o*ptn$|^TVUDzmypK2K3J1vDAaFyA3`V>p5*H(W0TXP zxR~ZqjLARn^Zq&>n{Hp9wY~(*rqif)w)=nt#k|nYre`RBqU5dJiABsh@^!x()lVRQ zI&dI$Mh{vWwNgw(B%qO4m$(v9klnlov<`cZJbr5p)w$cmo{o4y z!|GGQO2`B_%*wD0j68*e5(a}$yOW?zA-V9zsXb5=aGUY@wOJ%i-KBU!@eWH2UcpSL1sGR7NRd|E7bze5M%2%Repk2$| z%b&2|VewVQtrsx9f$q`Z(Ksxq0%wPgn}2z#>{iorT^NJlJi~YK4)*A^ylTvT`vHc? znEm;D;1+}$cBjX71!Inl%Iw(o4)E|H@KSo8#)A2{^EP7h5ZdDuNXu>jz700#*D_5Y z`3?0{*SlX(Le_utO{yIdXS>Enbcs-7#B=RO*C1qc=k3ybegYbH)9pCi)(f#&GWX2u z*no7{^oE!OEmSf6Sz@D30fMdaL)tgG;Cc$^20mK=G3K$hLe3v3mfpvHrqUb=9mBh> zA6Y>EFz>kmkC#ZGG&=mM+!zV3S=KZjp2DV#+kDLwFR}D<(4^$R5loeAO+5Nt6Ma4% z@KJtp5!+77EI&CWh}m5=J8Md7FCzee6kn@e?h(%&xYEN=(YVb>A-H#oaE#X%+bSq6^hO8 zW;wxj&5G~p_H>LLQh%IetO=nNZCBI-UBEGx6V->yA-7#z=|=QT~Z{Z{*sKc|5toS?&bAPP$Un04Lzwh1biR&VgU;K!8b&CMmo0+6Y6 z^Y-3cHz4LPF>n7a1;w3FR>?+t(V3NZB*u&ev)#tld>L~wtMcSkC#4arVN)KbjTgd# z1ga;lJg31i@`^u6^c?zII-)`2Ja*maaj|Vphe(q{ocHYCV5miOa3UX~jlQ9Y*HAoW znzx4QYaID|-es`{zN=lm4W%?`gVdFKfbcsj18P~J>Y;Uv~7MC%9kKvhjAg3niQ0XM?JDRcnyfIK2>AKZ$N3&?O8b84|PX+xE@-|fI~!R zO3|0!5YOyL2-Glyh!BI1$6p_S?B!z#%suXym%BJk;5Y;+(agb@=c6H-U-TQZ-+8bz zBNY99lMTTf&Vu7Rh>%eq6m&u(0fHV;n(P}Hft;A{^3hKcu`od5yY%BNSkl(>?l`w9 z=EwTpaeR9pJB!8Is|OsgCZ?x8mdgW+PFon7^xnrrfe~GYpAWH%?f1S#Iu*=aQ2u$j zJqZh%GXtz83Nd*LmHCxZ+p$bZxzsy_<8M4g@Oc^c6ZN3}@Y-*eg?H@CBhWC3yW8+9r*6(}8U`6Y42 z4#}LeOHvd%*uc@))KHXzbq`Jj6Bit@#ZLbI2M@9Jy!S@+ci;6qFW_`m=SvZ+S;_mj z#nlbFvLbo|d9<+Ru`|WYNdc_SH9XLn@?kxHofejSNx-JEo3Fe%694j4WE!da2d1Is zoEgI&85gi!WoOUMSI6+rhj{566o9leHg)h4Gv@n{O>#3+ph!cDZf`pa6#u$hKRU7% z682C}oV}X_NzKKZ&+F#`@y_1Jc%gR4`pNwB`xRTrvo3E9TdjbqXJ5t#AD#!&!GyTa z+uI@U5Q!!|gBNo5G=|w%bV0RK+<4@*JCGg!L1>$q{d#^_-1Fq)FO1M9eEHNug&|5@ zC-$eSLXGaxxOCljQ2sWodjI)gBu(Ba+@WTSRR-rpybPnUw0&H>Ub6|?-aMXuvHt{C zEt?yx_9$b)O~vO1cX^P=bYG6|X+CzHE~^!hD#N%>!3S$T@M0N>Kbb|c3k%k?x68#f zWBKC4aaYB^bugA3uGpIJ)((n=R9;-H-Ui{hOiDlAy@Rk@TkZ}JxS?c;ynp8DFLd-U zUmT!bhJuq`g`@Us5L!6k>s|B~;$G|)KUk##ImaG)C;xZ}#E)L`+Xly=M8UXXfvydr z?NXCoQWQgVoF^}=xk2Xqjr#&C(omjxW-Mi%9rW1kg1#_@L)BBeXyNxfQ1FNJ@z8Dy zAk2oyF(~L^1U<{ym}jbxoHtS?JlzYWkGIdb24_N1+MR6+x%rs$HcH8IZ5&C`gF>Ek zu~>P1*7lO{18nV5_g1;I6Dy82iLP8&zs~!2xV%+3hcSU|XB6ALv8#+2;!Eix zF_4qz$#=;72^6bjvHQQd1C@gchsRr0fi&dc&u;Mzs=ns#k@n?=vP&6(((&O?RQLSD znQRWIuMi-V8M{D@{ao!?xh1HS`lCBkGzXP(Dew1p>OgfPPgb*r8dRPPyi9wF2(`RC zl*W!pSoWJqD5pXQTi|k;z1s(@z3wY0e%}a5cQ)g_`c`b&Y@e;dl!~pQ0&Wc-&SFh} z#frOw3AWBk7CFwnL*hq+;vK2SvCb%so#sy>HmY$ku-eUFYktJH9#et8JQeuzoe$jX z1JTyO@`GJ!m@SZcPo&op(rfz-i+oHV;`e|t&Bb|0+RETco16)uj{MIiIJThgZLNYQ z73p9ME$+C|oz!`hMq@)PZ95cc&Zhgi=U zteX9~xX8Nc0)ylsPy`RCx%bhB# z-7{D}QFB;&>y^JeWs7)Cbd7F@g4;}O;|p|{DHk%u615k7&x<8ck+(urWTD-n;U=^# zx={S=votsJ^YP>ubPucV|}Fhp4}~Ke{nH8B)8&DnE(lA&KE)oKuiJ zR8S9RB);~6l3dBROX300*dV%4W^DmE?5-?N{iuP|{8ZA!@dPyKPJqlx&U$_}7hbtv z7m5ka6WeddK!xXhm%dCP$TQh~v%>QpG_=$tWW7Pq2x4PP0tS+eANjZH7JscA+ zDt)8F9P^OW@=I=5cW-rN&?p!Q79-z+_No5mNspAHPrvgK)DG=9p{x`EIXTbIB--{s zoVA;|^!`O~dtoC{FmF~kQBvOo>xM1S=H-B=D9O1l6 zh=mG@Bk5!B*&#W)Om+NfEs*}~kF*rM0d=E^oi2^RKo-vnX%yWLm3k>h&z&uYf-vC# z`g(){;}5c(cH7q5+ew(^1w%F0RA5+AI27;Up{qFR44yexTuWu5GEakX~%Kk&w@-ab43S59;0$ z(}D^TAnWAxu^clCNN$uqDPvg(`Gb5b3Z4Z}1l`xkdP7jyWv_nci2$^Ujj0b^=zt1* zDf@Cd8)|kcbzY*22J(!&$@nfs$USoWnrPA%DA1HY^1P!Q9WAv_8GNh7(%L=LWAklT z_hV@LRzp*~)`^!_2<+r>eAC6&M-G1J$iipt$bvEyNtKe3Z z{DlK^is08Wkldl94sN^mO-{sHq2KGF;?f<#m{ygnwAi7Dc^nU!9Qc+YgTBE%`^_o_ zTI}|4#aS4LY2X?sL+6TpE_KSJH4InP^gLLg6 z1*Gg!PoP`K2U70w?*07TkYF5YsYg)&bqb?u-nLUvkr()8po16UE;#qUW%&(-PhVbI zig^ZwizX;xmx;m0A2uevs)t}hvoqXp*)VfJFQ?#iAGY>Z86WZ-!`4X0L!}}TSgCgL zXzG&(SaEfL;tQb@!|dOh$gc0pe4}3ut#>^!ay3!-qVzZts?#j`3EqEs%0ARpQ)cxA zt)J-kJJwI&sj|;nDoxpt^TlA1_4yNs?pQrg`*Ru!qW7kkCjrFTSdR@fjez4hn$pnc zBbYOk$C?u=0+j;t_ZTVYAvgb^dG(Mc1fB_FZL^((n)}<&@o80J9-(AK`~kgDN=O;y zwY$w%0zs-P_r8-opafhbg5x|F8#+I-Ocqh&$mRN}^J0PwMU_|es6K9nT#xBZMUm~0yZA;K zS-l|p<65s$bq#nXXM9tzQN#p#^&Ojbmq87$5T5|QI#h~Zu9INrTCZ~hb4I8Jpy0E% zt>nXOOeKEU;lfmhVNF_$Dv{TqKy9c)zEBa1s|2AdWsza+*(*zLCu_cm<3D3{O9k~G-6td^}Ucv zZY+0dygHKUfl(5xDzeXhuIu?NW9P?D{@wTOF1&k4)*rH7DzNA!%s@bQ+j5;^0b~o^ z+}z6a33AfuQ*~CUF(P_WZs5;V$T5!Zyx}kaQB$X%38h$p?owsoCAKp_NVbi9t4;~& zTpX0zMpRHAFLlYMW)BqUv`$m)`3V6x=Ba0O^fBXWp211iNGMPDp|`hrgVCRLMOigI zVyN?Oqjv(yK&I8It}DueG}Rd?Us-mDQ~Gq;@M|+9b5lj%_U;1fNaOC;XM)hptmW)7 z%K{X>Dooc+-iHxQDe{%iX4J}m zoV=ulRdSnbsB4d4er`D*p!!{VVCyswElV?#;YnZb_hr%>TGwMEg|3i@{2vEPC$Nq9_5QR5+oN} zE1%1e1gB>96w5>(h&k!-Tl;zw#&V^w^P9F{3Hy_Md5qsM_iG=U_J|$Uuk9Utkz9pX zcB3gOp$t*^6FUfWBk51$8xNX;*hbpPBNq{mB)a+nzc67W_ja;Z8xLXLcv-P-rx}tR zg&K-Zc>d)n$f5p;lcG0Rzf1P$9j<`*y2gPBo0s6;zMqQl$O+7|r&nHKsyxbxgWXg3tN8Z%RxlI~oNCGX99S1JS3jOhkNL5F2U}gZuqCIKPCq3FbFlM%6*&m= zH=VshGkX_9`s`Vp+&3d}_OXKaOZLA!C3LFPhWfAT(@s^=7T*0(X#ashb@NsrTUsnu8Xo@K<`TSK56Z;_H=M;ZOYyakgT7Gk zs^gIx?Oh)hiuAs8rOP5_<3z zeR@yk-RkGUswMt}?`Ec$zpLPzu=!7HdQd;LD$tJg-8)kPdz7)&(CXuYw*ppvJpKJ> zT_v_qzdFG1iVn-!sT~h+dm_=ZNZIba6jm8)j*1=1#ro*uHMnc_FHb2U7_qgZ3gT%4 zm<NO?PQbmaCe3n_(d<~NA1=UNQyFr?a+SHoZbMU)MDRtJ30V4Fpe6L+x*P}5$lEQH- zNPp}6U@4D;jbS-lZ)>t4PPTSG=e9mj^S_!rczw(IzNBr8BRFCHca_ndvRg4+F;H}^ z*%9J=?Ua-6+(La5x^3@UEWm&Hn$1!DIn1>jrT^S(iFtR$>U8YL*rY%aJ7C=i$vf$U zc(=@Jm*l1qRtxpZf8*5KEO(4{Z%t#+uD+sWZG*m~xtq>wDSYU8Bp#otgjM!{C=XlKkd8G*(wg#3t}Sy>;D7-Kb7zxFt|6 znluJYi<9an6GWl@m@;YiS#M~NSVR98<51J+HFS}`7HTK#Z|_*EffkOH8)rU6LUTr( z>9$!Hzdi>WWN&(g5E*y8B(VsWgNp>a%d-HODi6SvI7 zE@H=rG8*!mY9u|5*4SzsjMZnq(5F20z$y{_N@CU~tk)c*bZR*DFP_#m`3dCL16j8> z$GY$?G&p@Db6nj6WZof#D>p`<<<+cW8igZd>y?iXgbzbgsvhkbZePe_d8wHzH2{Si zTd!aFY7V(wSLYI+m_fNRTc@_BH#GLe=u5|K#_a8)Ps`%EFiTIuJ#l~<3(dqsvc!%- zve@x0<2Tu`>ZUB?%8fn@Qofv#Ot}Yx&tExd)aeYV6F(FE6nCxH$2?ZwzMDgItEo%g zBUP-JDe)T_o5h??uBK46V_10GHuz6tE|PT_k0gvuU87n3 zkna{sxu_>YzD(1S^vxn5Ug6psc`z2p)0aM0@-joMN&nGr9!=1=S?E-4b2z5ZdAj5i zq98Tg*roSe(t18f(fQlM8qzl#@n`LQigjO>3s(+%fcs85vi}EpDCp989_zsY=_%(E z@5FA00<&?Wnpdr86+*oh{W2O01(Uw@e?N@AWj~Y~iUqOYaMj&L+6b(6eDgW{!U|?R zUXe^+&x_hvB^I+M!?9K;YnJeS6^Y#SmTYw_n3ExtvW>YPJG>PvoOoxkgv8M=eoYO@ z@pLJ&3$Kv;WZOr&T@{$wsi0Gn;{KN>f5&~-yZW0kr;l|~=i?CgWAW=JbJ7{Vgq3W~`$NK6AU8AyugXZ}w@R|}zfkLv(WST^KGSxfDe3xg7`T_U`jr=coj zkx6az5yZaqdvtp8BDy5sxzD7j2jvmkTU&isp!US0&nJaNA;n#3^l7oldYtXKN98q( zcT?vdznLInPQ3a1LkAi#$)!A#^G_HAY}@^NOKe|u2jFHZzU(Sg#QB_#b4 z*_8xW!TTa&@# zUPAE?8*hxc;r&{B+aD}{c;bxFIZiAp+U9b{!x4N`xXJac%UJ1>BYDC42X_xxQ?11DO-`?@A>uOco9dE1PVemXFQ7e1vJ{n!lI`B|qC-_oXeik;Tgg9N( z^95qRP;Z+}I1|-YNVcos8=F3db&s~p>`-n4ldxfhA5kli|1IK3RY^D|c=-q?7|?>@ znO<9yKJC9eHEvy)e{ie=syL*)-A96;e!u$OXKUl1=^?jERH$%C&3{Bt8mpz~U2=$99`?LKl>;cuu)@3 zR=w;oHZP=bh#xV-YWCCiuB5|QGyE%J6Qd!vD0O>t~~<)E&`Lxm0P z-3dnn6|gjLI9KiKQ><;*c*V^~i;avmqkEH`kt|}X)nfDmQe-diT78>_EF*%vP+lz* zaCRJS=eEEsrb|lN0(T(!zD3g)rV_}wcgD5c#uMTQ4k7i;EJ%JN^oCRM8bnM4TUnO; zhD1|)qbH-(5SrVxB%9p~@lxh9-8I+#@>HuT_@M603YBc9FpZVzQar=JpCf^o4^{ElwZ6aqqR@iMw-CPF*}?758F2iESo3E870^uu*{oLnQ@`*8^+P8GBvCj+Dx=?8^oH3+y*1Q zL?~T2-L0|~mxfZs7YeZP{3)Y)ZeFY(HgZsZV}jZAw!bb*e#g?n z`L|_yepr9A#LSO_4T*8b(w&d4pBo2dG&?FPpw2NuqNL*}G_Y+tA$npGOblmk_C|6; z*~e1moad5Icl}54y;tKqs~fd}&^1o)-ktefB+Sw^t~$Nts9BIPSCnH^YJ#Whj+$pcD@^5n|gD(sipt{*}EHKOs(%a~qP zCK51Kjb+ptw321oNKjF^%HS0NSBD~=HW^wY>AYZ!z^_59D3|-~DX;^JAEx&f3$C9Z zAA^_QbyY*+W^0ccu1ip2a!dL{#59sj$F54u=wkM)pe)5JS0JCbu-D|mWvJ@AbxGuy zDOAoKpQ4W)!{Y1mK^#;*5ZdFSTmAeGl-)k+zpG9fs$5bTH7P%16<@;hnROk?q4YZW zzJduzabs=AYPUeeB>$MO!G5Sb9p-BoD)IL@Lh1Hl)14iDNSv{LAv?GaQ-3rRH68kh zxtA$P36(`yJ};KjzI+?=BwAQ0UYBA11NBQ)hMTbIu+$1hBIX{sp_$Wh7wdCMFK0}6VAJjtKCY8$SozsF$BFVF*5pjoA9FLo7X4zv z=w^E?IGQ*?%`Sk|<9GL;q9$PdL6tNnbw(_4zG11|Bn7cSd~K?$z7VG~_aNo#MF{Lw zjyZ8r1d5^pxde2FAR%W^I84|bBJPJ_pxqb9jqL8fHfaprKXdnMSeQXl%8cR<|3s*f zJyy-ceG1I9|IAAXib6=atIqvbg+QLL;!$@khUzb6N-FIvKvZUZxr=KY^XfFWZcdGd zyfqhkC;bONRybVv>CYA@yS{tGecl4|=cT13zGOjp_`qD)XA!8K6MGRTO$pUksVcXJ zZUs_yK<0J7!M{8us%^)p`aZ0=8Dw}PJ`f`~uT)O;hONgX-RVLlZzR3ayM9XH2ohA; z4P3SSu_DU(`{!IetUL6`F7n0*W^?S=^YY_77C$gO`{?KRt)V9Z8}4( zirF`|967&N8eKd1?uXu&h4gM^8P=I`EXZ<5O7y&i8DABp*#qB09mV(RRgWXkV*N&l z4O5`L$bVvMOeexX?=Uvhs#Of>ThC#fyF|GA#j%=v| zW1?DqQz~S3f#_NFw<{>n^ofME8b4I`IcPJ2Z@5smLr(wpa(Ob)I-!O0R zwK&aVIUrHw`+vScgvR%#=7W4UfjBbReWAJ>OLM2?utDYuH@yn7#4!*W&h)mVcI~Fjw78IGkwkiLNksnlzvuG2rK!!7v)bs_*X)2nT zZ#H6ns$7YM;UkP;Te5#V>4O!ghuxZ!?ZL@zOQTRND-;{n$TzS9`YL1)KYu-ssSE_| z8`ca^@_=@H*UvY2Su&YmS^os1diMzqq&|QOlWw7!Y94UxlHf^lp8&GDzxeZlRH$Ax zbgo_5iy<}P(L=nYaH_)l+E=|pK(c--Ik$raN@sSH{#Z;x;c1mC_jbR>^m>m+9j3lm z6e!&@Ve=b{w8nSe+hGWfLump^PuJ^T$;I|Rxm(cYrm<*2?=Zx5W%32|SV5_@Ig8%6 zd2m+TygTHl>|dU8J5TvLm-#}X)(*0K?`tTNzcnF{mI3*rTbn|c3?X%UnySXpDhP>v zmbPPI0CLPbwtgqHVPv;C`Ms+MWaZ}^F}lJF70-qL?A2O>e8V;OmS4j5Hr3ye@34@q{UggR zHY^Q0BA=G~4xK3X+D(a2fqN0~g)oz`_!rgdGZuy?R z%TthRxudFs?Ibp+$6qonw8L0kRl1tCZ;;@VW+{EG4)S)+gNS3=Cw@X?DiXY&uL)RZejMGdHN1xeWeTsK_G2rA&!d!|yJEhltgUnk`x)VFu zFzw#8OFt54FuAvMX7qh77Ln$%I>&ZnQo_iCuPq)(9t!NV2txGRM#VGVwuHJDLMr>) zV=(^0nHIwnMv%SKW8q3VhLPN71&czJAzI7oi{j!Q41AGrYwE>CASgak>YdfY_~SUqDFc7~eukFcfaLvl5ed4OKT;FZy#NL26(v;mXr3e|bu3S1G0J z+k!DtC`i3Qm9bos%22ZO7d;3}=0Pr|-@bMD#uOVSf7oG!*7K6eiVC{AzXBMOhh) zuV|Qj%WjWejju}<8ThcMBYbUH{SrFO?B()mQiSv^yf*C0Y~Wnr{Of76GE^x077D+p z0Y4Hu-E+#LkeK8f#AL#Nch8-nTNX5i%Ev$4+CGY4u2SRK*2-{*-xTw!y=Vps?!~fX zn9ZS1Lb{|x$|1<>xt?+DmKaopFaKK57y)a~!@C_XU4%@M7Jcg#0_YUW*!5GHfYlh= z{;N9U5Y!=G>HWhP8-_XdW$(O?RUDrVi9|?a-k0LL)FDs7C#H3N1f{Xg+jW<;mo-Km z|Ek~O(geEfQEl_riI9~;sv~`Ph4J(JQBy{|f8!~#HT{4{M+#Oj-(|nDlLq}WJ}~4y z=f!d>Q__=AYa|2;$(C)K4@FYx37Kd)SA$?#TghtWg+EVAQTs zaK&7*!JB7PYUmZydi`4$a@-cW zg5xELjN8zxNs@e+!IQoy)v!X5}-;=dfVrlmJnD<^P{9|2M|>E zp0V^^1X6flCWCh_M3SGFoi4b!p094*az~*JNjIfGHO{7BmA$iO;*b$$xQOd%dh?>~ zgA^$_iclmi#H(1`xeIw0Z@9OZ`9acauAY-EyC9(oKdT5&qBpO)o3Yp5`@R@lz4vRU z`z6f#An@%yr8^khS+aig-V5`MJ9u>8&S0Jp_-lMg#;hOjlglZdV}^VVRs4iER!r~C zHPK!Lg99Cwb^|z&FU`4fHdz9Y-8o$}4dY z#$ftcdSY%HnpUh&`iKNx%(rygA6sFQiu1X&EPwmD^f(rm@E-r*iVUm1T(8Y{garGVfRLV zsS&fNUUe&7J&D~%na_CLVM2mf#E!Q8WGufnuPN7KhS^08ri9s%2BU=9 zILvkUiBcccAl!?FyL+M!O3eBjdFZsDdd}Q&_+SR4@X7yZPzr@=dYua{ZMsmZtvLo>yxgYReb1n%wIkD_C>otvsrZtm|Gr; zG1Yu2#zz=k{Uf)$2c(TCE_YdtgCAFbpeLubgHZ9>m1=v!Ff50=8z;D*d%R zm_h9!?)%acldt4xmYWq~R+Dx9ucP6R&G?I*)o1;8zjwd+C?8a8Rk&oY;EU16ThG-V z$_1jQ_t}%VPLN!^Z^G>TCh-3pp|({Z2*{&xrD;;$P%I=ZQhxmyq=*0NRTeda%BkHC zdJgY^xHP^;KOUH&?;WDwOUp%YFV11LY*2!t8rL(!SA)TRb|OfsnTT!XhW&#I#bBLq z%}_vOy$&j(K-<{+=%v)1mR0rw+nIY32$i;2c*tvNkj5R2btj{LaUFtU11=7hCu``E zb%OoU;jRDjC*!@NFNLPD@axcpw0(h?Hra98@98|=%#ZC`&Zom7%|$oOcY2sOC1_G1 z@d}g2zTAJM7lQfYS~LrHYcQw0u~%9^2{p%eZ|fd$hm>K@YOn2=!ROWN-L$Rib&q%F z6W&Y5FwL^~nu~@S#z?dHR8=QHg)X5bT4)~zu zeAGgB;Uu`Ml-|8}-4z;@-aWs*=>g_HjJuO^myB^g#U{mdrJzVa^f^PxGNdu>@F~A~ z91~)G66iecL$*gJ#ZiVpC^159k0&4S z^Lq~wE6hXJ)yu$JT|jvH**sdzH{A>OfB%ms`{lcL_AH>~yZix~+I!&j=0I0%@_OA^ zmYYy9r396##fttD2|((8y_73y1x1-t!ZK+y80XG$`ff)3(Ksb zuz-epNL(1|l{&a?ur8o%y}%celO~XoL-d2+JAl+{c^LlRiAn6rO5q|3 zK(Omz3Z3}~`4J~*dwoSA)q5hcUV;-6Z0PaympQ2U^C~mWDhmSnNoFNwk0Co+kbGlC z3ge$eHuc;{!Wy+lZ1zu_u=ZN1mb!K~R!+T15Zl{bOpoMj(t3_BArT+3%UKllUVCn-V za48zR36IAZ%lWjlB_8yVe)LZ6`Adwy&eHPGxEa%rimb9p6kv+r*yljsn^?fdaEr#+ z8M6;;`>Ci<13nFv2JN@Mf=i`A-Is-K$h*?$aYVNt2ur@Z1Vjpfw9lpILvaO|?2`Ag z%NvDe1FB>K%l-Ab$jR1xYd)0Z|M_w0;t4<*<{ETm0pfGZU7e?YgQ@RCpOT~okVfcl z76}s}yG&#_^@a-M8`Rh3a->5+fMWh0G6R;ET#x-Kh{E;W5*LRWZ z;=bhb^zgs;$1}CU@0zJ&qPxnEy<@|Wxa}})rM3i=+sW!bvzkL&r8hLn7mq>8h5HY$ zKT(C~gy9TduUX6pTqND2-v+S{*4L8#v}}vze0ZA{nW(wgV<@a@6v?YXUuzHt>dBu zSY{C$sG}>5N$FE^YmVwj;9e9IIr&e=3z@~=J>GgBD-7G3t{Kat_VI7}yqj-9MEbR% zF1bOh%-qzuoq-3fBQpg&!}~B^<)9+Z2Meq|bX+ZWJRhS?4hTjDX<~R}qfpqW6b2T} z4tK?IL(mGPa>Fu=M=yyBm4NveZTL#XxIaxWg^Z)j1{abAG1 z(88vil~Ewb`RgP)%t4I$dDf-maCAfO!}RJ2KwcKoGkRo=AxwsW#fQ$p9s3{&vH0gu ztvxrMx}6%mZj2Gv2z=;uc>C${#|d~&%!{gCXFZ>nevvX-`Uso{m3~VF$N({YHFMcN z6k8+)%3d7ILIVYkWS1f)&@W!OF8f6Q5>!u$?-OM}yNK8#xdhFB&x83KZXO?!$CJtr zWNKf(hJq^Mez8syC<#jl4|2W%@nwd_EFYY}nUv7~vq~IY^nbsv>IvFeeSl%r!4 zgp7399Oe3skqY?~5x&E(pVkjq$v()w7JHF(|h=@z7ObdHkRsKImpgZKS(zFhy{oNG^6 z41^lc4u?I33 zIwW@KzXUcu^oc_W41_2hv>gH%x-9vB&Ufb6Lnbm*S}db5nKs61)x@1rb(MT4dlHZ~ z;Q;##y)o~s`B;{_eW&1Iobs95%=N3WS4x1*&R-Kp15!d4&S8D5J9@>@_^*ZB3I*{KPg;k(JNTF%hqLf(>Rx# zkX_xIK@dEz)G%$T*tDByqSv2qgRE|`gTqIhL{&aRE#mxe{J>=D{`}K{f)d#)#Z+}RJ&^|-M#d%>Y z{#kAW4?2C~Q;6fNMVU7UJ1SlMqzbYLgjY=km@wWHW~?*%-Czy-Wm-MJc>D&%^Y4tf zuno_lNUw#U<<5DcLEF_-X*tmCSl9DC$Cl5Wgc;Kd=SGvq`l zU%uDl6Ud2y(3X;paRdmDdv*GnHTSun`S9xKuh#z>4iQ7k6=vR@C$oK*ez%5Q=fshl26$RRy_P*Ekf(G{NfOc?VLG- z!i3f)^3to8JLo1=^$I%Uvw_g<8314m+tqxgyCcd_O`m*0!e7gzL7>od9ep;X5JskFGhXC zc`$>QbB#dL+ww^y|CVQDO@T*M{zBLb(~Av&9SEjXizTh>BB`(k{QF|s#Z=^~CX2rF z>El`uWhZg1pH(69<9s+ZM{B~XWm}`P=g$QS z>FBK;v&{S3i5~jde%s`{&uC~oyB1;o(6aW$X|k_gZota~fFD{Vj+%5k28lkcwljhb z_uTs&R@h=%sjaznt=y+*%61SJI23sd97x$|wMr+YI5VHt!u6nry`U!fULf~xq&bM* z3$K^sYH;bsyLzclp#A;BbKuv*?OopmLJ4bkn*vrEfz6}*s9n72ZwT#_k1lFK{=_M% zY}zB;_c!V}=htmIKHIAs4i0%;kbmY`8_Q?E_vlihKYiEvMiQ~0nWQM@4t6?-RA@ex zU=9%pQhr{`f7dUgVnL`;PQO*flD1&tA-AF%{rAII1k+R{E*Hbb0FthSke39;qrW|A z5OYU|B;br;H_0#wuhtB}bq2Zqshpw4dj4>?A(8G6H#(X4AJ$8D;oUbA8K0j1&B$_9lrk@2bR{V z|Mj5cGySTUjS16v*Oi}y31LmiN%dOO*G%bv8m{>B08Za*WnZ8ix*s2Y_rK9jtlM-K znH%_`30Ef$@gO~np{vIzbr#>qX`#3p^QC6VlIIArh#5rC+%d+AjUEOYxOk#tdPh=C z)8t>ou5gudo__LiHI%%$=bWo8XK(&ur?|OPo;q~RLXw-0=@JtU#=|Y-pKI77lUFTD z%G3)x@SFFE*DsyC_CKfDim=;v*hVVuJ>ZW@Nh{6ttLXGy_r$a$lcwv6Wo=Ky2Sfb< zmLwu$U|KE7$6GAF{&bFl)@X2ur3(fGA1wT=%Ts~Qa_Ty(Tu|+PCQ$?>{w89sT+IMc zIOH7k6m*U)tHtgeU7SVr(R_U$pg#&mYh>9!>)4CLy5%k_YO{N^c-ngB zY&9Ki7d@%r6NXwLWRnyxSkJje{(cRqZZ7F#eD05m_=!Y?x7U&dF3-#gAZ_i>7{qw;64YKSPB+9MAqZ;-Dn37tm%5 z1uw29tYh?WO6yZoc(`Y>t+SJY7bPsDeYa7dW-qdr=M-q}{5s4m8@(re^sn!O*_^00 zSMlFK<_k1siZ-VyG-uka5TO54oF#MeaTPxxvEXC+C~5YX)&4k65G+=xd64(U3E4XC z>jop9`N@bDa@aM$%;Epg^6xIJDwb;@N;ijO=48+-QVHWv{Iv?Y58j&eX^lb^Yz!d9jLPukUq-d-qo> zQisI8p5e15@$Hv2&$ZrHjGYdM$gdDuysp%ISL4xT#dsx8{8sWr01+n}B_mJM8-X}6 zW^U}M-aaFFtLK0dR};OYm5hsBxQ@=j(P6AR zHk?T^r!#u)#h0$!^@&o* zTBC`q`T*&vEX$9|D#B(mATN0V3%zl+#Iwv{rjOz9jyNG91e1%u^nAReDS_*9{%SHpzL1*Wq zJ9~?l+_4DHG?!&fCHgf+r|Uk(1Wn}B0O(g>Zi{h}@-Kr}=y!kkZU}3_->*&QR<;2j z(MvI-6P{qxi=6BkA-mXlu_-P0OZH~NF~*59&o}m7^i^B+qLj4Wp&o`c!MUpR+~8ky z=_s^4=vosOFFooRV&ay;LGQ=DkM9f7Gxt9!`fYNd3hQ#5#~}8tMvHZl_z(|HJiX+j z{!~p3?ft+k+5SU~ETXAK`Ljj0I1dd7`!)ZWqvpCqv8z(lrwMZm8pHV_=TjaDKAEYN z>Y?KJxHVEVj#d9qdR3O8Q|CU!Y^n)HX5(QDz51u@sv6`KYrtNl+KE(CcBFz!!8Kedzy!Y zk--Xd1m(WqR^>KD)}HA>nUhX0W@jk1p4;!do;VIUQjH&5)~*#O=OWlz`{=aJWKUd zf#ue0vpu|>Asmx+2QYK~G&oXi5EoZG2AlwjUxpd2;w(-D2b!SLW5qk__y2yLP|}JFWOlWbjVn0<73$f!Zt!k_vwer2jy8U15iNbYSVFGuPbz|%CWhz zh)%O{ID_P=plr`HO-21M1ll7Uw>>IO;#tZ>M^zG84c1UC^gvSU81b|WNE=`+suylm z&0mbZj0)RxreH9RAiO7ZD9zrzI14$g7?1H90lnGm0Op6hv|bz<_evo=jp-IO7H_}_{*bmjwAiel_xjv&~z<3 z9;6);8v8ENGl7(BamrkorbEj#5C_3`*v`TTzY;tbaoT>@h04<>W{34SM;tnPHZ09Wf45P0nx^}aV4U&Y{53FxcQJ3?=3 z%3T(o?8}3^m=Z9IbkY=QmH>qvmy!qFW*XHvYGP!CvWS^ov!ZLJWWsfmw3Bul_l(mj zptVZmv}TQ5pPoKg^>=YLMU-L2lc&xH^gLfqD%SiwPd(UNYwyZyN~`6lr`hzb+xC~! z>#C}Td52pHu+B(r{yV-8Y|^&&ZN;I_Ulk?ea#qHLfaV7CJXMM-ZLTfK@2M&AA)-_d zfpv7%@GG(bE(1HS4jXJAAL23@pYcGCh}5AyyBx>3KYM1|MLbZq zjZ7)+f4Ew~iRv}0qrEFr1R@Osg~={&q*>}T@Saa7h<({~?`u5LWKKDN>kPSV{i|z( zFlcHSCz~CHD6`EoA3P?Q<~1*$EJu9?+XY@k;-^XBY8JLz6S9~b!1W1MN4z#(Os<}O zy}4X-Dd`xb8eFx-F#;~cVF!1Mkx%PAQ;t0yAZ;igec=3z6aXD|u z7<`4f6(sRnpFLRXS3icfz&Pih`4sdE)`-C>#aQOhTiq(E|LGE-=S{wLF||Pe?@@>L zJ(f>uQ0J-JO!vK=ytlu7J!mb{egm+IP8~nDLwIu4E;ia-Ei}y*Li=48B#m}*tq=MptEM^g%l6`P^10_F_IwFsy}(^vqQdIsiPNZ5 z!AdIu5!(r+#3#hNnCkgcw@=Vq|01JHZxHbi)MoMC^pFurL{7yyLfl<7n|^jv4FTHL zZkQjC^!WEySd(55R5xedT6ubM?LOBdI2BFr|0v&Iw^@O%k#hTRb=TygQ)zhF(OpVB zI-@N88qL}@P@D;Ek|0+t$URxoGV1py<3Acw@2ZVHXyDR&K}k82a42P$n_YMfnJE5R ze#uSxIFV%iK#!gHPzBf1eY1ZF(brnf;^9V304bN&*Nh`w0Fxi zm>8asx?m=fwMK|&Q$ua5gTzSqT1{>kO&?qB9f|E;eTQ9;Dk+GS8=CcU_KUv%X1$e6 z8yjIO%id^mMH@`MBD-5cI!eC$6FC!H(iiwZT9?zZAIvI88%(vhV6jzCYs=8>D{tRk zPy3TN>tWuk7hBYDAxd1dOQyjoffQ9jN}09v_p_NiZFf|}8Z|d9I}p(YUl;%{{YB3N zi~7^r=0u@_ub$e2HO&HenoURM_Mf`XwnKfT15-djFj zVrVK!Xk?9`*Nr34hqAFK1?)eZU5B6wMAPHcKr90s%C#nQxwpB5`J2|RBy=X+2t_t4 z7%IT<*G=dJ#FEqYmK+F?h+l^cK}YH9NJ}ocvOBc8l9A3IreW!!ms@b899J@yQ9=t9 z(P_pYZZXc9jKc%`5o-dp006|g03-^|0sw#y0N6bOfOI%J06+p_BeG&W?AKjCB{iKN zm$@bBIwE6~M2(0MM(3|%6CJqIojBG3)x#xspWZVS+2%65{3)G%J5(Q=L$pBx`o)2IZP|$1?FvJ1w zxF0*FZsHRDx&kfdG}#e;GBA36e)Y1%PMCClo}x%ePaROi03RDd>oSlk`4x}CWgHIY z|2H2fStMoU0?dZjL?QM@D;t^@1y4Aotoz2)xSv-USHyc(v?m)W5NyYuSOe)0*^ z4;|llrG}vqIKP2e%?n=1N936}U-bMaz9HledC8O5)4 zar-rz^IUBi4+OD69fUHn^_+5SL=es0Wa_x`QF(;v03|SQD?~r+#rjnM_#^80$ASmm zQV6RF3(t`3u`*EB@LOH&Bn`}L)Rpsmxv>}n4$BGqtqx{As^kB`qyBKB=Lf>G&NMVj zFI%&C_Hzu#aC8#RwtnF5kwWlx=}!9v`!WHk8!jK~soc&x6jWVA2^W)ymclwGe21mF z#R<;kF6vEn;U&jeN&_~#ly?Lp1Qe9qV?-0vUV479D!hj4BGulU2)50DuZ@2Nl8dtH zoDD`PJpyPtu>bn)tS}5-Lt60w_^nG)&eGJGd5NSEXY+G%gw>&%TRW4FcT_Wd-nk{> zDr4*c{g=YP8M6{(>F(r4Q*bh+70iLnFt5?ktOb| zLcjj5XdP4)Vlx8U0MKGb%6d9kGs>1~)* z*E-65ncO#>Ys|ODEkqrYh1+{KIaI~5i@SI*AYsMY+x~Pg(wA4e`FrXxka;D8RR5{yrwyQGwEXH^{pL0uBe z@GiB=*EaQQg{jnE_?U@+x);F)ziDP?az@pP-n$@N{>0#kWui1sJJ{3uzGs($JX9Y1 z=n+1d1`h+07e_w_b^e`^0_=+|SJxc7`JFcq41&DEqA^kqcY}{~16SH4vSwoxOY{2a zehWi?Nr((o=xX|(N=Kc(2Rw{+Gdfz*!ixG3+r~Q{^nN9x_fi>{lnGRtSs zg815p!?jg4GSrjxeh1luMF@e$NsYMtaY^w3nCci0Qc9K~dun=xG`hdT;MauT(uyeF zhkOImjMQ)3^fCM9mte$96ui-Z{dJ&g8mwhuqUT91%xLy-d6#-BKbS^n=^kGYc8O2- z-l}2aF{GyLkk;||IZasoy$q$AJJ?l!ok0WL@!gn;=;~gZ4Lkv)LAGY4O~4i^5)+CI zRo%nex6dya%p$vZ#8Mp7=7}4et1()2XeZL__sQ#^8mIlq!_=;sVt83zTOz1cG+VvB zJ^S&?%Lz-X^zzBmKq_B9$7xmX#$DSGmu&wGHOmx+Tjk*wTdSj?Xyc0lCH=nN?w5mO z^A~D7q`X2AEpfHV&2Vn~mxxnJx&2BmN$r_SEfJPm%&a~gO=o0nXRrBE*2~iv-Ee1I z^e#0}5=$2@X|!t-lLu1^4TV?GODWie?FNT?eiAtjjSchMJDB+7veKSgDeK+!_r=Wz zuN?W~oRFrKJ80jM+K#uKs!paam0Ji}^X3WQ*O$#3s}{%X#=q(|co1Z+)qjc!l$jV8 z1Xad{C)fI0`nPV`#da1pu=Qtsr|e7%HC4@<2T#xlB**8qL^VfvbjM3C(7pmIfeIAM zq-xGF;vY?1>&GxT0Jz@nBZ@LHWHT1|=25ufJz@n)P91hQ2m}f9X>uDm#tz0hKR{zT zts0-reqEV^`Vk%a!R+`8FXbYON#1KMJN;RX>rHi%-+4L8UXNSpE7IEq@gd1=%+D+P z^$7jLM?C^bG`B^#C3PO*G?V?v=|fYLqyn?aIg6n!Duw Hd_evMA++@l literal 0 HcmV?d00001 diff --git a/data/Transcriptomics.rda b/data/Transcriptomics.rda new file mode 100644 index 0000000000000000000000000000000000000000..8d356b51cdf9189b123f9def302515c1dc3b943c GIT binary patch literal 22187 zcmagkGq5NMv?S1D+qQL&ZQHhO+qP}nwr$(CjsK= z&Bpv@tE(sW<@k29hwAo?`)%goZOx~sXD$A1_H9M;{dUjOt!6j(_RV!?agDWQkL4nF zPDif(6Q3Wz8Q}i{Z*}T@?y9TlV#6xMlW7`-Q&HK3c4HZ<^Kxsm+3LByY3uRqxpK4Z z+2k}{U?}wm87S081_05Qp_e8@W?o)m2FwSL$YfrUxbUa)BFn&tp$A!r2%dD8VPKw! zVFu%81XBW5k}2La$0`ru6OtCF@MliXG&dg9mrf@GnV+xb^uWrr&LoqDjEG3#)(40VB^0X21dL}!CNoDimj;l|poi$st+?=y zVIGiMNGKA1UMS5(W?*DaU&sdtB9qRv2wr}X(?=|l&Qg*Fu}CN-TQ**ZCUd%HhVaJ6nteot z3$`*iEKu5v-O|3%Uu7YA;snOP%Rv|}U0ON$j9Wa-A3?ma;$Aw(F@Hkf;Q#9fK=6k! zO#BuELPG%%wXZJ?1|R^S9fW|P!@PoGC}vnUiKw(1wNW)*?u zfIR|9fDiz#7F)Kc&#=NFE(Cv+cdT-{Oi?z^oSQtE-=>>{-Y(zgCZw+KM_y5vEBgmv zK9}op0%EYq(op@*2hUEI&&j(kwMA+bD*AGj2Sv_Z3A}KtvRfi=Y30j2SVp5TmZ#nj zU39)NLXto@)&Zbr7SaU*yZ6Mm2w+1$z41LzpScbv{;J$~sAzMJC)l3rQ{lMZA7`0N z{}Gweb8Uq1Et1VswK+!%hb2NzmZr7wUF{qxO(>7oHh|c&S!X%b;}`uESKqDg+&qpv z002zkwQ>AP;C}i&3B>7u2HKT2_zBRP)S)$Gy##HHN1g}3@oK+c;;DDGyZn4*MbNR- z>i$wLD4Mi;{Mq%&JS(%$YrfWbf|C#c0B4om3#~B7>d~2StirqD-wp8a&5#LP znrb#t(Ro&orF*vpqs8PfzIo=vvV#yHEq1L_ST0;|V%b!O6Dy+tn-;T-(hHp8iImY381BiXomN>?L$(%lhDnWty6sdKO`dcJ$Qb|ScI1fFAqB!`X_GSX zEmwAA#w7C<>lJk8{8qO_xfMW3S(J@_2ZN#GA-L#dC57W!6DkS5ax<&QH2SfANwm)& z=xW!aza527=3%(A?G(hvJ3*1X@wGHRKP_B;H(_$*nK*>=5ggK|sw zzcM2;L}vh6jw)zu2w+QjFCIS>zdyNp?Y*@c`*9o{{;1*>nzAucW2wEQWJHq#W!CS_ z4D%OJc7h~yIT5=Le?ZU!zyM%iFb52|MXFxcBKvKG$eMM2Sl$&d&=DEXS6m7t{QZZ@LC!08y<@f`t1JL{U?BusD8XSo@*BTs_alw4tfT#;8l_DF*SPMBXEBW;z0Vfve zZf4;}I6a^>Kgb(2X7FhfnkcM{k8_fmc%)hovj*vuL65eCog3%Ph)7#5ZVMz_Ek&*V z-lSamU&&4dDrCm1IYh*A%S_8orLmsJeSa7r<7BO3hmR%Mri=EWj=I%Gow?Tg zAB1d9)Fa|ysdf-DLyI-IihiL*W{48pxt8)#$CP&eMat^ytiQ=vZdG?S14<~65;0ZR zN&&;w#R_&^4r{Xx-%I#=<@9IJIC!lIR!rBtl#g=c67wp6=@`R}x73>h00#I2AUR?I zDn3k$9I8m%G^kUekxtOHVEL9q^&Ht1p`?-e!QfsevT=Z`JB(xXz1V}+#{pS;5{*=& z7#sIuh|HE?Ze%d^Gg`{!k^78*Ab%3YvfecHg1kHmW2Lup31gfqDREm)lJhvR1N8On zls!Rxaf8be38P#%vRQ8vf9oqs$Zo zBhbf}CaXzIoeS^zz5D2P7o%`ZNO=dEw0ru`M|+P1#d?(l0jcQ^O>{LT8;_T!NP{;U6X2ak49!(N@Bd>LX?~Egf$z)Hn2n$?F*Ht8ZkWqW| z*?&pbp^^Z&-BNY*4v_wW|M6f0``r8%ljHVS#v-QF+%wL#8(~OOG@A6ExtCJv+rC^%w< zhR_j&NX~hoC6^5`^1DFFTFl0o`LDy|6I@syCI0;q$T8GKR*XD&+ zM4bzr)e*iO^!xechUXyBp5%_X3z2vlDSDrG8#C7Frsq4i5tkUP-L~@4=SP3d-;Sg- zenlQQN)9;kwjU;!Tj9qVr3F5^v*n;mX7s1F8g%pk_Us#VU>pqY((E!_eFG{GL!r5H z_fRJ;{ho8zEPmXM$0(qMc^cqrM;$kJ)Vj(%DwXl?T`CVp?C<;kKzr+-oJKYr2Mz&A zQINnGLnh=!LHmZTP1BgC5U<{KZf)c49;;%)RV)hX{bqfMUNcE>wJ`SAYmVPM-i){; zls_$t-YWAgDpY3zXGCfN9N5_i3WK~`sI^(aL7mV+2ei2HYb_~rggxbFIrj3u*-c5# zku6G$xxE;OhBvG9yxs>ivw@I8hIo#C@GxI${^K%T>t9*-kZ8Z{oitmxF6$Gi>%WAN zc7dH*$?*d4T6uVdW8Zm_XC7KT4t!LXX_@X+3Y>Pk8B=5GOpqy}IY}VKOAGj##PZ4e znJUCS3NXCxSrBZrx-_bPW)T#it8c{X1W@J0xs7QHXTtbXe1{=}zAbk5sg?>7hLOJw zDsp>H4jG27+rCoUVH8^#o^`^gOin?NFdTRE%fJ}_O|9zq4*z?pTcs4#grzhtEDF9;kXb62@Vd2|wHB3aVv#WsuH%S8K0sz*IO)i< zgL^lpYAiSP6K|ikNJ!Nb7PUC;_mPh*%1V;&j_aXvG5K{qhTTXD zRf0W+vb82AiPOxUNX07DBJ||dCAJnbLN^4t66|=`p zIO)n$XP`R1ST3j3{s}0WFj%`VE4Y;E$M zY^03Xq@pVD*b>l!APTQ303&;&iTaxxt@QUWf*WkMPlvy^*-x?S2fg1p9e?Hq5ksx9OS|9*^Q?*; z+7o))*a*j_H?%^ge(3W(s3YNyhmz_-w9k^3mCx;adkjg6Ip%w#U9&D-25aM=9TD;# zp)n~!HyLkjYE@xsrbV9#d0N&6KAd7pwVMeV#mTJ9d+)=j4p*My7IB)TbrOUJ#5u5%* zdzMrCL91hR;(=F_4LeYo#uAHQF_t zhHKf%3(2n)J4OaR)wl&Ev}Ds&C!aF$s~=C);DdmBv%JqOlYFJ{D+g`Hv*LBk51gqJ zA*tEdt6Qrj5Rn=hj{-n7qB!07S(VfnnHqsa8Ge5_ji-s2ZY+6QR~Ej`eCt#QT$KLF ze_yrnoX%Zyo_?#?(Wg~6EfejJBKGBI+`$8zQNJuSjFqF5@c={J_IujDPJPjQ;%mFC zQqRVWz+xf4?A)*Jmfn4_o)Y&_8cdDrm~tB*=c)>ee8c%9CnxdJAX}?=@#GPtSNLVK z`0O#ekug`J@J~VI90tPHuFFje0+z@~+;z0G<%4mFF@MOVP|OY@>g49+6jq-j5*AM? zeMS?o6QK}c3s_j}S=PSYJx%OSeN!}gZkYoyo&!Xz?e5>X$@)bMrXYKo*O@j5-qCM; zoKpPzQRGw`7+lvatn;^bb0fRHKEJ&MfSURdjOK=0pNev!p4kuXFVg2B6xl36wu}$= z@YI6$a&7G4^Cx@&K#};q%hUa3{YyhLJUWR@*pgM-l8LNF>flBLlq0&jGJtLBswKPL ze|}b|InPfq(dU_$EU-{tl>q{34yRhpG;+qn%rsQ*2RN~L^am(v{N}QTWjDFqB638? zX!P)v14ZZ2IPnTWsZCpNEHMA2!#qc(Y8(`gpuHeHgZaz(+KQjFD;(TM-uDMU?O2zA+Jz-i6MqGaO~I*%iPPSXp zYeL*dp%U*P8j0H0@C0t=$7vad0svm{2W-HAg7||4jtUb zV^tDbg@4=TbmX{*9EB+X)6vn^@kuY`MK>iPsga(nB=lIyl?T9jIUfM~?7<8+Hfs&v zF)UA&2Ad2SSbJLDZo6vp@=SmN43PYo{KsD+VPm67c_!gYPFH2?_6Q|n2h`qFRACrf z5~g_A&X8$$4z(cW3Y2 z*yi+;LF%O=W0m$3DEN%NGh_#;AXbFEe%-{xR8Q`BxJh94Omwx$nV&w!durm;EJn=% z5jLKjRV~?dRmXyLm8#R-6VeA1A-NvOV;1s!Td;e-!;}88Lz|}X%T0W3Yl=p$tIGx* z0)?Vo%|ALfh>=pw6(}e%lV_dQ-gh$g>e;&RPegMmwwDIQg9}u}MY!YFbNvjU1*Zl6YQ%2!ImR_KjWSttc3>j-Qy9F=@lnKlVK_ew zE;7|Q&KXLrC>@c39p&vgajl?Bhsd-ybWx79;8_Q`PlIgg!p_kkQMElt)4(LTg|z#- zSY*0+>kGwhc4ZA_>UHtbMT#4I9W8i~_AG8W#SiVIj-!T?+@N=)u?pcAO1k1Yw67I% zId-0p$)FW>6=IA^Md4!|7g|{m-|H9m-*KQdAF64!OeVVEK#{#=#*`a^)uX*1tnvAp zO{uZ}XmKDKRSFWd(EP7pV8N48zRb zMa=X*a|Rp+0^d~oP+VJub-lKO$|6iD|EX<9RRvCgPcZ8FE>j*Bxt;LFN8~4&j{{|j zOc|?thUyX0Sz_p zPXwJtGwwiZ^dC)Kr@I+yh`OP>r7k1Nrtj+_QXxY%A7%W-7fGleDQ}$s`@XGV+)2bm zQtA-V`E+gs*to5u`vgiglN|Sncq`K`T(q(0H3AECWl{--T26f@TCwePf=^ZzyAlZB zMU_R#xx+BV3+G>2$?;L7`mCBIAx4ynbS_V0M`-BMAYndQ4qO&q$<2O%;i3y$PZ-8- zfz8_B`$|gdOf4c=T58)GK%yOOV_U#`ncV#3elF}&JMU-BLK!nTe33Hs-jw${ zu%WLN%=}1Ext5g^)?M%qRSrVD)tNU2Rcw;#vyr9y{s%inredz{5=@-Y5y=W@7PiPbCCuWWp$nfk zIo~MRi%Q_5&+$|p_Fs#*N)r5-TN&CWK&B>ti}V-15x&k~4TK1)GnDsNJrUM+icBew zdO3gUDedJz?2+CEj}b*W49oA(9_u@)U*w5#;*m3;CxxEiOBSm8xQMSYUL^Q`W^7tS{Eja1F(xCl%=}L)S?}4>E;xzC%fI%V= zt$5YcNSFQ}Y1i~fpMtOoody|hy;;DQsF0UWNCpInLnz0Lk~VjDRrcl3Z&1gvmLsuz zTDCTmhJR1b(RM}I1DOmtPw-?;i~Y?W9S3u}+O0`6eJhiyEE%+&W-B;Pn^ul z0z?Xf|>a<#o zJoj7aEamgx6cK)7Z|Tk&f#_^UCPzznF@=@D6(uqTehpLRjcczuT&@Uyg6G(dR%0@X zaQl}*y!1bUkjmp9{mV^kAfk23!AA`pJ210`yBzOgCz%(SHmDF>+gOJKAO8r(y@Y?x zcAZer1@9P`#^p|aGdN}r;8?Vte+M?0b1PgPV6hi)-jgNpCxV`;?Hr7D;zXS9-MH+# zazFGI?pcXaYPXU&gYarhGt zng1|)iR*gY-;z>%DnH^~=LDQ)5|Q3L{UCF4|FWs!o`@F`J`9A_13x^UYBBchNtxXP^)t z`%nU~J=t#B-|J>I(6g^hY|7c68TvvUm~tX3#s5_P!47#p$tzpH)n+(gz9g0629`hI z(d|#7+76tjYTo9TVp}Q4AjZBe8b+tiO4XOGR<$28v>Deqv@+~G!^oPj1~ zPJw*~%#7tST)N8AD82KkbKY3b<WiL$|`X zOA>Xur)^lrb>c?-qk|Ee8eGMOg|JE2<64Qsr+WMLPX>kEi5|Re3F%a5o5;J*l!8Uyq3Qj6!lW2)+ zdw_V3zb9NSS#Wsg8wk=vMQ*G}h8Q zF{D(KqV7!6uK4t&7P~UJ9>J^_37*!J0YhKqbURyW^#czN zrsacZm8Jl{*>06MLQbvE+u9fnmLnl>hRU~L9B6nv1^4m^$UkP+pjWQ9OGd*f_}L4l z&aKe9>nxQP4m)_@ZeLLdX~SP>BXkwXHiegaCg<>ZT3)(}QDKGd*b&UWsHkC?$5G-P za14w0cxEJ%10qeg_m zpSSpl<%8r71sOkXfJ0~QAY5HLaB^LgLve^R3w?0u& zlgXHa{;{PrEV2GNX#^%;6&Hw0F;L(QgsB$oz0(X-|G+s%!UllBs4*>N)x8_< zOx^)cNx_nRny_)o8e$K-w0QTyVYpALDkhy1x|PZ!cj~UKit}&Ypt4(4OZ)Wp-8_+L zTmS<2hnmk15eOiR07T%A2m%mJ03sg%9rlzJ>@=ye!4KAz`{IA9YQu4@$C)F4I^*25 zGB>g*G<;FNqwkC6H=;<-wEwy)*=gN4jWNe4hAl! zr`=P#DQf1<1iYx+8KX){LyEp_OkfRkncZ<@c!26ufx+e9dDW@_ti4u8Qc5bI(?%Bb zqvd)u;L2@oNPH*-gLy8MJ<4}!{S=`k-<#=BKPsJ>HW%!a0{4sK#5gFaAfS6kx5iva;H~U6o-CBiTvbYmvWGi_sXcz98 zFpi?c;*{By)WT=1rFYa$?1Brr>sz#ARaj4$NUN<$!oaSt$S0XTK@JyB^0+WpO!_em z;D+f7T$v6ZnV;f$lW_azZ9pwYS0mzJ<4A5npgGd1bz4r01S2{ep9aV+)gE(q8& z%4oB$ySl+;W;dhNV3v~N=kiHbO{N^1eM!6v4uoc!gtou)*_Z7n0N^kHJbIjoB8uw5 z_6!m?4nKhJT`Xzug{^65qYZ&Th->qbb^BzveTE4nD2Dm!O##tvs_#9@Dsv%||1BjYO8mlu9ZgbNyL z=r}bkL%FbS*}60w;_HrI7u7SJX$ykwg#uaw`ujqvn-Ilw^=u9UudKXnuJYE~TlHy^ zN6alhm8Mo>W5Alzkha+nqvX2|dEj*kMy(&Jx8heIkX75z&_K=Lh!Yjx8}1blay5rR zhX{pX@569BBevJp_Ul$~E1w+D@>hYlP!%=w;0I|QUa+Z**f>;l4!<16nM4vR>d`DN zEU!awyV9bnl83Kn961R&W&H0Jny&NfCUhk%g(0Df_9Y#VaWQhsx_$KA1hfZrIhXJU z+{bFu{PsyEA=6DByrEV?JpJa_g!xO88RD&k1>s8<;8Mhx5*RS(kO_TC_?*OY@_q$t z&d4#39Z|xxcR8Zhh^c0_Wh}#EztGDv@G^Z@sldlqW_w?(%~x( zLx;D+Ql#0{m2aD|tZGF2bP{I_z})c;UD=3@+v3`lau27B`cMPN5eGWLA&Qn-}R?J%$XvqrDX5OJNy zA8W4U+;;f?22st3@t&2#@uc-?zEF|z|D&2l;AjD}im6RD6i5tGkE zUdr@y5%@{!+WL%)E|S^msG6ZGij4vZ<|X>7oJEsoZ=UEG%m!52l}Qz(gz(W zR}dD%AGnSvR1<}$f=pjkC)RbNvhV`)Q0a-w{E0b^5m`Dt)~Jf|*rY>2?6|oMx(>QJ zdLXny|0duT@bIsOk_`PGL;cuZ^s`qjpzGvbQwA8beXs7zD=bNIZ>z0Rd}AXG_MxxMS^>TKoOF) z*GY)^j9s>-AAM1ECSSU=_-waGLf2UspTHbP>W<-V!DUalK)Q>3QUTv+$#}S+hHgyb&jHl5`j&U6QounjKGyEo<{<7YhUPR zoS!Wi5av%@Wnjr9kBx6t|0^5m&LAp2tIEfvv%Mx&_+8HB+Ey)-BsF~FW8Z9E#PjM; z4^~(YM*RpfOo(dM!H|v+%Cczy*Br!PTuT1wLL8qGu)5Q=8swVncR^W6LFEyI5}(I6 zmO+kJR{q5a&nB8r9nqP2K;ksK@4<&6#u#2UL zz2@rT5vRNlo)M+gv+QME`@(nzS*e*(zA^M=}Bl(Nbvz4jG|O^J=7*9*+n%?G~!(lsFt^ zFmFuwgaW+?Cx$?6yh zc`QC`e_GRbqE8fnd1vPoisnl0!;LIgjvMJ}C-!X;ayu*JU>0XZ+h5$&92+-WtDX+ZVM*m}}*$3rl zNlIN$^Gk~X`01B_1Z@UO9`b@@K8<1AfYZFMO*whK277F_#<>tt8=DLDAiWfZZk2~Q z#~&!7{|$e1YMu{My!{N*W=OmC5R!cFBpwp`qJ$Bbd-)O<+*EiB050ABr8Np^Bbnee zo&}yQH-PwZ!H7ZvN{UXej~aDxYTL=E6Yd) zH25VMoDTY-a2(z*;tW2zb-shnf9|^%Mzn?3xGp)48Jp$Gt=NK)KQBspKvbe)p}kV> zlWbdScq6xhr?dVlD9h~eqd2)MF=}&pJ+pVp^yoJv|Aa%pByukG2d0X#V3waf`1M!sX3a5xTqfEsd0I9yYmS>S;6C1}iM(%dK`E>2Hnt zN9e1ONh^Bq)j`|HvIi2ekfnAsa9DZpff9cuvkf95J&C?i3*?*ag2p$SdX~nmDP{|2 zK~il5{G&&=Jn{?x3;+OhKHY@mh;iE(djV4^*$`FazwITH+#aM>Wka(Y7p?_{J9R_G zbV`zu7{LEoyBng2^?I9pXVstyA9rhVH)(}U-9f-bqxB{T^&*CX^McjO5VTm>&EQ*J zPEg0W9h*@EWhmQ-HV4&>wsT}m;%Eduiz5Yhwe1OzvBMkiE!-}`#pCBh9Gn>5tgC;6 zeRrb}!J9}I0#I{rs{A@i5nD>e-=EKbbmP54KCbqYce59+5TSFCi&qOYg};2aHBa_& z27Hhv5k5!xn~-iLy*v;VGJEx1K5ns%P4iKpdbzWq$)lJ9BA|`GQS{|UujSdDOR4jk zCA8lwSpPD~%f}51B*E-dl(qO4loyx4FCMG%U-_$@ADaNBKcY@t>AA#A=|N|8JeZ_6 zd3?jdHuPKUoYMqAk$idydVj)3urtfS7G&}QI znsL#&l5qb4`~lxs3!@}!jGa_MTu%bvK6JK+lmgjLJOwuKB^Bl+wu6nR4@YlRH|QHT zvFK+ZcbA_`Xeaff4pMl*`XeNHQv1^P&_kPp8X?q#D6ons)Tl}kc8&( zpVK4z!$-hmHLr*qBsm+#^RC!)8~%bP0y&q(Qr38j)L1Pz5_I-^ATGtkB2EF5SkrO( zNN-YqV)h=TIE-74^tz|ep|J5OTG{&GGL5M?fy7UPTjbvAWd0Drj!8&cz9J=cgl zV;VvXyvoU0GM;LaPdFgp*QY;rV^f&^xgronwIm~<6I?u;5gcif6En!esY7&Ybr9}B zV`+=NWGQsm`nHJF_W8;J#fyMpkoW`$0UZFGY(>d8V?SK~n z$4{?x&!L+cA1cmWuYo{tPGcfu8&{~}DgF!PX#?WQj9LQjmg&#QE>9z#U)xAeQ?Ty> zN10&N{s{oTP)#EmpWlo?AkBf?sa&-$>GWbx>qV`eJOqX%C-dYCTGywX8`M*T2$BAQ z=NX>bYqL25aTpy>4X1m%7=^5fCi68~w;3LEY+{>}_+4o0xX(`MELT^5ozR_gwC^o< zlVHg^Epz2lKf81jgX3*=aKO~*ygk$r7YPlM_bhC(t&jt?^1gj7W8=95rKj?ewPGne zE$@<+!adbcgZzVZ&?t<7ulpLNUAfBD!E#yJa19Tf!-K&kJY-Y0{Y3$@S(vg<0=7sDhG^pbk;;E8c2W$SALEZ@fl9iX76#2H??nQ zA#-gTp23OKm$GtZ8jx@EU0l(^13=j0P7(uWUuqtWoA_a&R&OD6U&AWJef=M_33#7T zeI(BI4j6r7N)}BkZT5~aYri#O$XkZp1AH)Z?)g03qTg!p2&Xx{R!IftT1QrhK4`IO zJ`zM&6G9t1l(F2X&r2IW7;;@$D!j-rPE%QOONgj6xD(>25;@YiX#>`&4Gav8Ho^27 z&;zS6$=R$@cAZe+VE#H5g#_~}jaLf!&Vm1|=mt1U()8h>_nfL4FM*cBIUB-C$tGKWDQ$B`n+@|XqB#*slXj96TkcC2{Fhx*ERt+AuP&993r%fjODE4$XtggvC@xZi)>SWY_r-B^v`O zsUZ?Xay_GnL}yq1G_Lw`vgdwX>XTn=JJ-~ zXFO91ond2ZDswpfTmowsX{*^gi%wU!FOVHd43T*Z>UF1VhV7B5cDJl4(Q)tk1EiK) zdTYVof9jI1+$8mp2MtUZBJSv-a;KLjXIaV7Gza(?GD zmKBr@FdPHj_Q%S2SRCxIWlZON?){#u++?rqwO14!=uMc>DS(1ZHu}JvVL&Se2*?RI zTj2K1W+8gJdm>5>HZX?IBcUEn92XU78^sY|=dE#-18yX9tYjRawZM9^gf8x>sIW@6 zXUUwu5=fxtQ8i1Gy3WUkRIcb*zh%-d6sWCmp4a}STZO)!M0whdv*e-@C4;5hg^5LF zPsJ~BHr=_n`GmgAcK?b#H0xEQhcFpJvkF%izw!Bsg-K2+?fZnQ;ejn!m5tYO0?IyS3NpVIB#h5sLx{0K=G*etdjimMKN+P8OLR%reu>sPudU^)ylS zJ;hY~<$5|qc}ubm+2?L@LgHi5f@Yo*LJUiVz8q$2>$S#uR0whc2y_qG^gb9VyK)2H zE++))uJF@)qlk;G1%jkfN|*z2!xaDAo@aPlm4^xyZ;Ug%2AUg_S^ z5VXnoUcL5fX6kZoot-6z*}8ORK6nexV~Zf&N?3X$RpE5kR7zsT83y+m4Zs=xD!%@% zQZ(1EGYKonc}XyGRZo7atA$wLZs}{SDH7@Xc;>GE+sQ#vkLK@ zBBsd3ilv4lsMnb}y+ofFaBZ!3j80v*wxGZQtZY`OGO^Kvo9i)Rt8f7aAn14r2@ZVY z8!w@)8&o=QK6PU{CT^V87|~@;U|bwrnt^?IC;SCt>%&)Oejq)j}J2XM+Olz_R6Eh&B z6Um1_q!m3mo(_dQh)|Dl<-(O7$43CT<@>J_(IwVmuV8VPKr6FLBG67$wv|jpi$9Wr zhPI$Ov9O-`N3V1${kXk83#jsdp?~pqHKwz-@h_GEu}-%{M**se?P6RRCI#33DDC1> zgotgJnEL$5>6Ex{XEwwj1qPFyOSsRz>Uo)2iEpht)-*XNNS!oJJEq(7RCHoO4RtR2 z%<%Jmb2%!qK8kEj;9T51&NhbfXm?s)dblp?4QU62O3vy`!-7#{o?>kE&Yi4I?-pfi zEJt77r9kG=#rU1(Yj4MmRI?VCHH>|Yv(a}2xx`y~*j5C%)O*PC+ToU3lwRo_Vte=? zi{+#J+BoH|6P!#qTUnIV>Jr;{Zrrg>`!v3p|rBSh;hbI#Oa z+{H01oQ=Jkrh$zgM)Qa;B+T|>;6!s4CpLIBV*d2t+1XM zs=F8S^9)8HSNY7K7r5iGhrLu+&5)ML3hn1C3vO-iX)h4x&|NWiikS(25A)sMWY8&1 z98^jU_$V6^zY@(Lu@P6x$(!1R5WPfgl^wBtVz7xJ4_wRYRk1mU(*QBnx76~E*sc|- zvFImVX=e>_2rc_c8DQNA7G}5ae^fOrO)nQeQIE|Xblene0ZPEQcFzCgD?SUW#Klg& z4J$5REl3v!DK*dMHL>-wXO}2!vF~_*PJzAg4%-<1h~T|dRImbA;~e&-E2NvdNAy5> zjK~+vF<;Sr5}ILIKg#6)ycp~AAchU~Eb1uczrs}u#CN0m-_%vaiSxjEa38#oJ^Vb- z<46ZV=#?}?(2GYauwuQ12hz)dDJNzIH9;#Np&QZGC#%mwutSM~^SgD>5G@+jyd|SG zSKL4dX(;oZ@U>W%xBShPo)I*# z>Tm@F)x*E!K*DrwW6L3a&{}aaa^yN*w|aeSufRGGC#FdR`l_Je0R;#q`dYZ&)@r2$ z8&424#Jwa407YoZ_yx`*opv&5)#X?!+u_s`+mPY;+&~&cK#1k>s99#;Kc@9JD1-jo!Ant{?jMw zg-cxxW*b)(-d!K_aR&Tvkv(}chTs~5K`x=}2GH;V@g56q$9eS(i|n|{sYiKiSB?&7 z#pr7&0k5>7@QS`8`vQ)ksQp&jg!7rvlsgU^~t zozj#C(hWI8)2@tSP~8B%FeZQwYfyg$ygmb9T3(wkZ1(LjJFfix!Zv}s4?u1-n~xK& zjA3;0u2T=ppPjssU>z}-p(}YL9M<_AzS>@UZU4GPoO{=y zbQD(0e{+6O#W@nbAdlTthRaCZC&wujgAYC=&M7qPB?E=Ja_jr6Tf=G&na`noYsFXf zYW7MZ(gs2{nyk?lXx`HWYe6(_8bU9zR;w7VAV7VKw&YIF4|r2qG#*cw$MZVef6%sX zhGMlbpx@S{KidsH*;l;{p@MH-==Vrzh?A3vJCSd2u4Ebyy#AQ9VwF_FxMk5KO292o z)qV@)6`I*^i?^Fxtt8Z6S2IW0_y$&$QF}*+k-P=PBPFz-R)9MsFbCTCNfXg&y(Kte9o$FOUo2pQuSjAlrp z4PHKR&(YfHZ)|G(m;nClGr>+591HdjC~D1}68?b@?Gi$d_tMUAmEJH5MK;{nsgQ{a zo?vN^TUCU_{fnlo`8uQZf~5sTA{|FU{MDL%E3xM$Z6Tk2L=4xrCe0}y!M$Xffu zK3HRitWgNVZ2)L|0_>IaSLjR6%uPKrBwinDHuhd?!K1oMTzav;0^*p@OM*>*1MSKs z)I!jk%x-D))2`Y9%ZKF`7FtRx6DUS~BeE1GAPN9) zIF1MpAHRs|=6cwG>L*quQP7k~sB|hNKD&p;cn@vUCxGM-yHpew=155txNvIRDj`e& z4OS%}bUWusH$)M~ix2C>w;a=hJ(;9m{O5U`t{@*Cz61fmeKLhSVUn$srNAC<_--7J zq6`twN+hy$8o$RyPL>kXWwN@m7`LNOmmXhOO#rjzw{_SL6-t88e4h#5k=|gc3(`$n<0d}o9z<^@O3>qJ(qQ#zHW4r8$NRWZ#NQ8ikFyk&k zK!oJr_z($7z%O*Np7d)={gMu?zO9J_B9?JRK#x;W-Y$ER$dwRU5Z+LX;y`633*{83 z1wCQcQ3)0ypcb4!ztzJS?>VzWUiRrXC6PLA&f6Di{&881kh3KC7Xh#t1;7A=Ex_>P zcF@@eYWCnl%f&CRS1bPGKLkRZQ?%Ey(&nSPQ9ndH{pO8`ya6 zdWy*swBojgcVS5kx%7`2whVE}{_$%=h<6Lns___iY_O)G1DU^(zKgG-+}Q3!OJ<$_ z-cCd*k-)muqD`#gaMe>o->&1y9AZ1Gr;Kn~hxSYYcLw)vuBkb-0v$l0l8^M+9 z(%)RR?U{d?dm)%&nqB2~%0flM38AA`ud=u=9iE^jRuV8ni2urB(IOtv^2ajcu6 zN<61_VEb1ayzR4>y9uK^lGPk?(Dp(dm8XgL8YraJL8_G3sqxYxm=9gXt4a<@OMPpN zGa;3K=ziHbJpSG~M=WFW52&O1#nfsyK*~>|8eb(GHmsN{7r~U}N7Q{RB%r9IXUCS% zKR?y|&8M3u%aP@?HPSam`bn@-k#l2|HCHa1QUc)7c5jBg(N=!QFf%*uPxk1DvoU7h zG7V&(s?uwP(<8LIvA!~8P3Qei1ov0&C|J7!$|ikx+C_PKA#(4sTUyo5)gR@?+2e$V z_qiZTn<}6<;Ij*%5Kc`hlHXa|#!t7DB@E*MPa#4H9olPwhVBRCP&uQp=Fu4NqrUHz z@eNY8`umuG>fslb*|`v?)cSsdn^q1yK{CUlnp;IW2MwvsEtGBiZkYN%S6TW3`!i>7 z5-}eUtyKX*%Ad`n5dRZc7N_Z%V3YF*I1~g~zn@r}x#@&+B zvK!pB@40eU40cj4R+m@WcvzT*K7h6qjW550FwAUfue-V_49j<@LlG;JZK)8ZY%&!t zTTLY-g=d3vS-2?2;qq~W!P~rdR*Z-3$~Y~}zheY3s(=>wIsnB&*xpellEpUxTZ(GF zl?-OJl2qiby=bx!j#LN`4+t5cHZi727nDmhw+=LJrZ!WnLr$`_ z0iR0Bnew_wAR{`m%r=FVk@5wr#st;HvTDz&N9d6d41X$U8fbt7C)>|T1~B%%M}B0@^plQV+uT*i|GO}+_FEOygN4#M@2d~b>P6a(t-M-7kK>Fj{x*9k z$J$~3*iYoxWj_4tqP0n|WlBH%3ZKr8>i~4E?9Z3-&~=s475+vXMF0izyoLent}CST zqKu!q`poA#O+bn&bmB|u$5|uy5j{jEHFfZelew=FmjKtvmweuCMy3iZbCPRUW>;NkIwXS^o36dC zrEedO>;I6KWBIvPt8e06=aQzkYStP~) zkc@G~HD$>fy8P*z5;5-!496vYf^Ky&C)ZD-;iWp2U_4?+4F5#`TD6KSZ^jokC~y2l zz!EGWMpnvMj0h0ix_6@M{za|njr8BG7mNd)x$K1gry zg>@r8x+|%@(tI27+jE!A*#^3XT8Y3&juO1cv_dQEEAA4muR}%Tcn$9~Szmm_I%bG0 z9=IH50-orota4sc_My>uDel(atSM$&Md;BXge&0}|GhG;0Z!yjAG-dM&Dfljj)nnb zgZ`C=m6f)6ZD7El`qT{m&GYcrHa;z&Y!1mPgM1`s7+ySC1S&p{e6aM;w2&1T{jSMX zZL@3EiwolRsi;mHA7NE+Y8tOc{N~#0aWE!z@jAMKNKa+TVbf1R*Fq^@uA4GuxxyYP ziLnGZf6!vBUV|HGR~c+~8htT0XtRj=Fz(*o=T18C-U%-sN#lZ=?2o6|U*lf*eo#6L z9#Ppz|7*isO1_dFqrPSiYH}cvnPmd3`NhPg#>n>j(~wW%hffF9MFz1Gpm$+liwDah z2(XJuCJV(@URCB~H7K?C<&i|e`t+31!AmwJ8!tWD$^p(N$lOp|$Mi7l;7QXnQ*x7Q zqONSZd|r@-`Ik!l)Nhe~uty?c1n0kO=YiNI5?(Mlyo2B=_hCOAJ%;Wbh=~DR^(f545jE49Tlnav| zLM$_ibMHES>lN>OeSdakGj-j162Y++#+JzSu=7FrC*uEi z%NrG4^rN7PC8hd;4(ZD`{JOQK)m1Yaa);~Vn>QT{(zj6XTM(4j-h)2#ww=dRjj*Rr zpr3|XyV)@Enz~O!}3;@kqrlG2EaFB2)gk0kI;k}&Z+?$nJW_{dtl@%k7`f&j) zNzVc{(m$!26vop<%^UNiAE&cHQJ(*IZPH+_5b0Vzq^Ys&xe=zf(OWopG4uIMK)*Ci zr-xlO?(8K9(CIuLiIKWAEHk;ytLv+%V>{pS!x9}^ELCGAM zVE&D7ASqnQHnoO2OX;j&AAdudEq z`PAIa(sDD_bAJ-ZGbG}Rb*JS2kD(i^2}Yev;Zt1{hd(rp4h`{C`~A^L0zU|Zc^+E9@ zt^_44Nan3C{JEXu=8m+41i6egArL7NpgI!7HY~B+^FGltU&Jd(y-&#ySr91=-zj`2 zIs3%j{8Um4xH~-(5w^uwivDBY1s2RFC!_4y9g6BTWvk-_{3g#YUT1%sp3}88Yh?v< znii5w-`1-KxIEj(0L;)6yiooX?tnB9x{W>ldH+2Lk@|Fi8C86{-LR{ zk+iZXc5Q{Pp}4R$#gwlmJb%bP1~&We;I4o;9MR@2@c5=8@B4YQ%jJXlyZsx;9$!v4UsY@7!*@9Z>&A3FvymBH=R=(vNbH#PcYrJ8owc* zTuk%WF-%v{W{bl4;hAO<48e)`%`$j-suiS?8j!6hfjoIEK{i=U?PHzB5}M@z7eTKx z`di(sQ+gWem7$?+)$b|hdTJ6HJGsZ7cry& zSO+|_)wmIGnDVoHX8K4rYQIX0_KEI~Az>T!41-Uhq{Mk37A%{%ZS-4_{;vA${pE_91pKB{QVEl6U z(Qgf>Yfh1sFj9R_<%>PS8%O62W?6dsowxa%<;Ib;z3wTG%Mdkv`Hs)q0*{`ZsG-`7X>#_lwgsBnbSX0>``;pNTTvGIOt&7R4wd>Ja-> zu;_R-=sH>2DmZ^6hQhsyyxdx$RLkV8#LzzOvMKE=ev(b-Ji)kL`Ozht8~Op81F=b) z9h_Ql(Zlz<+uai#YvfHw&7Ry$QfOL~TvHKA?ZH6s*JlpR74|JEyiMWi_9xM)=`vNUMAvW9l%iq0t2?-A9GJ03pPx{0XE)ZzAkCKw%cudk1)Hzyn9eJwcviLc z+OAp3TbHw8Z~VPwX?U?`e=h2U(Y6*BN5GGmV2_Snwk@04zw!@tAU?YEdL4*8;JjzL z4qb;2$5ZTn*Chow(R?vSh?x^%f6n~QzO@Wp1rZzP>*r=47f?-RJ3O`zH&X;PxdiUl zPnnCYuMhNY${fK3IEM0*2mYbeL-etoB%1p!M{IO4q>ro;C{29`6peCZDuR-ixFcCg zx54MV$>96~98qK*O?5oXHxG+Z%i0xoLc51UKx-+WO5$q0pHycFk`gW}xkuZ>TpJ@~ zY<_%cuhe!=NU(FqL+QMyThnTh8$FtD1npuAK*$2C*Z2NB$gu7q6EkIz?V`fE2ZGv)m^rFQ-ivli|>3aeLDYc@L-KKeQu3G!udqeJ52ztQ^k~^(W z(vwmbtp~?$94v}<DePrSPFE=*?m184BkznUcrs7>t}y$pkdJJ} zz;T>cX(5PVsWL>OoWM^e|K?OIW+D+S)VikI4_@mnY~uQ;n5z!O8GN05@6~JVC?bCh zJxZ{~j8kZRsB0BLDH=i zyT^7#Ta_3%3m(thXOCrEZQ%2=jL)nM(DO} zXGN-pk0!JPNbL~ZsLY%zpALAEQjit{e=nD2P7mdGTEA_OG$9pdO)ijD`FKK#du!1c zVtCqJp-eZ3(~0C5gn8y2l%v%dO&(6CKlKzUYxeStnCOcXNQD{OnsrQU&ijH!@ z;@*g!k>W|7NDOtPjy)x;(}Al4Z9r5X!$BCLL7r;iM}+h-&qdYIaYfkPpZd&hSn=Ye z@q0u%jTjYi8!Ng^(z#XnB0BPFh{&eSV~uTfE&=S-C`0vQCvXv8^gCV6tMYIn{k1MJ z{B;sltrZ4VTH-n$AiQMla&7T$uQXu!MS9dYY5DO9Q0q)^;gUO2 zj03LJ8Vqwo6h!0h0xtKrl!^rawlGo7FyMQ zWw{F+zJh1;vdtK_#r+5?{&?|%9}Z`Tgg5s<5H@5r_v=QDeUg_p3k!-SA-dcMkZ!fJ zzBxx_y^u((7q?dLt$jN)s$nH0_Zc$r+g=h!F_taPt4p#hT?F|uv_%{7#DzdY)+5Vb?LU0ex_g@6=)No zE3_$>nX#}|PCRc*uL#{uJ$v5+ZwRfMrmDw>TNDme^bWq+=fG>iV>!EUBpJ?tt0_(? z1VrkLTV>Q@P~4=lSYckh`WCUpT8YqN#$>pyuPt|5n&4d!=~miTxVb#j@yQO48jx6c zTWeH8Oo((=4FdiUWkkIp8clVt<&ts^Q3lsvQN7=DaI^Bxkn?h5gco9Hl-7r_&?|po z2yi!^TK)y*KL%!h{(kh(uZZ!r{gmA^9M&GxCh{g_2QGjaY>HrK;Jw#)XtkH8Hc1rRf@|BK zx+OfejoX{J{1nvl!m?&yWE&?m~|Nz*m69v6|c{*zXl!90s~^N6F~{;E>RSs^+O z_zh44aWc8?G7=WKnm$lIyn1vTg9gvGmTtSGs z!~1;w#5P#S_>`=~8o0)Xj)g6Edca3*QUuB**u?!3(`^1s#hQ;2XV&5z>dtAp_!>|1 z59N?#dgDu>m|loc`=-+O`=o=(rtzbaJbu9>=qkAZpMr0K6yE*$$#-kT$v)FipTe=_ zFRuc|r3UdP5aOZs4U+W#h2*i|o8x>HQ7I`faqQuJHm*-B)xMWn7L<>CIySp$&wyVQ zKKAvL#v=t}ooU}Q)w}dA(!{dlArg)eH5t_6m0}x42uqL-rdLQpo|YxRW&Qu@qUG>A zW92wrI8E_W&e;$0ACaFnEu<-qC+y*g@6>{EFc^`g02Q?##JK(#Dq)w+<%BA!Z=a@z Wui6EVap>3o;_gVN3K9!3fNnqorL9~5 literal 0 HcmV?d00001 diff --git a/man/OmicsIntegrator-package.Rd b/man/OmicsIntegrator-package.Rd new file mode 100644 index 0000000..586bb33 --- /dev/null +++ b/man/OmicsIntegrator-package.Rd @@ -0,0 +1,15 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/OmicsIntegrator-package.R +\docType{package} +\name{OmicsIntegrator-package} +\alias{OmicsIntegrator} +\alias{OmicsIntegrator-package} +\title{OmicsIntegrator: Calculate correlation between two data frame} +\description{ +This function takes in two data frames df1 and df2, and calculates the correlation matrix between the columns of both data frames. The user can specify the correlation method to be used by passing the method parameter. By default, the function uses Pearson correlation. The function calculates the p-values for the correlation tests and adjusts them for multiple testing using the specified adjustment method (default method is False Discovery Rate adjustment). The function returns a list of the correlation results including the correlation matrix, the p-value matrix, the adjusted p-value matrix, and the significance matrix indicating the significance level of the correlations. +} +\author{ +\strong{Maintainer}: Shashank Gupta \email{bioinformaticssrm2011@gmail.com} (\href{https://orcid.org/0000-0001-8839-7003}{ORCID}) + +} +\keyword{internal} diff --git a/man/calculate_correlations.Rd b/man/calculate_correlations.Rd new file mode 100644 index 0000000..29f169a --- /dev/null +++ b/man/calculate_correlations.Rd @@ -0,0 +1,39 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/calculate_correlations.R +\name{calculate_correlations} +\alias{calculate_correlations} +\title{Calculate correlation between two data frames} +\usage{ +calculate_correlations(df1, df2, method = "pearson", adjust_method = "fdr") +} +\arguments{ +\item{df1}{A data frame containing the first set of variables to be correlated. The samples should be in the rows and the MEs in the columns.} + +\item{df2}{A data frame containing the second set of variables to be correlated. The samples should be in the rows and the MEs in the columns.} + +\item{method}{(optional) The correlation method to be used. Default is "pearson".} + +\item{adjust_method}{(optional) The method for adjusting the p-values for multiple testing. Default is "fdr".} +} +\value{ +A list containing the correlation results + +$correlation Matrix of correlation coefficients. + +$p_value Matrix of p-values for the correlation coefficients. + +$p_value_adj Matrix of adjusted p-values for the correlation coefficients. + +$signif_matrix Matrix indicating the significance level of the correlations. +} +\description{ +This function calculates the correlation matrix between the columns of two data frames. +} +\examples{ +# Calculate correlation between two data frames using default parameters +result_list <- calculate_correlations(df1 = transcriptomics_data, df2 = metagenomics_data) + +# Calculate correlation between two data frames using Kendall correlation and Bonferroni adjustment +result_list <- calculate_correlations(df1, df2, method = "kendall", adjust_method = "bonferroni") + +} diff --git a/tests/testthat.R b/tests/testthat.R new file mode 100644 index 0000000..2b656fe --- /dev/null +++ b/tests/testthat.R @@ -0,0 +1,12 @@ +# This file is part of the standard setup for testthat. +# It is recommended that you do not modify it. +# +# Where should you do additional test configuration? +# Learn more about the roles of various files in: +# * https://r-pkgs.org/tests.html +# * https://testthat.r-lib.org/reference/test_package.html#special-files + +library(testthat) +library(OmicsIntegrator) + +test_check("OmicsIntegrator") diff --git a/tests/testthat/test-OmicsIntegrator.R b/tests/testthat/test-OmicsIntegrator.R new file mode 100644 index 0000000..8849056 --- /dev/null +++ b/tests/testthat/test-OmicsIntegrator.R @@ -0,0 +1,3 @@ +test_that("multiplication works", { + expect_equal(2 * 2, 4) +}) diff --git a/vignettes/.gitignore b/vignettes/.gitignore new file mode 100644 index 0000000..097b241 --- /dev/null +++ b/vignettes/.gitignore @@ -0,0 +1,2 @@ +*.html +*.R diff --git a/vignettes/OmicsIntegrator.Rmd b/vignettes/OmicsIntegrator.Rmd new file mode 100644 index 0000000..40f818d --- /dev/null +++ b/vignettes/OmicsIntegrator.Rmd @@ -0,0 +1,19 @@ +--- +title: "OmicsIntegrator" +output: rmarkdown::html_vignette +vignette: > + %\VignetteIndexEntry{OmicsIntegrator} + %\VignetteEngine{knitr::rmarkdown} + %\VignetteEncoding{UTF-8} +--- + +```{r, include = FALSE} +knitr::opts_chunk$set( + collapse = TRUE, + comment = "#>" +) +``` + +```{r setup} +library(OmicsIntegrator) +```