实验复现、消融、多随机种子与研究记录

import java.util.*;
public class Main {
    static void near(double expected,double actual) {
        if (!Double.isFinite(actual) || Math.abs(expected-actual)>1e-9*Math.max(1,Math.abs(expected)))
            throw new AssertionError(expected+" != "+actual);
    }
    static void check(boolean ok) { if(!ok) throw new AssertionError(); }
    static void bad(Runnable r) {
        try { r.run(); } catch(IllegalArgumentException e) { return; }
        throw new AssertionError("invalid input accepted");
    }
record Effects(double a,double b,double interaction,OptionalDouble standardError) {}
    static Effects effects(double[][] rows) {
        if(rows==null||rows.length==0) throw new IllegalArgumentException("seeds");
        double meanA=0,meanB=0,meanI=0,m2=0;int count=0;
        for(double[] r:rows) {
            if(r==null||r.length!=4) throw new IllegalArgumentException("four conditions");
            for(double x:r) if(!Double.isFinite(x)||Math.abs(x)>1e100) throw new IllegalArgumentException("score");
            double a=((r[1]-r[0])+(r[3]-r[2]))/2; // 顺序为00、10、01、11，A效应在B两种状态取平均。
            double b=((r[2]-r[0])+(r[3]-r[1]))/2;
            double interaction=(r[3]-r[2])-(r[1]-r[0]); // 差中之差，不能与正负一编码回归系数混用。
            count++;meanA+=(a-meanA)/count;meanB+=(b-meanB)/count;
            double delta=interaction-meanI;meanI+=delta/count;m2+=delta*(interaction-meanI); // 在线维护跨种子交互项样本方差。
        }
        OptionalDouble se=count<2?OptionalDouble.empty():OptionalDouble.of(Math.sqrt(Math.max(0,m2)/(count-1)/count));
        return new Effects(meanA,meanB,meanI,se);
    }
    public static void main(String[] args) {
Effects a=effects(new double[][]{{10,12,13,15}});near(2,a.a());near(3,a.b());near(0,a.interaction());check(a.standardError().isEmpty());
        Effects b=effects(new double[][]{{0,0,0,4}});near(2,b.a());near(2,b.b());near(4,b.interaction());
        Effects c=effects(new double[][]{{10,12,13,17},{20,22,23,27}});near(2,c.interaction());near(0,c.standardError().orElseThrow());
        Effects d=effects(new double[][]{{0,0,0,0},{0,0,0,2}});near(1,d.interaction());near(1,d.standardError().orElseThrow());
        Effects e=effects(new double[][]{{4,3,2,1}});near(-1,e.a());near(-2,e.b());
        bad(()->effects(new double[][]{}));
        bad(()->effects(new double[][]{{1,2,3}}));
        bad(()->effects(new double[][]{{1,2,3,Double.NaN}}));
        System.out.println("PASS 8 cases");
    }
}