calculateMaxDD_simple.m

function [maxDD maxDDD]=calculateMaxDD_simple(cumret)
% [maxDD maxDDD]=calculateMaxDD_simple(cumret)
% calculation of maximum drawdown and maximum drawdown duration based on
% cumulative simple (uncompounded) returns.

highwatermark=zeros(size(cumret)); % initialize high watermarks to zero.

drawdown=zeros(size(cumret)); % initialize drawdowns to zero.

drawdownduration=zeros(size(cumret)); % initialize drawdown duration to zero.

for t=2:length(cumret)
    highwatermark(t)=max(highwatermark(t-1), cumret(t));
    drawdown(t)=cumret(t)-highwatermark(t); % drawdown on each day
    if (drawdown(t)==0)
        drawdownduration(t)=0;
    else
        drawdownduration(t)=drawdownduration(t-1)+1;
    end
end

maxDD=min(drawdown); % maximum drawdown

maxDDD=max(drawdownduration); % maximum drawdown duration