plot_ttd: Addfill=TRUEwhen binding dts
run_permutation_test_cor: New function to run permutation tests.plot_ontology_levels: Add more args. Make "info_content_binned" a separate col from raw "info_content".
- Update all calls to
facet_gridto not use deprecatedfacetsarg.
- Copy all core
phenomixfunctions used inMSTExplorerinto new packagephenomix.coreto reduce dependencies and issues installingphenomix(which requiresscKirbyandechoconda). validate_associations_correlate_ctd: Drastically reduce figure size by aggregating
data to ancestral phenotype groups first.
plot_celltype_severity: Ensure cell types are ordered by composite severity.get_ggstatsplot_stats: Return all stats, not just specific labels.ttd_check\ttd_plot- Make % prioritised targets a percentage of the total number of rows.
- Fix mistake where I was including non-gene therapies too (due to how grepping collapsed listed works in data.tables).
- Account for situations when there's no failed targets to check.
plot_bar_dendro: Pass "top_ancestor_name" toget_color_map=arg instead of "ancestor_name".
validate_associations_correlate_ctd: Adddownsamplearg to reduce plot size.
get_cl: new wrapper for getting the exact version of the Cell Ontology used in our original study.- Update README to better reflect the package's current state.
- Use
BiocManangerto install in README. plot_bar_dendro: Fix #22plot_tissues: Don't try to facet plot types that aren't requested.prioritise_targets_network: Fix examples and tests.plot_celltype_severity: Fix tests.
get_color_map- Fix
color_vectorassignment.
- Fix
-
ttd_hypergeoviattd_check: Missing column (NULL) made overlap always 0, and thus hypergeo p-value always equal 1. -
Tests
- Use
force_new=TRUEwhere some tests occasionally fail with cached files. test-load_example_results: Update test files.test-prioritise_targets: Remove unused arguments and change input size.test-prioritise_targets_network: Processtop_targetsto include effecttest-plot_differential_outcomes: Use non-specific plot name inpatchwork::wrap_plots.test-plot_differential_outcomes: Wrap p3 inexpect_errorto prevent test failure even if error was handled.test-report_plot: Fix detection for ggplot object. variable.
- Use
-
Vignettes
MSTExplorer: Update effect variable tofold_change.
-
add_logfc- Return
resultswith new column rather than directly modifying the original input. - Update references (
add_logfc(results)->results <- add_logfc(results)).
- Return
-
ttd_check,plot_differential_outcomes,prioritise_targets_grid- Add check for
disease_namecolumn before executingHPOExplorer::add_diseaseon input.
- Add check for
-
plot_ttd- Remove
fillaesthetic forgeom_text(doesn't exist anymore).
- Remove
-
extract_help- [DEVELOPEMENT ONLY] Look for help docs only in legitimate pkg installation paths.
-
subset_results- Add new
effect_varargument. - Adjust default
effect_thresholdto 0.1.
- Add new
-
add_symptom_results- Only merge
resultsandphenotypes_to_genesif required (prevents column duplicates with altered names).
- Only merge
-
map_tissue- Rewrite merge logic to fix error: attempt to replicate non-vector.
-
test_target_celltypes- Uncomment
add_ancestors(results)to ensure required ancestor columns are present.
- Uncomment
-
Add missing import: simona
run_phenomix: can now take multipleannotLevels at once (will iterate).
- add function
run_congenital_enrichment - add function
plot_celltype_severity
- add function
prioritise_targets_grid
plot_bar_dendro_facets- Add
hpoarg to pass HPO directly. - Fix
data_summaryso that denominator is the number of on-target cell types only.
- Add
- Rename
MSTExplorer. - Revamp to work with
KGExplorer.
gene_results- Return merged versions of
$resultsas well as$gene_data(when available).
- Return merged versions of
- New function:
get_bg:- Creates and caches background made with
gprofiler.
- Creates and caches background made with
- New function:
standardise_genes- May move this to
orthogenepackage as it's quite generally useful.
- May move this to
- Allow users to set min number of genes:
min_genesarg ingen_results/ewce_para.
gen_results/ewce_parabgwas incorrectly set to use only genes ingene_data.- Now uses
get_bgto create background using gprofiler.
get_valid_gene_lists- Throw error when 0 valid gene lists found.
get_unfinished_list_names- Fix example
- Use rworkflows@dev to avoid API limit.
MultiEWCEfinally gets a hex sticker!gen_results/gen_overlap- Check for existing results and import if already there.
- Name all results "gen_results.rds" or "gen_overlap.rds" to avoid rerunning duplicate analyses on HPC.
- Update rworkflows.yml
- Fix unit tests and examples to use "hpo_id" instead of "hpo_name".
load_hpo_graph:- export
- Regenerate and update "hpo_graph.rds" file.
- Drastically reduce time to run examples.
ontology_plot- Fix function and add test.
- DESCRIPTION
Depends: R (>= 2.10)-->Depends: R (>= 4.1), to ensure|>function available.- Rewrite
Descriptionfield to reflectMultiEWCE's current purpose.
- Update to coordinate with
HPOExplorerupdates. - New funcs:
plot_ontology_levelsttd_checkttd_plotttd_import
- Require
HPOExplorer (>= 0.99.10) - Switch terms:
- "HPO_ID" --> "hpo_id"
- "Phenotype" --> "hpo_name"
- "FREQUENCY" --> "gene_freq"
- "Onset" --> "onset"
- "Modifier" --> "modifier"
- "Aspect" --> "aspect"
- "Gene" --> "gene_symbol"
- "DatabaseID" --> "disease_id"
- "LinkID" --> "disease_id"
- Update all "data" objects.
get_data- Add
tagarg.
- Add
load_example_results- Change
tagto "latest". - Update colnames dynamically.
- Change
load_example_ctd- Change
tagto "latest".
- Change
load_hpo_graph- Change
tagto "latest".
- Change
map_tissues- Fix docs.
agg_results- Add "hpo_id" column before aggregation.
correlation_heatmap- Ensure row annot order is correct.
add_ctd- Fix example.
prioritise_targets- New arg
evidence_score_thresholdto utilise GenCC evidence scores provided viaHPOExplorer. - Merge symptom-level driver genes after unlisting the "intersection" column instead of using indirect approach that on nested "intersection" column.
- New arg
prioritise_targets_network- Add feature to resize plot after double-click.
- New function:
add_tissues- Add tissues that each cell type is found in, using a celltype-tissue mapping file.
- Regenerate "DescartesHuman_celltype_mapping.csv" file as the one previously generated wrong (had every celltype in every tissue)
prioritise_targets_network- Fill screen better in saved plots with new default args:
width = "100%", height = "90vh"
- Fill screen better in saved plots with new default args:
prioritise_targets- Move filtering steps and arg docs inside respective
HPOExplorer::add_*functions.
- Move filtering steps and arg docs inside respective
gen_results- New arg
parallel_bootlets users choose which level to parallelise over.
- New arg
- New exported func:
gen_overlap- Compute simple overlap enrichment results.
- Overcomes requirement of >=4 genes.
- Faster than EWCE (but less robust)
- New internal func:
save_results load_example_results- New arg
force_new
- New arg
get_valid_gene_lists- Can now handle any
list_name_column
- Can now handle any
-
ggnetwork_plot_full:- Handles multiple celltypes by aggregation.
-
New internal function:
agg_results. -
New functions for generating/visualizing prioritised targets:
prioritise_targetsprioritise_targets_heatmapprioritise_targets_networkplot_reportcorrelation_heatmapfrequency_barplot_frequency_histogram
-
Added
example_targetsdata:- Speeds up run time of examples/tests.
- Remove unused function:
count_results - Make
get_unfinished_list_namesmuch more efficient. - Add leeway to
gen_resultstest, as the number of sig results varies from run to run.
- New functions:
load_example_resultsontology_plotterminal_celltypes- various support functions
- Add
ontologyPlotas new Import. - Was able to replicate Momoko's results!!!
load_example_ctd- Pass
filetopiggyback::pb_download
- Pass
- Remove redundant
get_gene_listfunction (now handled byHPOExplorer::get_gene_lists). - Make
get_valid_gene_listsmuch more efficient and consider the intersect between ctd/gene_data genes. merge_results:- Handle
NULLresults in list.
- Handle
- Added a
NEWS.mdfile to track changes to the package. - Added
rworkflows ewce_para- Harmonize arguments with
EWCE - Set default
list_namesarg - Output named list of saved files instead of
TRUE. - Add messages that appear in parallel. Make EWCE message silent.
- Harmonize arguments with
ewce_plot:- Simply make this a shallow wrapper for
EWCE::ewce_plotas the latter function has since been fixed.
- Simply make this a shallow wrapper for
- Speed up with
lapplythroughout. gen_results:- Allow speedup by skipping writing to RDS and returning merged results directly instead.
- Use
@inheritParamsthroughout. - Use
\linkthroughout. load_example_ctd:- Cache file
merge_results:- Now handles lists of file path AND list of results lists.
- Remove
RDA_assign_load- Not used and supplanted by
EWCE::load_rdata()anyway.
- Not used and supplanted by
- Remove unnecessary parameter descriptions in MultiEWCE vignette.
- Change
results_dir-->save_dirthroughout. load_example_CTD-->load_example_ctdfor function naming consistency.- Add README.Rmd.
- Update to use
EWCE>1.0.0 - Make more functions internal:
is_not_analysedget_valid_gene_lists
- Make paths with
file.pathinstead of hard-coded "/" - Write files to
tempdir()instead of current working directory. - Remove unused packages:
cowplotewceData
- Fix all unit tests (not set up in correct format originally).