Full-Likelihood Auxiliary Representation Training
Implementation & benchmark of arXiv:2608.26610 — On efficiency gains via augmenting a tiny sample with a massive auxiliary sample
Source paper: On efficiency gains via augmenting a tiny sample with a massive auxiliary sample arXiv:2608.26610 ⓘ · analyzed Aug 29, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Train a neural model on a tiny target dataset and a massive auxiliary dataset using a joint full-likelihood objective: fit the target predictor on target observations while modeling whether each observation came from the auxiliary or target population. Share the encoder between the target model and odds head so that auxiliary examples improve parameters identifiable from both distributions, rather than merely reweighting the few target examples. Freeze or strongly regularize target-specific prediction heads to reduce negative transfer.
Formulas
Mathematical statement
Tukey's factorization uses the conditional odds O(x)=P(A=1|x)/P(A=0|x), where A=0 denotes the target population and A=1 the auxiliary population. If p_0(x)=p(x|A=0), p_1(x)=p(x|A=1), and pi=P(A=1), Bayes' rule gives O(x)=pi p_1(x)/((1-pi)p_0(x)), hence p_1(x)/p_0(x)=((1-pi)/pi)O(x). A neural odds model q_phi(x)=P(A=1|x) has O_phi(x)=q_phi(x)/(1-q_phi(x)). The paper's full-likelihood strategy jointly estimates target-distribution parameters and odds parameters instead of treating the odds as fixed inverse-probability weights. Its Fisher-information identity is I(theta)=E[Var(U_data|A)]+Var(E[U_data|A]), where U_data is the data-dependent score and theta denotes shared parameters. The second term is between-population information: it is nonzero when target and auxiliary populations have different feature statistics, and can be estimated more accurately from the massive auxiliary sample. Adapt this with shared encoder z=h_theta(x), target predictor p_psi(y|z), and domain model q_phi(x)=sigmoid(g_phi(z)). The joint objective combines target task likelihood and domain likelihood, with a tunable coefficient lambda.
Implementation notes
(1) Exact integration point: add a domain-odds head to the encoder used by the target network. For each input x, compute z=h_theta(x), target logits f_psi(z), and domain logit r_phi(z); q_phi(x)=sigmoid(r_phi(z)) estimates P(A=1|x). Apply the target loss only to the n0 target examples. Apply binary domain cross-entropy to every target and auxiliary example. Do not use auxiliary labels in the target prediction loss unless they are known to be semantically compatible. (2) Pseudocode: initialize shared encoder theta, target head psi, and odds head phi. At each step draw minibatch B0 from the target data and B1 from the auxiliary data, form B=B0 union B1, compute z=h_theta(B), and set L_target=mean over B0 of negative log p_psi(y|z). Compute q=sigmoid(g_phi(z)) and L_domain=mean over B of [-a log q-(1-a)log(1-q)]. Optimize L=L_target+lambda L_domain+gamma||theta-theta_anchor||^2 with AdamW. The odds estimate is O_hat=q/(1-q), or log O_hat=r_phi; clip r_phi to [-c,c] when using it for diagnostics or optional importance weighting. (3) Mathematics versus estimation: the odds identity and Fisher decomposition are population identities; q_phi, encoder gradients, and covariance terms are empirical estimates. Estimate the information diagnostic from per-example gradients u_i=grad_theta log p or grad_theta log q, reporting within-domain covariance and covariance of domain means. (4) First cheap experiment: use CIFAR-10 or a small domain-shift benchmark, retain 1%, 5%, and 10% of target examples, and create a large auxiliary pool from controlled corruptions, styles, or related classes. Compare target-only fine-tuning, ordinary domain-adversarial training, inverse-density weighting, and this joint objective at equal target updates and equal total FLOPs. Success means higher target accuracy or lower target NLL at each target size, with gains increasing as n1/n0 grows. Failure is degradation versus target-only training or rising domain accuracy accompanied by falling target accuracy. Ablate encoder sharing, lambda, odds clipping, and the Fisher-inspired between-domain-gradient diagnostic.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Not tested
Stage 2 · Practical benchmark result: Not run
Methodology: Toy-system gate first; the benchmark stage runs only after a pass. How verification works
Stage 1 — Mechanism check agent confidence 7/10
Built a shared-encoder target classifier plus auxiliary/target domain-odds head trained with the joint full-likelihood objective. The odds identity and Fisher total-variance decomposition passed numerically, with relative error 2e-14 and decomposition error 0. On the tiny synthetic shift, joint training improved 16-target NLL from 0.2763 to 0.2169 and accuracy from 0.9475 to 0.9535 at lambda=0.01, but gains were not consistent at larger target sizes and stronger lambda caused negative transfer.
- Agent confidence
- 7/10
- Baseline
- n=16: NLL 0.2763, accuracy 0.9475; n=64: NLL 0.0832, accuracy 0.9920; n=256: NLL 0.0230, accuracy 0.9978
- Idea
- Best joint setting lambda=0.01: n=16 NLL 0.2169, accuracy 0.9535; n=64 NLL 0.0864, accuracy 0.9933; n=256 NLL 0.0268, accuracy 0.9965. Domain accuracy was 0.87-0.90.
Limitations: Only a small synthetic two-dimensional covariate-shift benchmark was tested, with one seed and no CIFAR/MNIST comparison. FLOPs, multiple-seed confidence intervals, anchor regularization, odds clipping, and comparisons to domain-adversarial or importance-weighted baselines were not tested.
How to run: python3 experiment.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Artifacts
- experiment.py 4.3 KB View Raw
- report.md 1.3 KB View