A Julia function to calculate area under the disease progress curve (AUDPC)

audpc.jl

function audpc(evaluation, dates)
    n = length(evaluation) - 1
    disvec = Base.zeros(n)
    datevec = Base.zeros(n)
    out = 0.0

    for i in 1:n
        disvec[i] = (evaluation[i] + evaluation[i + 1]) / 2
        datevec[i] = dates[i + 1] - dates[i]
        out = sum(disvec .* datevec)
    end

    return out

end

Example from the audpc() help from R's agricolae package.

dates = [14, 21, 28] # days
evaluation = [40, 80, 90]
audpc(evaluation, dates)

1015.0