Skip to content

data: deposit canonical identity-cluster artifacts (recovered from permea-h100) - #3

Open
Permea-lab-admin wants to merge 1 commit into
mainfrom
fix/deposit-canonical-cluster-artifacts
Open

Permea-lab-admin wants to merge 1 commit into
mainfrom
fix/deposit-canonical-cluster-artifacts

Conversation

@Permea-lab-admin

Copy link
Copy Markdown
Contributor

What

Deposits the eight identity-cluster assignment TSVs used to produce every Condition-B number in the paper, with checksums and provenance.

They were never committed and were not in the v1.0.1 Zenodo deposit. They existed in exactly one place — a working directory on permea-h100. This PR gets them into the archive.

Why this is urgent

Regeneration does not substitute for the artifact. MMseqs2 18-8cc5c on the same input at τ=0.6 produces 2918 clusters, not the 2913 used for the paper — a different partition, and therefore different Condition-B numbers. The MMseqs2 version used at generation time was never recorded.

A checksum of a file nobody kept proves nothing. The files themselves are now archived.

Contents

File Backend τ Clusters
clusters_thr0.6.tsv MMseqs2 0.6 2913 ← headline (ROC-AUC 0.853 / PR-AUC 0.539 / MCC 0.479)
clusters_thr0.{4,5,7}.tsv MMseqs2 0.4/0.5/0.7 2907 / 2908 / 2919
clusters_align_thr0.{4,5,6,7}.tsv Biopython global align 0.4–0.7 2752 / 2823 / 2863 / 2887

Plus MANIFEST.sha256 and a README.md recording the identity definition (matches / shorter sequence, global, gaps free), the MMseqs2 flags, and the provenance gap.

Each file covers exactly the 2959 sequence_ids of the canonical dataset — verified.

Environment for bit-exact reproduction

Now recorded, because it was not before:

platform : x86_64 Linux
python   : 3.10.12
sklearn  : 1.7.2
numpy    : 2.2.6
scipy    : 1.15.3

The published numbers are bit-exact only on that platform. On arm64 (Apple Silicon), ROC-AUC reproduces exactly, but PR-AUC and MCC deviate by ~0.001: floating-point ordering differs, and the F1-maximising threshold search then lands on a different grid point. A follow-up PR will correct the README's "reproduce every number on a laptop" claim, which is not true on arm64.

Scope

Additive only. No existing file is modified; no result changes. The harness reads these from results/ exactly as before.

🤖 Generated with Claude Code

The eight cluster assignment TSVs used to produce every Condition-B number in
the paper were never committed and were not part of the v1.0.1 Zenodo deposit.
They existed in exactly one place: a working directory on permea-h100. This
deposits them, with checksums and provenance.

Why this matters: regeneration does not substitute for the artifact. MMseqs2
18-8cc5c on the same input at tau=0.6 yields 2918 clusters, not the 2913 used
for the paper, and therefore different Condition-B numbers. The MMseqs2 version
used at generation time was not recorded. A checksum of a file nobody kept
proves nothing, so the files themselves are now archived.

Also records the exact environment for bit-exact reproduction (x86_64,
python 3.10.12, sklearn 1.7.2, numpy 2.2.6, scipy 1.15.3). The published
numbers are bit-exact only on that platform; on arm64 ROC-AUC reproduces
exactly but PR-AUC and MCC deviate by ~0.001 through float ordering and the
F1-max threshold grid.

Adds:
  results/clusters/*.tsv          8 artifacts, 2959 rows each
  results/clusters/MANIFEST.sha256
  results/clusters/README.md      provenance, identity definition, environment

Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants