Return to rogers.exec CVS log

Up to [Research Unix] / researchv10dc / cmd / matlab / exec

Annotation of researchv10dc/cmd/matlab/exec/rogers.exec, revision 1.1

1.1     ! root        1: exec('d.boug');                        // reads data
        !             2: <g,k> = size(p);               // p is matrix of gene frequencies
        !             3: wv = ncen/sum(ncen);           // ncen contains population sizes
        !             4: pbar = wv*p;                   // weighted average of p
        !             5: p = p - ones(g,1)*pbar;        // deviations from mean
        !             6: p = sqrt(diag(wv)) * p;        // weight rows of p by sqrt of pop size
        !             7: h = diag(pbar); h = h*(eye-h); // diagonal contains binomial variance: p*(1-p)
        !             8: r = p*inv(h)*p'/k;             // normalized covariance matrix
        !             9: eig(r)'

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.