/*
** Example 1.4
*/
qdata=randn(48,4); @ generating 4 variables with 48 quarters of data @
adata=qtoa(qdata); @ converting quarterly data to annual data, 4 variables @
print adata;

end;

fn qtoa1(x) = meanc(reshape(x,rows(x)/4,4)');

proc qtoa(x);
    local r,c,y,i;
    r = rows(x);
    c = cols(x);
    y = qtoa1(x[.,1]);
    i = 2;
    do until i > c;
        y = y~qtoa1(x[.,i]);
        i = i+1;
    endo;
    retp(y);
endp;