function mndn = grt_dev(pM,oM);

[r,c,k] = size(pM);
[s,d,h] = size(oM);
if s~=r || c~=d || k~=h
    error('matrices should be the same size')
end
if sum(pM(1,:,1)) > 1
    for i = 1:k
        for j = 1:r
            pM(j,:,i) = pM(j,:,i)/sum(pM(j,:,i));
        end
    end
end
if sum(oM(1,:,1)) > 1
    for i = 1:h
        for j = 1:s
            oM(j,:,i) = oM(j,:,i)/sum(oM(j,:,i));
        end
    end
end
% hacky code depends on 4x4x5 confusion matrices
pred = [];%zeros(r*c*k,1);
obsv = [];%zeros(s*d*h,1);
for i = 1:k
    for j = 1:c
        pred = [pred; pM(:,c,k)];
        obsv = [obsv; oM(:,c,k)];
    end
end
mndn(1) = mean(abs(pred-obsv));%mean(mean(mean(abs(pM-oM))));
mndn(2) = median(abs(pred-obsv));%median(median(median(abs(pM-oM))));