Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
BDS
GeneDER
GeneDER_core
Commits
5be415ac
Commit
5be415ac
authored
Nov 11, 2021
by
Leon-Charles Tranchevent
Browse files
Corrected garbage collection.
parent
f293ccca
Changes
3
Hide whitespace changes
Inline
Side-by-side
06-Data_integration/compute_gender_specificity.R
View file @
5be415ac
...
...
@@ -303,11 +303,10 @@ for (i in seq_len(length(config$integrations))) {
# We compute the delta in both cases.
FM_enriched
<-
compute_genderspecificity_scores
(
FM
)
MF_enriched
<-
compute_genderspecificity_scores
(
MF
)
##FM_enriched_FDR <- compute_genderspecificity_scores_FDR(FM)
##MF_enriched_FDR <- compute_genderspecificity_scores_FDR(MF)
rm
(
FM
,
MF
)
FM_enriched_FDR
<-
compute_genderspecificity_scores_FDR
(
FM
)
MF_enriched_FDR
<-
compute_genderspecificity_scores_FDR
(
MF
)
# We save the full data.
FM_ofile
<-
paste0
(
output_data_dir
,
integration
$
name
,
"_"
,
limmas
[[
5
]]
$
name
,
"_max-avg_genderspecificityscore_fulldata.tsv"
)
...
...
06-Data_integration/merge_and_filter_rankings.R
View file @
5be415ac
...
...
@@ -151,7 +151,7 @@ for (i in seq_len(length(config$integrations))) {
# B: saving the tables (for manual exploration of the results).
#################################################################################################
# We first add the pi values to the resuts of the male and female analyses.
# We first add the pi values to the resu
l
ts of the male and female analyses.
res_females
$
pi_value
<-
compute_pi_values
(
res_females
$
log_fold_change
,
res_females
$
P_value
)
res_males
$
pi_value
<-
compute_pi_values
(
res_males
$
log_fold_change
,
...
...
06-Data_integration/refine_GD_rankings.R
View file @
5be415ac
...
...
@@ -17,8 +17,6 @@ options(bitmapType = "cairo")
config
<-
read_config
(
config_dirs
=
c
(
"../Confs/"
,
"./"
))
output_data_dir
<-
paste0
(
config
$
global_data_dir
,
config
$
local_data_dir
)
input_data_dir
<-
paste0
(
config
$
global_data_dir
,
config
$
local_input_data_dir
)
input_edata_dir
<-
paste0
(
config
$
global_data_dir
,
config
$
local_exp_data_dir
)
input_rsedata_dir
<-
paste0
(
config
$
global_data_dir
,
config
$
local_rs_exp_data_dir
)
message
(
paste0
(
"["
,
Sys.time
(),
"] Configuration done."
))
# ================================================================================================
...
...
@@ -93,5 +91,5 @@ for (i in seq_len(length(config$integrations))) {
# We clean the workspace and log the session details for reproducibility.
rm
(
i
)
rm
(
config
,
output_data_dir
,
input_data_dir
,
input_edata_dir
,
input_rsedata_dir
)
rm
(
config
,
output_data_dir
,
input_data_dir
)
sessionInfo
()
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment