Skip to content

Commit ef92299

Browse files
committed
editing index file and pushing new visualization images
1 parent 3b7374c commit ef92299

8 files changed

Lines changed: 56 additions & 11 deletions

PROMPT_ACTION_LOG.md

Lines changed: 42 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -377,3 +377,45 @@ Updated `workflows/code/build_graph_inputs.R` to use the user's original SpiecEa
377377
Action
378378

379379
Updated `.gitignore` to ignore generated workflow outputs under `workflows/output/**` to avoid pushing large local artifacts to GitHub. Added an exception for `workflows/output/.gitkeep` so the output directory can still be tracked when needed.
380+
381+
2026-05-14
382+
383+
Action
384+
385+
Recorded the active project goal in `workflows/README.md`: build a graph neural network workflow to predict community coalescence from donor and resident species co-occurrence/community matrices before communities are joined and compare predictions to observed final communities. Also cleaned up `workflows/README.md` formatting and expanded `workflows/code/README.md` with the current graph-input workflow and rerun command.
386+
387+
2026-05-14
388+
389+
Action
390+
391+
Added `workflows/code/train_coalescence_gnn.R`, an R `torch` workflow that reads `workflows/output/graph_inputs/` files, builds taxon-node prediction examples from `coalescence_triplets.csv` and `combined_sample_taxon_edges.csv`, trains separate bacteria/fungi multi-task models to predict final presence and abundance from donor/resident features, and writes model outputs under `workflows/output/gnn_model/`. The first-pass model intentionally ignores the empty SpiecEasi edge file and uses self-loop graph blocks plus learned taxon embeddings so it can run locally or on HPC before association edges are available.
392+
393+
2026-05-14
394+
395+
Action
396+
397+
Renamed the GNN training script to `workflows/code/train_coalescence_gnn_CMKR.R` so CMKR-created code files are clearly labeled. Updated the model to read `taxon_taxon_spieceasi_edges.csv` as a signed, weighted taxon-taxon association scaffold, build positive and negative normalized adjacency matrices, use those matrices in separate message-passing channels, and write `spieceasi_graph_summary.csv` plus `spieceasi_edges_used.csv` with model outputs.
398+
399+
2026-05-14
400+
401+
Action
402+
403+
Added explanatory comments to `workflows/code/train_coalescence_gnn_CMKR.R` without changing the executable code. Comments now describe package checks, configuration, input loading, taxon feature construction, SpiecEasi adjacency construction, data splitting, model architecture, training losses, evaluation metrics, and output files. Confirmed the script still parses with Rscript.
404+
405+
2026-05-14
406+
407+
Action
408+
409+
Added `workflows/code/visualize_gnn_performance_CMKR.R`, an R script that reads model outputs from `workflows/output/gnn_model`, visualizes kingdom-specific and combined GNN performance, writes plots and summary tables under `workflows/output/gnn_model/performance_visualizations_CMKR/`, and documents the generated visualization files. Confirmed the script parses with Rscript; local execution was not run because `data.table` and `ggplot2` are not installed in this Windows R environment.
410+
411+
2026-05-14
412+
413+
Action
414+
415+
Revised `workflows/code/visualize_gnn_performance_CMKR.R` to be fully separate from model training/compilation and fully tidyverse-based. The script now explicitly loads saved performance CSVs from `workflows/output/gnn_model`, uses `readr`, `dplyr`, `tidyr`, `purrr`, and `ggplot2` workflows rather than `data.table`, and still writes visualizations and tidy summaries under `workflows/output/gnn_model/performance_visualizations_CMKR/`. Confirmed the script parses with Rscript; local execution was not run because `tidyverse` is not installed in this Windows R environment.
416+
417+
2026-05-14
418+
419+
Action
420+
421+
Copied `workflows/code/train_coalescence_gnn_CMKR.R` to `workflows/code/train_coalescence_gnn_presence_absence_CMKR.R` for a presence/absence modeling variant. The copied script now binarizes all observed abundance values before feature and target construction, disables abundance loss by default, writes outputs to `workflows/output/gnn_model_presence_absence`, and replaces Bray-Curtis dissimilarity metrics with presence-set Jaccard and Sorensen similarity metrics for predicted versus observed final communities. Confirmed the copied script parses with Rscript.

docs/index.md

Lines changed: 14 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -186,37 +186,40 @@ People who could use this:
186186

187187
## Findings at a glance 📣 { #findings-at-a-glance .oasis-report-out-section .oasis-report-out-day3 }
188188

189-
Headline 1 — what, where, how much
189+
Headline 1 — This type of data is adequate to successfully compile a GNN.
190190

191-
...
192191

193-
Headline 2 — change/trend/contrast
192+
Headline 2 — The GNN performs better on presence/absence data than on abundance data, but the predictive capacity on our initial dataset is still poor.
194193

195-
...
196194

197-
Headline 3 — implication for practice or policy
195+
Headline 3 — Looking at a different dataset (See fecal transplants above) will both broaden applicability into the medical field and might provide cleaner data to train the GNN on.
198196

199-
...
200197

201198
## Visuals that tell a story 📣 { #story-visuals .oasis-report-out-section .oasis-report-out-day3 }
202199

203-
![Story visual](assets/figures/main_result.png)
200+
![Second visual](docs/Train_vs_validation_loss_MS.png)
201+
202+
*Visual 1: Model loss showing convergence*
203+
204+
205+
![Story visual](workflows/output/gnn_model/performance_visualizations_CMKR/validation_metrics_by_epoch.png)
206+
207+
*Visual 2: The GNN performs poorly on community abundance data, but still performs slightly better than random chance on presence/absence.*
204208

205-
*Visual 1: the main pattern or output we want people to remember.*
206209

207210
## What’s next? 📣 { #whats-next .oasis-report-out-section .oasis-report-out-day3 }
208211

209212
Short term:
210213

211-
- ...
214+
- Prepare a cleaner dataset to determine if the issue is with the data or the model architecture.
212215

213216
Long term:
214217

215-
- ...
218+
- Finalize a model architecture and test it with multiple types of data including medical, restoration, and agricultural community coalescence experiments.
216219

217220
Who should see this next
218221

219-
- ...
222+
- Someone who understands GNNs better than we do.
220223

221224
## Cite & Reuse { #cite-reuse }
222225

143 KB
Loading
143 KB
Loading
134 KB
Loading
117 KB
Loading
111 KB
Loading
253 KB
Loading

0 commit comments

Comments
 (0)