comp perf

comp_perf.jl

t1file = "/tmp/after.txt"
t2file = "/tmp/before.txt"
NS = 5
TP = 2.31

data1 = readdlm(t1file);
data2 = readdlm(t2file);

testnames = data1[1:end,1];
mean1 = data1[1:end,3];
mean2 = data2[1:end,3];
sd1 = data1[1:end,4];
sd2 = data2[1:end,4];

function tvalue(mean1, mean2, sd1, sd2, NS)
    N = length(mean1)
    t = Array(Float64, N)
    for i in 1:N
        num = mean1[i] - mean2[i]
        den = sqrt(0.5*(sd1[i]^2 + sd2[i]*2)) * sqrt(2/NS)
        t[i] = num/den
    end
    [abs(x) for x in t]
end

tv = tvalue(mean1, mean2, sd1, sd2, 5);
tvf = [tv .> TP];
println("failed: $(testnames[tvf])")